Probing the effects of one-electron reduction and protonation on the electronic properties of the Fe-S clusters in the active-ready form of [FeFe]-hydrogenases. A QM/MM investigation

A QM/MM investigation of the active-ready (H(ox)) form of [FeFe]-hydrogenase from D. desulfuricans, in which the electronic properties of all Fe-S clusters (H, F and F') have been simultaneously described using DFT, was carried out with the aim of disclosing a possible interplay between the H-cluster and the accessory iron-sulfur clusters in the initial steps of the catalytic process leading to H(2) formation. It turned out that one-electron addition to the active-ready form leads to reduction of the F'-cluster and not of the H-cluster. Protonation of the H-cluster in H(ox) is unlikely, and in any case it would not trigger electron transfer from the accessory Fe(4)S(4) clusters to the active site. Instead, one-electron reduction and protonation of the active-ready form trigger electron transfer within the protein, a key event in the catalytic cycle. In particular, protonation of the H-cluster after one-electron reduction of the enzyme lowers the energy of the lowest unoccupied molecular orbitals localized on the H-cluster to such an extent that a long-range electron transfer from the F'-cluster towards the H-cluster itself is allowed.     

The structure of the hydrated electron. Part 1. Magnetic resonance of internally trapping water anions: a density functional theory study

Density functional theory is used to rationalize magnetic parameters of hydrated electron trapped in alkaline glasses as observed using electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectroscopies. To this end, model water cluster anions (n=4-8 and n=20, 24) that localize the electron internally are examined. It is shown that hyperfine coupling tensors of H/D nuclei in the water molecules are defined mainly by the cavity size and the coordination number of the electron; the water molecules in the second solvation shell play a relatively minor role. An idealized model of the hydrated electron (that is usually attributed to L. Kevan) in which six hydroxyl groups arranged in an octahedral pattern point toward the common center is shown to provide the closest match to the experimental parameters, such as isotropic and anisotropic hyperfine coupling constants for the protons (estimated from ESEEM), the second moment of the EPR spectra, and the radius of gyration. The salient feature is the significant transfer (10-20%) of spin density into the frontal O 2p orbitals of water molecules. Spin bond polarization involving these oxygen orbitals accounts for small, negative hyperfine coupling constants for protons in hydroxyl groups that form the electron-trapping cavity. In Part 2, these results are generalized for more realistic geometries of core anions obtained using a dynamic one-electron mixed quantum/classical molecular dynamics model.     

Exploring the thermochromism of sulfite-embedded polyoxometalate capsules

The Dawson-type polyanion [α-Mo(18)O(54)(SO(3))(2)](4-), with two SO(3)(2-) templates embedded inside a polyoxomolybdate(vi) cage, exhibits thermochromism over an exceptionally wide temperature range (～500 K). The temperature dependence of the cluster structure, established from X-ray crystallography, IR and Raman spectroscopy and DFT calculations, is related to a decreasing HOMO-LUMO gap in the near UV with increasing temperature. We postulate this is due to geometrical changes that affect both the occupied and unoccupied frontier molecular orbitals of this cluster anion.     

Theoretical study on structural, electronic, and optical properties of ambipolar diphenylamino end-capped oligofluorenylthiophenes and fluoroarene-thiophene as light-emitting materials

Ambipolar diphenylamino end-capped oligofluorenylthiophenes and fluoroarene-thiophene show great potential for application in organic light-emitting diodes (OLEDs). Here, we provide an in-depth investigation on the optical and electronic properties of OF(2)TP-NPh ( 1a), OF(2)DTP-NPh ( 2a), OF(2)TTP-NPh ( 3a), OF(2)QTP-NPh ( 4a), and 2,5-bis-(2,3,5,6-tetrafluoro-4-trifluoromethyl-phenyl)-2,2':5',2':5',2'-quaterthiophene ( 5a). The geometric and electronic structures of the oligomers in the ground-state are studied with density functional theory (DFT) and ab initio Hartree-Fock, whereas the lowest singlet excited states are optimized by ab initio CIS. The energies of the lowest singlet excited states are calculated by employing time-dependent density functional theory (TDDFT). The results show that the highest occupied molecular orbitals, lowest unoccupied molecular orbitals, energy gaps, ionization potentials, and electron affinities for the oligomers are affected by the thiophene chain length and the different end-caps. The absorption and emission spectra exhibit red shifts to some extent due to the increasing thiophene chain length and the enhancing electron-donating property of the end-caps. Furthermore, the large Stokes shifts ranging from 58 to 80 nm are examined, resulting from a more planar conformation of the excited-state between the two adjacent units in the oligomers. All the calculated data show that the fluoroarene-thiophene has improved electron transport rate and charge transfer balance performance, and all the studied molecules can be used as ambipolar-transporting materials in OLEDs.     

Generalized hybrid-orbital method for combining density functional theory with molecular mechanicals.

The generalized hybrid orbital (GHO) method has previously been formulated for combining molecular mechanics with various levels of quantum mechanics, in particular semiempirical neglect of diatomic differential overlap theory, ab initio Hartree-Fock theory, and self-consistent charge density functional tight-binding theory. To include electron-correlation effects accurately and efficiently in GHO calculations, we extend the GHO method to density functional theory in the generalized-gradient approximation and hybrid density functional theory (denoted by GHO-DFT and GHO-HDFT, respectively) using Gaussian-type orbitals as basis functions. In the proposed GHO-(H)DFT formalism, charge densities in auxiliary hybrid orbitals are included to calculate the total electron density. The orthonormality constraints involving the auxiliary Kohn-Sham orbitals are satisfied by carrying out the hybridization in terms of a set of L?wdin symmetrically orthogonalized atomic basis functions. Analytical gradients are formulated for GHO-(H)DFT by incorporating additional forces associated with GHO basis transformations. Scaling parameters are introduced for some of the one-electron integrals and are optimized to obtain the correct charges and geometry near the QM/MM boundary region. The GHO-(H)DFT method based on the generalized gradient approach (GGA) (BLYP and mPWPW91) and HDFT methods (B3 LYP, mPW1PW91, and MPW1 K) is tested-for geometries and atomic charges-against a set of small molecules. The following quantities are tested: 1) the C--C stretch potential in ethane, 2) the torsional barrier for internal rotation around the central C--C bond in n-butane, 3) proton affinities for a set of alcohols, amines, thiols, and acids, 4) the conformational energies of alanine dipeptide, and 5) the barrier height of the hydrogen-atom transfer between n-C4H10 and n-C4H9, where the reaction center is described at the MPW1 K/6-31G(d) level of theory.     

