Structural and electronic properties of polyacetylene and polyyne from hybrid and Coulomb-attenuated density functionals

The bond length alternation (BLA), the highest-occupied-lowest-unoccupied (HO-LU) orbital energy gap, and the corresponding excitation energy are determined for trans-polyacetylene (PA) and polyyne (PY) using density functional theory. Results from the Coulomb-attenuated CAM-B3LYP functional are compared with those from the conventional BHHLYP and B3LYP hybrid functionals. BLA values and HO-LU gaps are determined using both finite oligomer and infinite chain calculations, subject to periodic boundary conditions. TDDFT excitation energies are determined for the oligomers. The oligomer excitation energies and HO-LU gaps are then used, in conjunction with the infinite chain HO-LU gap, to estimate the infinite chain excitation energy. Overall, BHHLYP and CAM-B3LYP give BLA values and excitation energies that are larger and more accurate than those obtained using B3LYP. The results highlight the degree to which excitation energies can be approximated using the HO-LU gaps-at the infinite limit, this approximation works well for B3LYP, but not for the other functionals, where the HO-LU gap is significantly larger. The study provides further evidence for the high-quality theoretical predictions that can be obtained from the CAM-B3LYP functional.     

Density functional and spin-orbit ab initio study of CF3Br: molecular properties and electronic curve crossing

Quantum chemical calculations of CF(3)Br and the CF(3) radical are performed using density functional theory (DFT) and time-dependent DFT (TDDFT). Molecular structures, vibrational frequencies, dipole moment, bond dissociation energy, and vertical excitation energies of CF(3)Br are calculated and compared with available experimental results. The performance of six hybrid and five hybrid meta functionals in DFT and TDDFT calculations are evaluated. The ωB97X, B3PW91, and M05-2X functionals give very good results for molecular structures, vibrational frequencies, and vertical excitation energies, respectively. The ωB97X functional calculates well the dipole moment of CF(3)Br. B3LYP, one of the most widely used functionals, does not perform well for calculations of the C-Br bond length, bond dissociation energy, and vertical excitation energies. Potential energy curves of the low-lying excited states of CF(3)Br are obtained using the multiconfigurational spin-orbit ab initio method. The crossing point between 2A(1) and 3E states is located near the C-Br bond length of 2.45 ?. Comparison with CH(3)Br shows that fluorination does not alter the location of the crossing point. The relation between the calculated potential energy curves and recent experimental result is briefly discussed.     

Ferric complexes of 3-hydroxy-4-pyridinones characterized by density functional theory and Raman and UV-vis spectroscopies

Deferiprone and other 3-hydroxy-4-pyridinones are used in metal chelation therapy of iron overload. To investigate the structure and stability of these compounds in the natural aqueous environment, ferric complexes of deferiprone and amino acid maltol conjugates were synthesized and studied by computational and optical spectroscopic methods. The complexation caused characteristic intensity changes, a 300× overall enhancement of the Raman spectrum, and minor changes in UV-vis absorption. The spectra were interpreted on the basis of density functional theory (DFT) calculations. The CAM-B3LYP and ωB97XD functionals with CPCM solvent model were found to be the most suitable for simulations of the UV-vis spectra, whereas B3LYP, B3LYPD, B3PW91, M05-2X, M06, LC-BLYP, ωB97XD, and CAM-B3LYP functionals were all useful for simulation of the Raman scattering. Characteristic Raman band frequencies for 3-hydroxy-4-pyridinones were assigned to molecular vibrations. The computed conformer energies consistently suggest the presence of another isomer of the deferiprone-ferric complex in solution, in addition to that found previously by X-ray crystallography. However, the UV-vis and Raman spectra of the two species are similar and could not be resolved. In comparison to UV-vis, the Raman spectra and their combination with calculations appear more promising for future studies of iron sequestrating drugs and artificial metalloproteins as they are more sensitive to structural details.     

Theoretical Study on the Optical Properties for 2,7- and 3,6-Linked Carbazole Trimers by Time-dependent Density Functional Theory

