Excited states and electronic spectra of annulated dinuclear free-base phthalocyanines: a theoretical study on near-infrared-absorbing dyes

The electronic excited states and electronic absorption spectra of annulated dinuclear free-base phthalocyanine (C(58)H(30)N(16)) are studied through quantum chemical calculations using the symmetry-adapted cluster-configuration interaction (SAC-CI) method. Three tautomers are possible with respect to the position of the pyrrole protons; therefore, the SAC-CI calculations for these tautomers were performed. The structures of the Q-band states are discussed based on the character of their molecular orbitals. The lower energy shift of the Q-bands because of dimerization is explained by the decrease in the HOMO-LUMO gaps resulting from the bonding and antibonding interactions between the monomer units. The electronic dipole moments of the nonsymmetric tautomer were calculated, and the possibility of charge-separated excited states is discussed. The relative energies of these tautomers are examined using density functional theory (DFT) calculations for several peripheral substituents. The relative energies of these tautomers significantly depend on the substituents, and therefore, the abundance ratios of the three tautomers were affected by the substituents. The absorption spectra were simulated from the SAC-CI results weighted by the Boltzmann factors obtained from the DFT calculations. The SAC-CI spectra reproduce the experimental findings well. The thermal-averaged SAC-CI spectra could explain the observed substituent effect on the structure of the Q-bands in terms of the relative stabilities and the abundance ratios of the tautomers. The SAC-CI and time-dependent density functional theory calculations are also compared. The CAM-B3LYP results agreed with the trends of the SAC-CI results; however, the CAM-B3LYP calculation overestimated the excitation energies in comparison with the SAC-CI and experimental results.     

Singly and doubly excited states of butadiene, acrolein, and glyoxal: Geometries and electronic spectra

Excited-state geometries and electronic spectra of butadiene, acrolein, and glyoxal have been investigated by the symmetry adapted cluster configuration interaction (SAC-CI) method in their s-trans conformation. Valence and Rydberg states below the ionization threshold have been precisely calculated with sufficiently flexible basis sets. Vertical and adiabatic excitation energies were well reproduced and the detailed assignments were given taking account of the second moments. The deviations of the vertical excitation energies from the experiment were less than 0.3 eV for all cases. The SAC-CI geometry optimization has been applied to some valence and Rydberg excited states of these molecules in the planar structure. The optimized ground- and excited-state geometries agree well with the available experimental values; deviations lie within 0.03 A and 0.7 degrees for the bond lengths and angles, respectively. The force acting on the nuclei caused by the excitations has been discussed in detail by calculating the SAC-CI electron density difference between the ground and excited states; the geometry relaxation was well interpreted with the electrostatic force theory. In Rydberg excitations, geometry changes were also noticed. Doubly excited states (so-called 2 (1)A(g) states) were investigated by the SAC-CI general-R method considering up to quadruple excitations. The characteristic geometrical changes and large energetic relaxations were predicted for these states.     

Electronic structure and optical properties of chelating heteroatomic conjugated molecules: a SAC-CI study

The electronic structure and optical properties of 13 chelating heteroatomic conjugated molecules such as pyridine, benzoxazole, and benzothiazole derivatives, which are used as C–N ligands in organometallic compounds, have been investigated. The geometries of the ground and first excited states were obtained by the DFT and CIS methods, respectively, followed by the SAC-CI calculations of the transition energies for absorption and emission. For six compounds whose experimental data are available, the SAC-CI calculations reproduced the experimental values satisfactorily with deviations of less than 0.3 eV for absorption and 0.1 eV for emission except for benzoxazoles. For other molecules, the theoretical absorption and emission spectra were predicted. The lowest ππ* excited-state geometries was calculated to be planar for most of the molecules with two or three conjugated rings connected by single bond. The geometry change due to the ππ* excitation was qualitatively interpreted by electrostatic force theory based on SAC/SAC-CI electron density difference. The excitations are relatively localized in the central region and in the lowest ππ* excited state, the inter-ring single bond shows large change, with a contraction of 0.05–0.09 Å. The present calculations provide reliable information regarding the energy levels of these chelating heteroatomic conjugated compounds.     

Excited states of porphyrin isomers and porphycene derivatives: a SAC-CI study

Excited states of free-base porphyrin isomers, porphycene (Pc), corrphycene (Cor), and hemiporphycene (hPc), were studied by the Symmetry-Adapted Cluster (SAC)/SAC-Configuration Interaction (CI) method. The absorption peaks of the porphyrin isomers were assigned on the basis of the SAC-CI spectra. The X, Y, X', and Y' bands of the porphyrin isomers, which have weak intensities, are identified. The differences in the Q-band absorptions among the isomers were clearly explained by the four-orbital model. In Cor and hPc, the wave function of the B-band corresponds to the mixture of the four-orbital excitations and the optically forbidden excitation of free-base porphin (P), due to the molecular symmetry lowering in the isomers. The B-band character is described by the five-orbital model in Pc and the six-orbital model in Cor and hPc. Two tetrazaporphycenes and two ring-extended (dibenzo) porphycenes were designed, and the Q-band transition moment was successfully controlled. These examples show that the control of the four-orbital energy levels is the guiding principle for pigment design in porphyrin compounds.     