Nonadiabatic molecular dynamics simulations of correlated electrons in solution. 2. A prediction for the observation of hydrated dielectrons with pump-probe spectroscopy

The hydrated dielectron is a highly correlated, two-electron, solvent-supported state consisting of two spin-paired electrons confined to a single cavity in liquid water. Although dielectrons have been predicted to exist theoretically and have been used to explain the lack of ionic strength effect in the bimolecular reaction kinetics of hydrated electrons, they have not yet been observed directly. In this paper, we use the extensive nonadiabatic mixed quantum/classical excited-state molecular dynamics simulations from the previous paper to calculate the transient spectroscopy of hydrated dielectrons. Because our simulations use full configuration interaction (CI) to determine the ground and excited state two-electron wave functions at every instant, our nonequilibrium simulations allow us to compute the absorption, stimulated emission (SE), and bleach spectroscopic signals of both singlet and triplet dielectrons following excitation by ultraviolet light. Excited singlet dielectrons are predicted to display strong SE in the mid infrared and a transient absorption in the near-infrared. The near-infrared transient absorption of the singlet dielectron, which occurs near the peak of the (single) hydrated electron's equilibrium absorption, arises because the two electrons tend to separate in the excited state. In contrast, excitation of the hydrated electron gives a bleach signal in this wavelength region. Thus, our calculations suggest a clear pump-probe spectroscopic signature that may be used in the laboratory to distinguish hydrated singlet dielectrons from hydrated electrons: By choosing an excitation energy that is to the blue of the peak of the hydrated electron's absorption spectrum and probing near the maximum of the single electron's absorption, the single electron's transient bleach signal should shrink or even turn into a net absorption as sample conditions are varied to produce more dielectrons.     