Electronic properties, such as HOMO and LUMO energies, band gaps, ionization potential (IP) and electron affinity (EA) of 2,7- and 3,6-linked carbazole trimers, two conjugated oligomers with different linkages of carbazole, were studied by the density functional theory with Becke-Lee-Young-Parr composite exchange correlation functional (B3LYP). The absorption spectra of these compounds were also investigated by time-dependent density functional theory (TD-DFT) with 6-31G* basis set. The calculated results indicated that the HOMO and LUMO of the 2,7- and 3,6-linked carbazole trimers are both slightly destabilized on going from methyl substitution to sec-butyl substitution. Both IP and EA exhibit their good hole-transporting but poor electron- accepting ability. The presence of alkyl groups on the nitrogen atoms does not affect the intra-chain electronic delocalization along the molecular frame. Thus no significant effect on the band gap and absorption spectra of compounds has been found.     

Ethynyl Substitution Effect on the Electronic Excitation Spectra of Aniline

We made an extended study on the structure and properties of the low-lying electronic states of ethynyl substituted aniline and their cations. We performed these calculations using density functional theory method(B3LYP and CAM-B3LYP DFT) and the complete active space self-consistent field(CASSCF) approach in connection with the aug-cc-pVZ Dunning's basis sets and concerted ANO-L-VDZP basis sets. Our results included their equilibrium geometries, the vertical excitation spectra and the vertical and adiabatic ionization energies. The effect of ethynyl substitution on the electronic structure and the spectroscopy of aniline was probed.     

TD-CI simulation of the electronic optical response of molecules in intense fields II: comparison of DFT functionals and EOM-CCSD

Time-dependent configuration interaction (TD-CI) simulations can be used to simulate molecules in intense laser fields. TD-CI calculations use the excitation energies and transition dipoles calculated in the absence of a field. The EOM-CCSD method provides a good estimate of the field-free excited states but is rather expensive. Linear-response time-dependent density functional theory (TD-DFT) is an inexpensive alternative for computing the field-free excitation energies and transition dipoles needed for TD-CI simulations. Linear-response TD-DFT calculations were carried out with standard functionals (B3LYP, BH&HLYP, HSE2PBE (HSE03), BLYP, PBE, PW91, and TPSS) and long-range corrected functionals (LC-ωPBE, ωB97XD, CAM-B3LYP, LC-BLYP, LC-PBE, LC-PW91, and LC-TPSS). These calculations used the 6-31G(d,p) basis set augmented with three sets of diffuse sp functions on each heavy atom. Butadiene was employed as a test case, and 500 excited states were calculated with each functional. Standard functionals yield average excitation energies that are significantly lower than the EOM-CC, while long-range corrected functionals tend to produce average excitation energies slightly higher. Long-range corrected functionals also yield transition dipoles that are somewhat larger than EOM-CC on average. The TD-CI simulations were carried out with a three-cycle Gaussian pulse (ω = 0.06 au, 760 nm) with intensities up to 1.26 × 10(14) W cm(-2) directed along the vector connecting the end carbons. The nonlinear response as indicated by the residual populations of the excited states after the pulse is far too large with standard functionals, primarily because the excitation energies are too low. The LC-ωPBE, LC-PBE, LC-PW91, and LC-TPSS long-range corrected functionals produce responses comparable to EOM-CC.     

On the performance of range-separated hybrid in computations of dynamic quadratic polarizability of solution-phase organic molecules: a comparison to MP2(Full) calculation

We have assessed six range-separated hybrids (CAM-B3LYP, LC-PBE, ωB97, ωB97X, ωB97X-D and M11) in evaluating the dynamic quadratic polarizability along the direction of the molecular dipole moment, β _μ (?2ν; ν; ν; 0), by comparing with MP2(Full)-calculated data. There are 104 solution-phase organic molecules included in the set. Both geometry optimizations and β _μ (?2ν; ν; ν; 0) computations have been carried out by the same functional. Calculated results showed that simultaneous modeling geometrical and electronic effects improved significantly the nonlinear optical (NLO) property of interest. More generally, our numerical assessment revealed that the CAM-B3LYP functional could be the reasonable choice for NLO property computations.     