Investigation of the electronic spectra and excited-state geometries of poly(para-phenylene vinylene) (PPV) and poly(para-phenylene) (PP) by the symmetry-adapted cluster configuration interaction (SAC-CI) method

The symmetry-adapted cluster-configuration interaction (SAC-CI) method has been used to investigate the optical and geometric properties of the oligomers of poly(para-phenylene vinylene) (PPV) and poly(para-phenylene) (PP). Vertical singlet and triplet absorption spectra and emission spectra have been calculated accurately; the mean average deviation from available experimental results lies within 0.2 eV. The chain length dependence of the transition energies has been improved in comparison to earlier TDDFT and MRSDCI calculations. The present analysis suggests that conventional TDDFT with the B3LYP functional should be used carefully, as it can provide inaccurate estimates of the chain length dependence of the excitation energies of these molecules with long pi conjugation. The T1 state was predicted to be at a lower energy, by 1.0-1.5 eV for PPV and by 0.9-1.7 eV for PP, than the S1 state, which indicates a localized T1 state with large exchange energy. By calculating the SAC-CI electron density difference between the ground and excited states, the geometry relaxations due to excitations can be analyzed in detail using electrostatic force theory. For trans-stilbene, the doubly excited 21Ag state was studied, and the calculated transition energy of 4.99 eV agrees very well with the experimental value of 4.84 eV. In contrast to previous ab initio calculations, we predict this doubly excited 21Ag state to lie above the 11Bu state.     

The fate of optical excitations in small hydrated ZnS clusters: a theoretical study into the effect of hydration on the excitation and localisation of electrons in Zn4S4 and Zn6S6

In this paper we explore the effect of water on the excited state properties of ZnS nanostructures by means of time-dependent density functional theory (TD-DFT) calculations. Using these TD-DFT calculations we show that the effect of water on the optical absorption spectra is primarily a small blue-shift and that a secondary effect is that spectroscopic features that correspond to dark excitations for the anhydrous nanostructures gain intensity and new absorption peaks are predicted to appear. The effect of adsorbed water on the localisation of excited states is to produce small shifts in the values of the excited stabilisation energies but, more importantly, it results in the formation of extra minima when compared with the case for anhydrous ZnS. Finally, the effect of water on photoluminescence (PL) energies is predicted to be small but the appearance of extra minima induced by the presence of adsorbed water is expected to lead to a splitting/broadening of the PL signal.     

A quantitative prediction of the electronic spectra of thiocarbonyl chromophores: TD-DFT versus SAC-CI