Exploring the role of decoherence in condensed-phase nonadiabatic dynamics: a comparison of different mixed quantum/classical simulation algorithms for the excited hydrated electron

Mixed quantum/classical (MQC) molecular dynamics simulation has become the method of choice for simulating the dynamics of quantum mechanical objects that interact with condensed-phase systems. There are many MQC algorithms available, however, and in cases where nonadiabatic coupling is important, different algorithms may lead to different results. Thus, it has been difficult to reach definitive conclusions about relaxation dynamics using nonadiabatic MQC methods because one is never certain whether any given algorithm includes enough of the necessary physics. In this paper, we explore the physics underlying different nonadiabatic MQC algorithms by comparing and contrasting the excited-state relaxation dynamics of the prototypical condensed-phase MQC system, the hydrated electron, calculated using different algorithms, including: fewest-switches surface hopping, stationary-phase surface hopping, and mean-field dynamics with surface hopping. We also describe in detail how a new nonadiabatic algorithm, mean-field dynamics with stochastic decoherence (MF-SD), is to be implemented for condensed-phase problems, and we apply MF-SD to the excited-state relaxation of the hydrated electron. Our discussion emphasizes the different ways quantum decoherence is treated in each algorithm and the resulting implications for hydrated-electron relaxation dynamics. We find that for three MQC methods that use Tully's fewest-switches criterion to determine surface hopping probabilities, the excited-state lifetime of the electron is the same. Moreover, the nonequilibrium solvent response function of the excited hydrated electron is the same with all of the nonadiabatic MQC algorithms discussed here, so that all of the algorithms would produce similar agreement with experiment. Despite the identical solvent response predicted by each MQC algorithm, we find that MF-SD allows much more mixing of multiple basis states into the quantum wave function than do other methods. This leads to an excited-state lifetime that is longer with MF-SD than with any method that incorporates nonadiabatic effects with the fewest-switches surface hopping criterion.     

Electronic, optical, and charge transfer properties of donor–bridge–acceptor hydrazone sensitizers

The ground state geometries have been computed by using density functional theory (DFT) at B3LYP/6-31G*, B3LYP/6-31G**, and PCM-B3LYP/6-31G* level of theories. The highest occupied molecular orbitals (HOMOs) are delocalized on whole of the molecule and the lowest unoccupied molecular orbitals (LUMOs) are localized on the tricarbonitrile. The lowest HOMO and LUMO energies have been observed for Dye1 while highest for Dye4. The LUMO energies of Dye1–Dye4 are above the conduction band of TiO₂ and HOMOs are below the redox couple. The absorption spectra have been computed in solvent (methanol) and without solvent by using time-dependant DFT at TD-B3LYP/6-31G*, TD-B3LYP/6-31G**, and PCM-TD-B3LYP/6-31G* level of theories. The calculated maximum absorption wavelengths of the spectra in methanol are in good agreement with experimental evidences. The maximum absorption wavelengths of new designed sensitizers are red shifted compared to parent molecule. The electronic coupling constant and electron injection have been computed by first principle investigations. The improved electronic coupling constant and electron injection revealed that new modeled systems would be efficient sensitizers.     

DFT performance for the hole transfer parameters in DNA π stacks