Theoretical investigation of optical and electronic property modulations of pi-conjugated polymers based on the electron-rich 3,6-dimethoxy-fluorene unit

Poly(fluorene)-type materials are widely used in polymer-based emitting devices. One of the drawbacks of light-emitting diodes based on polyfluorene derivatives is the injection of holes from the anode due to the high ionization potential (IP) of most derivatives. Substitution by electron-donating alkoxy substituents or by adding charge carriers on the conjugated polymer's backbone produces a remarkable influence on its electrical and optical properties. In this contribution, we apply quantum-chemical techniques to investigate a family of pi-conjugated polymers with substituted dimethoxy groups at the 3,6 positions of the fluorene ring, namely, poly(2,7-(3,6-dimethoxy-fluorene)(PDMOF), poly(2,7-(3,6-dimethoxy-fluorene)-co-alt-fluorene (PDMOFF), and poly(2,7-(3,6-dimeth-oxy-fluorene)-co-alt-2,5-thiophene (PDMOFT). The electronic properties of the neutral molecules, HOMO-LUMO gaps (Delta(H)(-)(L)), in addition to the positive and negative ions, are studied using the B3LYP functional. The lowest excitation energies (E(g)) and the maximal absorption wavelength lambda(abs) of PDMOF, PDMOFF, and PDMOFT are studied by employing time-dependent density functional theory (TD-DFT) and the ZINDO semiempirical method. The IP, EA, and E(g) values of each polymer were obtained by extrapolating those of the oligomers to the inverse chain length equal to zero ((1)/(n)() = 0). The influence of the presence of methoxy groups on the fluorene moiety on the ionization potential is especially emphasized. The outcomes show that the HOMO energies of these systems under study increase by about 0.4 eV and the IP values decrease by about 0.3 eV compared to those of the corresponding polyfluorene. Both effects result in a reduction of the energy barrier for the injection of holes in related polymeric light-emitting devices and should contribute to the enhancement of their performances. Because of the cooperation with thiophene in PDMOFT, which results in a good planar conformation, both the hole-creating and electron-accepting abilities are improved.     

Computational study of radical cations of saturated compounds with sigma-type and pi-type N-N bonds.

Geometrical and electronic properties have been calculated and are compared with experimental data for three saturated diaza compounds and their radical cations and dications. The molecular geometries in the different oxidation states are consistently reproduced very well using the B3PW91 and B3LYP three-parameter density functional methods, with a modest 6-31G* basis set. The performance of the pure density functionals BLYP and BPW91 is less satisfactory. The Hartree-Fock method yields excellent results in some cases but poor results in others. Ionization potentials and electron-nuclear hyperfine interactions are reproduced moderately well with B3LYP and B3PW91. Electronic excitation energies calculated with time-dependent density functional theory agree very well with experiment in most cases. For 2,7-diazatetracyclo[6.2.2.2(3,6).0(2,7)]tetradecane 2 and its radical cation and dication, the reorganization parameters for self-electron exchange were calculated and compared with experimental and earlier computed data. The calculations allow a good estimate of the different contributions to the energy barrier, i.e., the internal and solvent reorganization energies and the work term in the case of 2+/2++.     

Predictive power of long-range corrected functionals on the spectroscopic properties of tetrapyrrole derivatives for photodynamic therapy

Porphyrin and chlorin based compounds possess promising properties to be utilized as photosensitizers in photodynamic therapy (PDT). However, the photosensitizers available on the market today are not ideal for use in PDT, which has emphasized the need for new photosensitizers with improved photodynamic properties to be developed. Computational drug-design can be utilized in the search for improved pharmaceutical compounds, provided that the methods used are able to reproduce experimental data. In the present study we investigated, by the use of time-dependent density functional theory (TD-DFT), the performance of the long-range corrected functionals ωB97, ωB97X and ωB97XD on their ability to predict low-lying singlet excitations (>600 nm) of a set of well-known photosensitizing compounds. It was found that ωB97X reproduced the experimental red-most absorption band most satisfactorily. The use of either B3LYP, ωB97XD or M06 in geometry optimizations has a minor effect on the spectra in most cases. Calculated energy differences between the optimized singlet ground states and optimized first excited triplet states show consistent and overall higher triplet state energies for B3LYP, M06, and PBE0 compared with ωB97, ωB97X, and ωB97XD. The calculated triplet state energies are, however, sufficient to generate singlet oxygen in most cases.     