In this contribution, we set up a SAC-CI methodology to evaluate the n → π* and π → π* vertical transition energies of a series of thiocarbonyl derivatives. We show that Frozen-Core SAC-CI provides accurate vertical excitations energies. Nevertheless, in order to obtain converged results, the R2S2 unliked integrals have to be taken into account in L3-SAC-CI calculations. In addition, we present the comparative performances of three computational procedures, INDO/S, TD-DFT, and SAC-CI, for the calculation of valence excited states energies and it turns out that: (1) no tuning of the exact exchange (α) included in TD-DFT allows to consistently reproduce the SAC-CI results; (2) SAC-CI and TD-B3LYP both evaluate the n → π*, as well as the π → π* transition energies, with a similar accuracy. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Effects of chain length and Au spin-orbit coupling on 3(pi pi*) emission from bridging Cn2- units: theoretical characterization of spin-forbidden radiative transitions in metal-capped one-dimensional carbon chains [H3PAu(C[triple bond]C)nAuPH3

Density functional theory and CASSCF calculations have been used to optimize the geometries of binuclear gold(I) complexes [H(3)PAu(C[triple bond]C)(n)AuPH(3)] (n=1-6) in their ground states and selected lowest energy (3)(pi pi*) excited states. Vertical excitation energies obtained by time-dependent density functional calculations for the spin-forbidden singlet-triplet transitions have exponential-decay size dependence. The predicted singlet-triplet splitting limit of [H(3)PAu(C[triple bond]C)(proportional/variant)AuPH(3)] is about 8317 cm(-1). Calculated singlet-triplet transition energies are in reasonable agreement with available experimental observations. The effect of the heavy atom Au spin-orbit coupling on the (3)(pi pi*) emission of these metal-capped one-dimensional carbon allotropes has been investigated by MRCI calculations. The contribution of the spin- and dipole-allowed singlet excited state to the spin-orbit-coupling wave function of the (3)(pi pi*) excited state makes the low-lying acetylenic triplet excited states become sufficiently allowed so as to appear in both electronic absorption and emission.     

Formulation and implementation of direct algorithm for the symmetry-adapted cluster and symmetry-adapted cluster-configuration interaction method

We present a new computational algorithm, called direct algorithm, for the symmetry-adapted cluster (SAC) and SAC-configuration interaction (SAC-CI) methodology for the ground, excited, ionized, and electron-attached states. The perturbation-selection technique and the molecular orbital index based direct sigma-vector algorithm were combined efficiently with the use of the sparse nature of the matrices involved. The formal computational cost was reduced to O(N(2)xM) for a system with N-active orbitals and M-selected excitation operators. The new direct SAC-CI program has been applied to several small molecules and free-base porphin and has been shown to be more efficient than the conventional nondirect SAC-CI program for almost all cases. Particularly, the acceleration was significant for large dimensional computations. The direct SAC-CI algorithm has achieved an improvement in both accuracy and efficiency. It would open a new possibility in the SAC/SAC-CI methodology for studying various kinds of ground, excited, and ionized states of molecules.     

Symmetry-adapted-cluster/symmetry-adapted-cluster configuration interaction methodology extended to giant molecular systems: ring molecular crystals

The symmetry adapted cluster (SAC)/symmetry adapted cluster configuration interaction (SAC-CI) methodology for the ground, excited, ionized, and electron-attached states of molecules was extended to giant molecular systems. The size extensivity of energy and the size intensivity of excitation energy are very important for doing quantitative chemical studies of giant molecular systems and are designed to be satisfied in the present giant SAC/SAC-CI method. The first extension was made to giant molecular crystals composed of the same molecular species. The reference wave function was defined by introducing monomer-localized canonical molecular orbitals (ml-CMO's), which were obtained from the Hartree-Fock orbitals of a tetramer or a larger oligomer within the electrostatic field of the other part of the crystal. In the SAC/SAC-CI calculations, all the necessary integrals were obtained after the integral transformation with the ml-CMO's of the neighboring dimer. Only singles and doubles excitations within each neighboring dimer were considered as linked operators, and perturbation selection was done to choose only important operators. Almost all the important unlinked terms generated from the selected linked operators were included: the unlinked terms are important for keeping size extensivity and size intensivity. Some test calculations were carried out for the ring crystals of up to 10 000-mer, confirming the size extensivity and size intensivity of the calculated results and the efficiency of the giant method in comparison with the standard method available in GAUSSIAN 03. Then, the method was applied to the ring crystals of ethylene and water 50-mers, and formaldehyde 50-, 100-, and 500-mers. The potential energy curves of the ground state and the polarization and electron-transfer-type excited states were calculated for the intermonomer distances of 2.8-100 A. Several interesting behaviors were reported, showing the potentiality of the present giant SAC/SAC-CI method for molecular engineering.     

The excited states of cycloporane. MCD spectrum,and CD spectrum of an optically active derivative

The MCD spectrum of gas phase cyclopropane and the CD spectrum of trans 1,2-dimethyl cyclopropane were measured in the spectral region 210–140 nm. The absorption spectra are also reported. Comparison of CD and absorption spectra of dimethyl cyclopropane and consideration of the anisotropy factors g = Δε/ε as a guide to the assignment of magnetic dipole allowed transitions led to the ordering of the states of the first excited configuration as A′₁, A′₂, E′ in order of increasing energy. The magnetic moments of the two observed allowed excited states of cyclopropane were determined from the MCD spectrum. LCAO MO level calculations of the MCD parameters $\text{A}$₁/$\text{D}$₀ for the lowest three excited E′ states were carried out and the results were used to discuss the assignment of the allowed transitions. It is concluded that configuration interaction is of considerable importance in the low energy excited states of cyclopropane.     

Magnetic circular dichroism in real-time time-dependent density functional theory

We apply the adiabatic time-dependent density functional theory to magnetic circular dichroism (MCD) spectra using the real-space, real-time computational method. The standard formulas for the MCD response and its A and B terms are derived from the observables in the time-dependent wave function. We find real-time method is well suited for calculating the overall spectrum, particularly at higher excitation energies where individual excited states are numerous and overlapping. The MCD sum rules are derived and intepreted in the real-time formalism; we find that they are very useful for normalization purposes and assessing the accuracy of the theory. The method is applied to MCD spectrum of C(60) using the adiabatic energy functional from the local density approximation. The theory correctly predicts the signs of the A and B terms for the lowest allowed excitations. However, the magnitudes of the terms only show qualitative agreement with experiment.     

Electronic and vibronic contributions to two-photon absorption in donor-acceptor-donor squaraine chromophores

Many squaraines have been observed to exhibit two-photon absorption at transition energies close to those of the lowest energy one-photon electronic transitions. Here, the electronic and vibronic contributions to these low-energy two-photon absorptions are elucidated by performing correlated quantum-chemical calculations on model chromophores that differ in their terminal donor groups (diarylaminothienyl, indolenylidenemethyl, dimethylaminopolyenyl, or 4-(dimethylamino)phenylpolyenyl). For squaraines with diarylaminothienyl and dimethylaminopolyenyl donors and for the longer examples of 4-(dimethylamino)phenylpolyenyl donors, the calculated energies of the lowest two-photon active states approach those of the lowest energy one-photon active (1B(u)) states. This is consistent with the existence of purely electronic channels for low-energy two-photon absorption (TPA) in these types of chromophores. On the other hand, for all squaraines containing indolinylidenemethyl donors, the calculations indicate that there are no low-lying electronic states of appropriate symmetry for TPA. Actually, we find that the lowest energy TPA transitions can be explained through coupling of the one-photon absorption (OPA) active 1B(u) state with b(u) vibrational modes. Through implementation of Herzberg-Teller theory, we are able to identify the vibrational modes responsible for the low-energy TPA peak and to reproduce, at least qualitatively, the experimental TPA spectra of several squaraines of this type.     

Theoretical prediction of valence and Rydberg excited states: Minnesota exchange-correlation functionals vs symmetry adapted cluster-configuration interaction

During this contribution, we present a benchmark investigation on the applicability of several Minnesota functionals from various classes like local meta-generalized and meta-nonseparable gradient approximations, hybrids, and range-separated hybrids for describing the valence and Rydberg excitation energies of some organic compounds from different categories. Furthermore, the performances of Minnesota density functionals from density functional theory are also assessed against a wave function theory based approach in the context of excite states calculations, symmetry adapted cluster-configuration interaction (SAC-CI) method. Pragmatically, the singles and doubles linked excitation operators are considered in the SAC-CI wave functions. With more or less different accountabilities of the considered methods, it is shown that the M06-2X, M05-2X, and M11 functionals have the best performances for valence excited states. On the other hand, for Rydberg excited states although the SAC-CI method outperforms others, the statistical analyses reveal that the efficiency of some Minnesota functionals is also respectable.     

Ab initio CI calculations on free-base porphin

Ab initio configuration interaction (CI ) calculations were carried out on low-lying singlet and triplet π–π* states and ionized states of free-base porphin. We take into account single and double excitations from σ and π electrons in the CI calculations. The composite natural orbitals were employed in order to reduce the size of orbital set to be used in the CI . The calculated excitation energies were in good agreement with experimental values. The use of split-valence-type basis and the inclusion of correlation effects of σ electrons were proved to be important to describe the low-lying π–π* states, especially the Soret band.     

Electronic transitions in [Re6S8X6](4-) (X = Cl, Br, I): results from time-dependent density functional theory and solid-state calculations

Relativistic time-dependent density functional theory (TDDFT) calculations were performed on the excited states of the [Re6S8X6](4-) (X = Cl, Br, I) series. For all members of the series, the lowest excited states in the spectra do not correspond to a ligand-to-metal (or ligand-to-cluster) excitation but rather a cluster-cluster transition from the HOMO e(g) to antibonding t(1u) orbitals with only a modest admixture of Re-X sigma* character. These results lead to a re-evaluation of the role of the axial ligand in these compounds. The calculated excitation energies reproduce the experimental absorption and emission spectra. This work also confirms previous TDDFT calculations on the emission energies. Results for discrete cluster ions are compared with those obtained from calculations in the solid state in Cs4[Re6S8X6].CsX (X = Cl, Br) and Cs4[Re6S8I6].2CsI. Significant differences are seen in the relatively higher energies of the antibonding t(1u) orbital in the solid-state case, and an inversion in the orbital character of the two allowed absorptions is calculated. The e(g) (HOMO)-to-a(2g) (LUMO) orbital energy differences corresponding to the emission transition are quite comparable for the solid state and discrete cluster calculations, and both overestimate the observed emission energy by the same margin.     

SAC-CI study of the excited states of free base tetrazaporphin

The SAC (symmetry adapted cluster)/SAC-CI method is applied to calculations of the ground and excited states of free base tetrazaporphin (FBTAP). The electronic spectrum is reproduced in fairly good agreement with experiments, and unknown absorption bands in the energy region higher than 4 eV are predicted. The effect of meso-tetraaza-substitution on the excited states of porphin macrocycle explains why the compound is colored, as in the closely related compounds, phthalocyanines, used as pigments.