Recently, we showed that unoccupied Kohn‐Sham (KS) orbitals stemming from DFT calculations of a neutral system can be used to derive accurate estimates of the free energy and electronic couplings for excess electron transfer in DNA (Félix and Voityuk, J Phys Chem A 2008, 112, 9043). In this article, we consider the propagation of radical cation states (hole transfer) through DNA π‐stacks and compare the performance of different exchange‐correlation functionals to estimate the hole transfer (HT) parameters. Two different approaches are used: (1) calculations that use occupied KS orbitals of neutral π stacks of nucleobases, and (2) the time‐dependent DFT method which is applied to the radical cation states of these stacks. Comparison of the calculated parameters with the reference data suggests that the best results are provided by the KS scheme with hybrid functionals (B3LYP, PBE0, and BH&HLYP). The TD DFT approach gives significantly less accurate values of the HT parameters. In agreement with high‐level ab initio results, the KS scheme predicts that the hole in π stacks is confined to a single nucleobase; in contrast, the spin‐unrestricted DFT method considerably overestimates the hole delocalization in the radical cations. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Electronic structure of copper phthalocyanine: an experimental and theoretical study of occupied and unoccupied levels

An experimental and theoretical study of the electronic structure of copper phthalocyanine (CuPc) molecule is presented. We performed x-ray photoemission spectroscopy (XPS) and photoabsorption [x-ray absorption near-edge structure (XANES)] gas phase experiments and we compared the results with self-consistent field, density functional theory (DFT), and static-exchange theoretical calculations. In addition, ultraviolet photoelectron spectra (UPS) allowed disentangling several outer molecular orbitals. A detailed study of the two highest occupied orbitals (having a(1u) and b(1g) symmetries) is presented: the high energy resolution available for UPS measurements allowed resolving an extra feature assigned to vibrational stretching in the pyrrole rings. This observation, together with the computed DFT electron density distributions of the outer valence orbitals, suggests that the a(1u) orbital (the highest occupied molecular orbital) is mainly localized on the carbon atoms of pyrrole rings and it is doubly occupied, while the b(1g) orbital, singly occupied, is mainly localized on the Cu atom. Ab initio calculations of XPS and XANES spectra at carbon K edge of CuPc are also presented. The comparison between experiment and theory revealed that, in spite of being formally not equivalent, carbon atoms of the benzene rings experience a similar electronic environment. Carbon K-edge absorption spectra were interpreted in terms of different contributions coming from chemically shifted C 1s orbitals of the nonequivalent carbon atoms on the inner ring of the molecule formed by the sequence of CN bonds and on the benzene rings, respectively, and also in terms of different electronic distributions of the excited lowest unoccupied molecular orbital (LUMO) and LUMO+1. In particular, the degenerate LUMO appears to be mostly localized on the inner pyrrole ring.     

Syntheses, X-ray structures, photochemistry, redox properties, and DFT calculations of interconvertible fac- and mer-[Mn(SPS)(CO)3] isomers containing a flexible SPS-based pincer ligand

The lithium salt of the anionic SPS pincer ligand composed of a central hypervalent lambda4-phosphinine ring bearing two ortho-positioned diphenylphosphine sulfide side arms reacts with [Mn(CO)5Br] to give fac-[Mn(SPS)(CO)3]. This isomer can be converted photochemically to mer-[Mn(SPS)(CO)3], with a very high quantum yield (0.80+/-0.05). The thermal backreaction is slow (taking ca. 8 h at room temperature), in contrast to rapid electrode-catalyzed mer-to-fac isomerization triggered by electrochemical reduction of mer-[Mn(SPS)(CO)3]. Both geometric isomers of [Mn(SPS)(CO)3] have been characterized by X-ray crystallography. Both isomers show luminescence from a low-lying 3IL (SPS-based) excited state. The light emission of fac-[Mn(SPS)(CO)3] is largely quenched by the efficient photoisomerization occurring probably from a low-lying Mn-CO dissociative excited state. Density functional theory (DFT) and time-dependent DFT calculations describe the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of fac- and mer-[Mn(CO)3(SPS)] as ligand-centered orbitals, largely localized on the phosphinine ring of the SPS pincer ligand. In line with the ligand nature of its frontier orbitals, fac-[Mn(SPS)(CO)3] is electrochemically reversibly oxidized and reduced to the corresponding radical cation and anion, respectively. The spectroscopic (electron paramagnetic resonance, IR, and UV-vis) characterization of the radical species provides other evidence for the localization of the redox steps on the SPS ligand. The smaller HOMO-LUMO energy difference in the case of mer-[Mn(CO)3(SPS)], reflected in the electronic absorption and emission spectra, corresponds with its lower oxidation potential compared to that of the fac isomer. The thermodynamic instability of mer-[Mn(CO)3(SPS)], confirmed by the DFT calculations, increases upon one-electron reduction and oxidation of the complex.     