DFT study of conjugational electronic structures of aminoalkyl end-capped oligothiophenes up to octamers

The molecular geometries and electronic properties of a series of bis(aminoalkyl) end-capped oligothiophenes (BRnTs) were investigated by means of the density functional theory (DFT). The calculations were performed on dimers up to octamers in the neutral and ionic species using the B3LYP/6-31G(d,p) level of theory. The results obtained show that the conjugated systems in the p- and n-doped oligomers had more aromaticity, with expanded and planar chains. The calculated energy gap values between the frontier molecular orbitals for the end-capped oligomers were larger than those for the unsubstituted oligomers, in which with increase in the oligomer chain length, the conduction band gap decreased. The calculated first excitation energies of BRnTs at the TD-B3LYP/6-31G(d,p) level indicated that both doped oligomers (p- and n-type) had lower excitation energies than the neutral states, and that they displayed red shifts in their absorption spectra. Moreover, the results obtained for the natural bond orbital (NBO) analysis showed that closing the end-side oligothiophene chains with the aminoalkyl groups eased the hole or electron transfer, owning to better charge delocalization through the backbone structures of BRnTs.     

DFT investigation of endohedral boron oxide nanocapsules: Encapsulation of He, Ne, Ar, H, N, and Cl atoms

The electronic structure and stabilization energy of spherical and pyramidal shapes of boron oxide nanocapsules (X@B₂₀O₃₀, X = He, Ne, Ar, H, N, Cl) were investigated by long-range and dispersion corrected density functional theory (DFT + Disp) including CAM-B3LYP, B3LYP-D3, ωB97X-D and B2PLYP-D methods. Based on these calculations, the formation of nanocapsules is an exothermic process (except for Cl@B₂₀O₃₀). The spherical boron oxide nanocapsules are mainly stabilized by dispersion, while the stability of pyramidal complexes is mainly due to monomer relaxation energy. The theoretical results obtained in this work show that the boron oxide capsule is a good potential candidate for gas storage.     

Assessment of a Coulomb-attenuated exchange-correlation energy functional

The recently proposed CAM-B3LYP exchange-correlation energy functional, based on a partitioning of the r operator in the exchange interaction into long- and short-range components, is assessed for the determination of molecular thermochemistry, structures, and second order response properties. Rydberg and charge transfer excitation energies and static electronic polarisabilities are notably improved over the standard B3LYP functional; classical reaction barriers also improve. Ionisation potentials, bond lengths, NMR shielding constants and indirect spin-spin coupling constants are comparable with the two functionals. CAM-B3LYP atomisation energies and diatomic harmonic vibrational wavenumbers are less accurate than those of B3LYP. Future research directions are outlined.     

First hyperpolarizability of polymethineimine with long-range corrected functionals

Using the long-range corrected (LC) density functional theory (DFT) scheme introduced by Iikura et al. [J. Chem. Phys. 115, 3540 (2001)] and the Coulomb-attenuating model (CAM-B3LYP) of Yanai et al. [Chem. Phys. Lett. 393, 51 (2004)], we have calculated the longitudinal dipole moments and static electronic first hyperpolarizabilities of increasingly long polymehtineimine oligomers. For comparison purposes Hartree-Fock (HF), Moller-Plesset perturbation theory (MP2), and conventional pure and hybrid functionals have been considered as well. HF, generalized gradient approximation (GGA), and conventional hybrids provide too large dipole moments for long oligomers, while LC-DFT allows to reduce the discrepancy with respect to MP2 by a factor of 3. For the first hyperpolarizability, the incorrect evolution with the chain length predicted by HF is strongly worsened by BLYP, Perdew-Burke-Ernzerhof (PBE), and also by B3LYP and PBE0. On the reverse, LC-BLYP and LC-PBE hyperpolarizabilities are correctly predicted to be positive (but for the two smallest chains). Indeed, for medium and long oligomers LC hyperpolarizabilities are slightly smaller than MP2 hyperpolarizabilities, as it should be. CAM-B3LYP also strongly improves the B3LYP results, though a bit less impressively for small chain lengths. The present study demonstrates the efficiency of long-range DFT, even in very pathological cases.     