A simple extension of the external magnasco-perico localization procedure to the virtual MO-Space

The external localization procedure of Magnasco and Perico is extended to the unoccupied molecular orbitals of the Fock-operator. The formal correspondence between bonding orbitals and localized antibonding MOs is demonstrated. Localized occupied and virtual one-electron functions are calculated within a semiempirical INDO-Hamiltonian and are analyzed; the externally localized occupied MOs are compared with energy localized orbitals computed by the Edmiston and Ruedenberg procedure. Various applications of the fully localized (occupied and virtual) MO set are discussed.     

Theoretical characterization of four distinct isomer types in hydrated-electron clusters, and proposed assignments for photoelectron spectra of water cluster anions

Water cluster anions, (H(2)O)(N)(-), are examined using mixed quantum/classical molecular dynamics based on a one-electron pseudopotential model that incorporates many-body polarization and predicts vertical electron detachment energies (VDEs) with an accuracy of ~0.1 eV. By varying the initial conditions under which the clusters are formed, we are able to identify four distinct isomer types that exhibit different size-dependent VDEs. On the basis of a strong correlation between the electron's radius of gyration and its optical absorption maximum, and extrapolating to the bulk limit (N → ∞), our analysis supports the assignment of the "isomer Ib" data series, observed in photoelectron spectra of very cold clusters, as arising from cavity-bound (H(2)O)(N)(-) cluster isomers. The "isomer I" data reported in warmer experiments are assigned to surface-bound isomers in smaller clusters, transitioning to partially embedded isomers in larger clusters. The partially embedded isomers are characterized by a partially formed solvent cavity at the cluster surface, and they are spectroscopically quite similar to internalized cavity isomers. These assignments are consistent with various experimental data, and our theoretical characterization of these isomers sheds new light on a long-standing assignment problem.     

Binding interactions and Raman spectral properties of pyridine interacting with bimetallic silver-gold clusters.

The binding interactions between pyridine and bimetallic silver-gold clusters are investigated using density functional theory (DFT). The binding energies of pyridine-bimetallic cluster complexes indicate that the bonding depends strongly on the binding site (Au or Ag atom) and bonding molecular orbitals in a given configuration. The donation of the lone-pair electrons of the nitrogen of pyridine to an appropriate unoccupied orbital of each metal cluster plays an important role. The low-lying excited states and charge-transfer states of four stable complexes of interest are calculated on the basis of a time-dependent DFT method. In nonresonance Raman scattering processes, the influence of binding interactions on the relative Raman intensity of totally symmetric pyridine vibrational modes is discussed. These calculated relative Raman intensities are compared with observed surface-enhanced Raman spectra of pyridine adsorbed on silver-gold alloy surfaces.     

Ring-slippage and multielectron redox properties of Fe/Ru/Os-bis(arene) complexes: does hapticity change really cause potential inversion?

Bis(hexamethylbenzene) complexes of the group 8 metals (Fe, Ru, Os) show surprising diversity in their electron-transfer mechanisms and associated thermodynamics for the M(II) → M(I) → M(0) redox series. In electrochemical experiments, the Fe complex exhibits normally ordered potentials separated by ～1 V, the Ru system shows nearly overlapping one-electron redox events, and Os demonstrates a one-step, two-electron transfer with a peak potential separation suggestive of highly inverted potentials. It has been conjectured that the sequential one-electron transfers observed for Fe are due to the lack of an accessible η(4):η(6) Fe(0) state, destabilizing the fully reduced species. Using an established model chemistry based on DFT, we demonstrate that the hapticity change is a consequence of the bonding throughout this transition metal triad and that apparent multielectron behavior is controlled by the vertical electron attachment component of the M(II) → M(I) redox event. Furthermore, the η(6):η(6) Fe(0) triplet state is more favorable than the hypothetical η(4):η(6) singlet state, emphasizing that the hapticity change is not sufficient for multielectron behavior. Despite both displaying two-electron redox responses, Ru and Os traverse fundamentally different mechanisms based on whether the first (Os) or second (Ru) electron transfer induces the hapticity change. While the electronic structure analysis is limited to the Fe triad here, the conceptual model that we developed provides a general understanding of the redox behavior exhibited by d(6) bis(arene) compounds.     

Lithium reduction of the bowl-shaped C60 fragment diindeno[1,2,3,4-defg;1',2',3',4'-mnop]chrysene: an interplay between experiment and calculation

Diindeno[1,2,3,4-defg;1',2',3',4'-mnop]chrysene (DIC) (one of the smallest symmetrical bowl-shaped fragments of C60) and its tetra-tert-butyl derivative are reduced with lithium metal to yield dianions and tetraanions. Due to the high degree of symmetry (C2v) of DIC and its derivative, their NMR spectra cannot be assigned using the standard two-dimensional NMR techniques. A novel carbon-edited NOESY method was used to complete the assignments of the neutral and dianion species, whereas the tetraanions are aided by DFT calculations for their assignment. Experimental charge-distribution patterns were obtained and match those of the calculations. An extension of the empirical approach for estimating the charge distribution from the 13C NMR spectra enables a direct comparison between experimentally derived charge-distribution data and the computed electron density in each of the lowest unoccupied molecular orbitals. The overall picture evolving from the orbital structure of DIC is presented and reflects the surface reactivity of C60.     

L-edge X-ray absorption spectroscopy of non-heme iron sites: experimental determination of differential orbital covalency

X-ray absorption spectroscopy has been utilized to obtain the L-edge multiplet spectra for a series of non-heme ferric and ferrous complexes. Using these data, a methodology for determining the total covalency and the differential orbital covalency (DOC), that is, differences in covalency in the different symmetry sets of the d orbitals, has been developed. The integrated L-edge intensity is proportional to the number of one-electron transition pathways to the unoccupied molecular orbitals as well as to the covalency of the iron site, which reduces the total L-edge intensity and redistributes intensity, producing shake-up satellites. Furthermore, differential orbital covalency leads to differences in intensity for the different symmetry sets of orbitals and, thus, further modifies the experimental spectra. The ligand field multiplet model commonly used to simulate L-edge spectra does not adequately reproduce the spectral features, especially the charge transfer satellites. The inclusion of charge transfer states with differences in covalency gives excellent fits to the data and experimental estimates of the different contributions of charge transfer shake-up pathways to the t(2g) and e(g) symmetry orbitals. The resulting experimentally determined DOC is compared to values calculated from density functional theory and used to understand chemical trends in high- and low-spin ferrous and ferric complexes with different covalent environments. The utility of this method toward problems in bioinorganic chemistry is discussed.     

Local ab initio methods for calculating optical band gaps in periodic systems. I. Periodic density fitted local configuration interaction singles method for polymers