Computational study of static first hyperpolarizability of donor-acceptor substituted (E)-benzaldehyde phenylhydrazone

Various hybrid functionals (B3LYP, B97-2, PBE0, BMK, BH&HLYP, CAM-B3LYP, and LC-ωPBE) implemented in density functional theory were applied to give estimate of static first hyperpolarizabilty (β(0)) of (E)-benzaldehyde phenylhydrazone designated as (E)-BPH. Against those of MP2 computations as a function of the underlying density functional, good agreement was obtained with the BH&HLYP and CAM-B3LYP functionals. The LC-ωPBE functional and the B3LYP, PBE0, B97-2, and BMK functionals underestimated and overestimated β(0), respectively. The basis set effect on the calculated β(0) was also investigated. It turned out that the 6-311+G(2d,p) basis set provided excellent converged value of β(0). On the basis of the calculated results, we investigated the substituent effect on β(0) of donor-acceptor (D-A) substituted (E)-BPH systematically by using the BH&HLYP and CAM-B3LYP computations with the 6-311+G(2d,p) basis set. We proposed a Zwitterion structure to explain the calculated trend in the substituent effect and the enhanced hyperpolarizability of type II compounds (A-(E)-BPH-D) than type I compounds (D-(E)-BPH-A). Natural bonding orbital analysis carried out at BH&HLYP/6-311+G(2d,p)//B3LYP/6-31G(2df,p) level of theory substantiated the claim.     

On the performance of quantum chemical methods to predict solvatochromic effects: the case of acrolein in aqueous solution

The performance of the Hartree-Fock method and the three density functionals B3LYP, PBE0, and CAM-B3LYP is compared to results based on the coupled cluster singles and doubles model in predictions of the solvatochromic effects on the vertical n-->pi* and pi-->pi* electronic excitation energies of acrolein. All electronic structure methods employed the same solvent model, which is based on the combined quantum mechanics/molecular mechanics approach together with a dynamical averaging scheme. In addition to the predicted solvatochromic effects, we have also performed spectroscopic UV measurements of acrolein in vapor phase and aqueous solution. The gas-to-aqueous solution shift of the n-->pi* excitation energy is well reproduced by using all density functional methods considered. However, the B3LYP and PBE0 functionals completely fail to describe the pi-->pi* electronic transition in solution, whereas the recent CAM-B3LYP functional performs well also in this case. The pi-->pi* excitation energy of acrolein in water solution is found to be very dependent on intermolecular induction and nonelectrostatic interactions. The computed excitation energies of acrolein in vacuum and solution compare well to experimental data.     

Theoretical studies on the electronic and optical properties of two new alternating fluorene/carbazole copolymers

Poly(fluorene)-type materials are widely used in polymer-based emitting devices. During operation there appears, however, an additional emission peak at around 2.3 eV, leading to both a color instability and reduced efficiency. The incorporation of the carbazole units has been proven to efficiently suppress the keto defect emission. In this contribution, we apply quantum-chemical techniques to investigate two series of alternating fluorene/carbazole oligomers and copolymers poly[2,7-(N-(2-methyl)-carbazole)-co-alt-2,7-m(9,9-dimethylfluorene)], namely, PFmCz (m = 1,2) and gain a detailed understanding of the influence of carbazole units on the electronic and optical properties of fluorene derivatives. The electronic properties of the neutral molecules, HOMO-LUMO gaps (Delta(H-L)), in addition to the positive and negative ions, are studied using B3LYP functional. The lowest excitation energies (E(g)s) and the maximal absorption wavelength lambda(abs) of PFmCz (m = 1,2) are studied, employing the time-dependent density functional theory (TD-DFT). The properties of the two copolymers, such as Delta(H-L), E(g), IPs, and EAs were obtained by extrapolating those of the oligomers to the inverse chain length equal to zero (1/n = 0). The outcomes showed that the carbazole unit is a good electron-donating moiety for electronic materials, and the incorporation of carbazole into the polyfluorene (PF) backbone resulted in a broadened energy gap and a blue shift of both the absorption and photoluminescence emission peaks. Most importantly, the HOMO energies of PF1Cz and PF2Cz are both a higher average (0.4 eV) than polyfluorene (PF), which directly results in the decreasing of IPs of about 0.2 eV more than PF, indicating that the carbazole units have significantly improved the hole injection properties of the copolymers. In addition, the energy gap tends to broaden and the absorption and emission peaks are gradually blue-shifted to shorter wavelengths with an increase in the carbazole content in the copolymers. This is due to the interruption of the longer conjugation length of the backbone in the (F1Cz)(n) series.     

Low-lying excited states of C120 and C151: a multireference perturbation theory study

Excited states of two 7-aminocoumarin derivatives, coumarin 120 (7-amino-4-methylcoumarin) and coumarin 151 (7-amino-4-trifluoromethylcoumarin), were investigated using generalized multiconfigurational quasidegenerate perturbation theory (GMC-QDPT), multiconfigurational quasidegenerate perturbation theory (MC-QDPT) and time-dependent density functional theory (TDDFT) with the B3LYP and CAM-B3LYP functionals. The absorption and fluorescence spectra of C120 and C151 were calculated. We elucidated the characters of the low-lying states of C120 and C151. The absorption spectra calculated with GMC-QDPT and TDDFT B3LYP agreed well with the experimental data, while for the fluorescence spectra, the TDDFT calculations overestimated the fluorescence spectra compared to GMC-QDPT calculations. Utilizing active spaces with large numbers of electrons and orbitals for reference functions, GMC-QDPT showed a better performance than MC-QDPT with a complete active space self-consistent field (CASSCF) reference of active space with smaller number of electrons and orbitals. In our gas phase calculation, we found that the optimized structures for the first excited states have a planar amino group with a CN single bond, while the amino group is pyramidal in the ground state.     

Assessment of the efficiency of long-range corrected functionals for some properties of large compounds

Using the long-range correction (LC) density functional theory (DFT) scheme introduced by Iikura et al. [J. Chem. Phys. 115, 3540 (2001)] and the Coulomb-attenuating model (CAM-B3LYP) of Yanai et al. [Chem. Phys. Lett. 393, 51 (2004)], we have calculated a series of properties that are known to be poorly reproduced by standard functionals: Bond length alternation of pi-conjugated polymers, polarizabilities of delocalized chains, and electronic spectra of extended dyes. For each of these properties, we present cases in which traditional hybrid functionals do provide accurate results and cases in which they fail to reproduce the correct trends. The quality of the results is assessed with regard to experimental values and/or data arising from electron-correlated wave function approaches. It turns out that (i) both LC-DFT and CAM-B3LYP provide an accurate bond length alternation for polyacetylene and polymethineimine, although for the latter they decrease slightly too rapidly with chain length. (ii) The LC generalized gradient approximation and MP2 polarizabilities of long polyphosphazene and polymethineimine oligomers agree almost perfectly. In the same way, CAM-B3LYP corrects the major part of the B3LYP faults. (iii) LC and CAM techniques do not help in correcting the nonrealistic evolution with chain length of the absorption wavelengths of cyanine derivatives. In addition, though both schemes significantly overestimate the ground to excited state transition energy of substituted anthraquinone dyes, they provide a more consistent picture once a statistical treatment is performed than do traditional hybrid functionals.