We present a density fitted local configuration interaction singles (CIS) method for calculating optical band gaps in 1D-periodic systems. The method is based on the Davidson diagonalization procedure, carried out in the reciprocal space. The one-electron part of the matrix-vector products is also evaluated in the reciprocal space, where the diagonality of the Fock matrix can be exploited. The contraction of the CIS vectors with the two electron integrals is performed in the direct space in the basis of localized occupied (Wannier) and virtual (projected atomic) orbitals. The direct space approach allows to utilize the sparsity of the integrals due to the local representation and locality of the exciton. The density fitting approximation employed for the two electron integrals reduces the nominal scaling with unit cell size to O(N(4)). Test calculations on a series of prototypical systems demonstrate that the method in its present stage can be used to calculate the excitonic band gaps of polymers with up to a few dozens of atoms in the cell. The computational cost depends on the locality of the exciton, but even relatively delocalized excitons occurring in the polybiphenyl in the parallel orientation, can be routinely treated with this method.     

Parameters for excess electron transfer in DNA. Estimation using unoccupied Kohn-Sham orbitals and TD DFT

We show that the energetics and electronic couplings for excess electron transfer (EET) can be accurately estimated by using unoccupied Kohn-Sham orbitals (UKSO) calculated for neutral pi stacks. To assess the performance of different DFT functionals, we use MS-PT2 results for seven pi stacks of nucleobases as reference data. The DFT calculations are carried out by using the local spin density approximation SVWN, two generalized gradient approximation functionals BP86 and BLYP, and two hybrid functionals B3LYP and BH&HLYP. Best estimations within the UKSO approach are obtained by the B3LYP and SVWN methods. TD DFT calculations provide less accurate values of the EET parameters as compared with the UKSO data. Also, the excess charge distribution in the radical anions is well described by the LUMOs of neutral systems. In contrast, spin-unrestricted DFT calculations of radical anions considerably overestimate delocalization of the excess electron. The excellent results obtained for the ground and excited states of the radical anions (excitation energy, transition dipole moment, electronic coupling, and excess electron distribution) by using UKSO of neutral dimers suggest an efficient strategy to calculate the EET parameters for DNA pi stacks.     

Electronic structure and molecular orbital study of the first excited state of the high‐efficiency blue OLED material bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex from ab initio and TD‐B3LYP

Bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex (AlMq₂OH) is used in organic light‐emitting diodes (OLEDs) as an electron transport material and emitting layer. By means of ab initio Hartree–Fock (HF) and density functional theory (DFT) B3LYP methods, the structure of AlMq₂OH was optimized. The frontier molecular orbital characteristics and energy levels of AlMq₂OH have been analyzed systematically to study the electronic transition mechanism in AlMq₂OH. For comparison and calibration, bis(8‐quinolinolato)aluminum(III) hydroxide complex (Alq₂OH) has also been examined with these methods using the same basis sets. The lowest singlet excited state (S₁) of AlMq₂OH has been studied by the singles configuration interaction (CIS) method and time‐dependent DFT (TD‐DFT) using a hybrid functional, B3‐LYP, and the 6‐31G* basis set. The lowest singlet electronic transition (S₀ → S₁) of AlMq₂OH is π → π* electronic transitions and primarily localized on the different quinolate ligands. The emission of AlMq₂OH is due to the electron transitions from a phenoxide donor to a pyridyl acceptor from another quinolate ligand including C → C and O → N transference. Two possible electron transfer pathways are presented, one by carbon, oxygen, and nitrogen atoms and the other via metal cation Al³⁺. The comparison between the CIS‐optimized excited‐state structure with the HF ground‐state structure indicates that the geometric shift is mainly confined to the one quinolate and these changes can be easily understood in terms of the nodal patterns of the highest occupied and lowest unoccupied molecular orbitals. On the basis of the CIS‐optimized structure of the excited state, TD‐B3‐LYP calculations predict an emission wavelength of 499.78 nm. An absorption wavelength at 380.79 nm on the optimized structure of B3LYP/6‐31G* was predicted. They are comparable to AlMq2OH 485 and 390 nm observed experimentally for photoluminescence and UV‐vis absorption spectra of AlMq₂OH solid thin film on quartz, respectively. Lending theoretical corroboration to recent experimental observations and supposition, the reasons for the blue‐shift of AlMq2OH were revealed. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004