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.     

para‐Phenylene‐Bridged Spirobi(triarylamine) Dimer with Four Perpendicularly Linked Redox‐Active π Systems

para‐Phenylene‐bridged spirobi(triarylamine) dimer 2 , in which π conjugation through four redox‐active triarylamine subunits is partially segregated by the unique perpendicular conformation, was prepared and characterized by structural, electrochemical, and spectroscopic methods. Quantum chemical calculations (DFT and CASSCF) predicted that the frontier molecular orbitals of 2 are virtually fourfold degenerate, so that the oxidized states of 2 can give intriguing electronic and magnetic properties. In fact, the continuous‐wave ESR spectroscopy of radical cation 2 ^.+ showed that the unpaired electron was trapped in the inner two redox‐active dianisylamine subunits, and moreover was fully delocalized over them. Magnetic susceptibility measurements and pulsed ESR spectroscopy of the isolated salts of 2 , which can be prepared by treatment with SbCl₅, revealed that the generated tetracation 2 ⁴⁺ decomposed mainly into a mixture of 1) a decomposed tetra(radical cation) consisting of a tri(radical cation) moiety and a trianisylamine radical cation moiety (≈75 %) and 2) a diamagnetic quinoid dication in a tetraanisyl‐p‐phenylendiamine moiety and two trianisylamine radical cation moieties (≈25 %). Furthermore, the spin‐quartet state of the tri(radical cation) moiety in the decomposed tetra(radical cation) was found to be in the ground state lying 30 cal mol^?1 below the competing spin‐doublet state.     

On orthogonality constrained multiple core‐hole states and optimized effective potential method

An attempt to construct a multiple core‐hole state within the optimized effective potential (OEP) methodology is presented. In contrast to the conventional Δ‐self‐consistent field method for hole states, the effects of removing an electron is achieved using some orthogonality constraints imposed on the orbitals so that a Slater determinant describing a hole state is constrained to be orthogonal to that of a neutral system. It is shown that single, double, and multiple core‐hole states can be treated within a unified framework and can be easily implemented for atoms and molecules. For this purpose, a constrained OEP method proposed earlier for excited states (Glushkov and Levy, J. Chem. Phys. 2007, 126, 174106) is further developed to calculate single and double core ionization energies using a local effective potential expressed as a direct mapping of the external potential. The corresponding equations, determining core‐hole orbitals from a one‐particle Schrödinger equation with a local potential as well as correlation corrections derived from the second‐order many‐body perturbation theory are given. One of the advantages of the present direct mapping formulation is that the effective potential, which plays the role of the Kohn–Sham potential, has the symmetry of the external potential. Single and double core ionization potentials computed with the presented scheme were found to be in agreement with data available from experiment and other calculations. We also discuss core‐hole state local potentials for the systems studied. © 2012 Wiley Periodicals, Inc.     

B3LYP,BLYP and PBE DFT band structures of the nucleotide base stacks

DFT crystal orbital (band structure) calculations have been performed for the nucleotide base stacks of cytosine, thymine, adenine, and guanine arranged in DNA B geometry. The band structures obtained with PBE, BLYP, and B3LYP functionals are presented and compared to other related experimental and theoretical results. The influence of the quality of the basis set on the fundamental gap values was also investigated using Clementi's double ζ, 6‐31G and 6‐31G* basis sets. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

A semiempirical long-range corrected exchange correlation functional including a short-range Gaussian attenuation (LCgau-B97)

We applied an improved long‐range correction scheme including a short‐range Gaussian attenuation (LCgau) to the Becke97 (B97) exchange correlation functional. In the optimization of LCgau‐B97 functional, the linear parameters are determined by least squares fitting. Optimizing μ parameter (0.2) that controls long‐range portion of Hartree‐Fock (HF) exchange to excitation energies of large molecules (Chai and Head‐Gordon, J Chem Phys 2008, 128, 084106) and additional short‐range Gaussian parameters (a = 0.15 and k = 0.9) that controls HF exchange inclusion ranging from short‐range to mid‐range (0.5–3 Å) to ground state properties achieved high performances of LCgau‐B97 simultaneously on both ground state and excited state properties, which is better than other tested semiempirical density functional theory (DFT) functionals, such as ωB97, ωB97X, BMK, and M0x‐family. We also found that while a small μ value (～0.2) in LC‐DFT is appropriate to the local excitation and intramolecular charge‐transfer excitation energies, a larger μ value (0.42) is desirable in the Rydberg excitation‐energy calculations. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Theoretical study of CH…O hydrogen bond in proton transfer reaction of glycine

Density functional theory (DFT) calculations are used to study the strength of the CH…O H‐bond in the proton transfer reaction of glycine. Comparison has been made between four proton transfer reactions (ZW1, ZW2, ZW3, SCRFZW) in glycine. The structural parameters of the zwitterionic, transition, and neutral states of glycine are strongly perturbed when the proton transfer takes place. It has been found that the interaction of water molecule at the side chain of glycine is high in the transition state, whereas it is low in the zwitterionic and neutral states. This strongest multiple hydrogen bond interaction in the transition state reduces the barrier for the proton transfer reaction. The natural bond orbital analysis have also been carried out for the ZW2 reaction path, it has been concluded that the amount of charge transfer between the neighboring atoms will decide the strength of H‐bond. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Comparing ab initio density-functional and wave function theories: the impact of correlation on the electronic density and the role of the correlation potential

The framework of ab initio density-functional theory (DFT) has been introduced as a way to provide a seamless connection between the Kohn-Sham (KS) formulation of DFT and wave-function based ab initio approaches [R. J. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005)]. Recently, an analysis of the impact of dynamical correlation effects on the density of the neon atom was presented [K. Jankowski, K. Nowakowski, I. Grabowski, and J. Wasilewski, J. Chem. Phys. 130, 164102 (2009)], contrasting the behaviour for a variety of standard density functionals with that of ab initio approaches based on second-order M?ller-Plesset (MP2) and coupled cluster theories at the singles-doubles (CCSD) and singles-doubles perturbative triples [CCSD(T)] levels. In the present work, we consider ab initio density functionals based on second-order many-body perturbation theory and coupled cluster perturbation theory in a similar manner, for a range of small atomic and molecular systems. For comparison, we also consider results obtained from MP2, CCSD, and CCSD(T) calculations. In addition to this density based analysis, we determine the KS correlation potentials corresponding to these densities and compare them with those obtained for a range of ab initio density functionals via the optimized effective potential method. The correlation energies, densities, and potentials calculated using ab initio DFT display a similar systematic behaviour to those derived from electronic densities calculated using ab initio wave function theories. In contrast, typical explicit density functionals for the correlation energy, such as VWN5 and LYP, do not show behaviour consistent with this picture of dynamical correlation, although they may provide some degree of correction for already erroneous explicitly density-dependent exchange-only functionals. The results presented here using orbital dependent ab initio density functionals show that they provide a treatment of exchange and correlation contributions within the KS framework that is more consistent with traditional ab initio wave function based methods.     

Local orbitals for the study of the π→π* excitation in polyenes

Localized orbitals have recently been employed in large ab initio calculations, but their use has generally been restricted to ground‐state problems. In this work, we analyze the molecular orbitals of the excited states, optimized with a recently proposed local procedure. This method produces local orbitals of the CAS–SCF type, which permits its application to the study of excited states. In particular, we focus on the π→π* triplet excited state in polyenes, calculated using a 2/2 CAS space which includes two electrons in one π and one π* orbitals. In small polyenes, these two singly occupied active orbitals are delocalized all along the molecule. The extent of the delocalization is analyzed by studying polyenes of increasing size. Different polyenes have been studied, going from C₁₄H₁₆ to the C₇₀H₇₂ polyene. The relation of the π→π* excitation with the cation and anion systems is also discussed. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Electronic spectroscopy and computational studies of glutathionylco(III)balamin

We have studied glutathionylcobalamin (GS-Cbl) by optical spectroscopy and with density functional theory (DFT) and time-dependent DFT (TD-DFT) electronic structure methods of truncated geometric models. We examined the geometric structure of the models by comparison of DFT calculations with recent high-resolution experimental X-ray structure data ( Hannibal, L. et al. Inorg. Chem. 2010, 49, 9921) for GS-Cbl, and we examined the TD-DFT excitation simulations by comparison of the models with measured optical spectra. The calculations employed the B3LYP hybrid functional and the nonhybrid BP86 functional in both vacuum and water (conductor polarized continuum model (cpcm)) with the 6-311G(d,p) basis set. The optimized geometric structure calculations for six truncated models were made by varying the chemical structure, solvent model, and the two DFT functionals. All showed similar geometry. Charge decomposition analysis (CDA) and extended charge decomposition analysis (ECDA), especially with BP86 shows the similar charge transfer nature of the Co-S bond in GS-Cbl and the Co-C bond in CH(3)Cbl. Mayer and Wiberg bond orders illustrate the similar covalent nature of the two bonds. Finally, absolute optical spectral simulation calculations were compared with the experimental UV-visible extinction spectrum and the electronic circular dichroism (ECD) differential extinction spectrum. The BP86 method shows more spectral features, and the best fit was found for a GS-Cbl model with 5,6-dimethylbenzimidazole at the BP86/6-311G(d,p) level with a water cpcm solvent model. The excited state transitions were investigated with Martin's natural transition orbitals (NTOs). The BP86 calculations also showed π bonding interactions between Co and the axial S of the GS- ligand in the molecular orbitals (MOs) and NTOs.     

Parallel implementation of γ‐point pseudopotential plane‐wave DFT with exact exchange

Semi‐local functionals commonly used in density functional theory (DFT) studies of solids usually fail to reproduce localized states such as trapped holes, polarons, excitons, and solitons. This failure is ascribed to self‐interaction which creates a Coulomb barrier to localization. Pragmatic approaches in which the exchange correlation functionals are augmented with small amount of exact exchange (hybrid‐DFT, e.g., B3LYP and PBE0) have shown to promise in rectifying this type of failure, as well as producing more accurate band gaps and reaction barriers. The evaluation of exact exchange is challenging for large, solid state systems with periodic boundary conditions, especially when plane‐wave basis sets are used. We have developed parallel algorithms for implementing exact exchange into pseudopotential plane‐wave DFT program and we have implemented them in the NWChem program package. The technique developed can readily be employed in Γ‐point plane‐wave DFT programs. Furthermore, atomic forces and stresses are straightforward to implement, making it applicable to both confined and extended systems, as well as to Car‐Parrinello ab initio molecular dynamic simulations. This method has been applied to several systems for which conventional DFT methods do not work well, including calculations for band gaps in oxides and the electronic structure of a charge trapped state in the Fe(II) containing mica, annite. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Hole localization in [AlO4]0 defects in silica materials

First-principles calculations based on cluster models have been performed to investigate the ground state and the optically excited states of the [AlO(4)](0) hole in alpha-quartz and in the siliceous zeolite ZSM-5. The structure and spectroscopic properties of this defect have been studied using the recently developed Becke88-Becke95 one-parameter model for kinetics (BB1K) functional of Zhao et al., [J. Phys. Chem. A 108, 2715 (2004)]. Our results show that the BB1K method is significantly more reliable and more accurate than the standard density-functional theory (DFT) functionals at reproducing the localized spin density on one oxygen atom and the hyperfine coupling constants associated with the hole. Furthermore, we find that the BB1K results are in close agreement with experiments, and with the self-interaction-free unrestricted Hartree-Fock (UHF) and unrestricted second-order M?ller-Plesset perturbation theory (UMP2) calculations. For the first time, we present results of the ground-state paramagnetic properties of the Al defect in ZSM-5. Similar to the theoretical work for defective alpha-quartz, we find that the BB1K, UHF, UHFLee-Yang-Parr, and UMP2 calculations show a localized hole on one oxygen neighboring the Al, while even the best to date thermochemically derived hybrid generalized gradient approximation density-functional, B97-2, predicts a different model where the hole is distributed over two oxygen. We have further considered the optical transitions of the [AlO(4)](0) center in alpha-quartz and ZSM-5. In both systems, our BB1K time-dependent density-functional theory (TDDFT) and configuration interaction singles (CIS) calculations predict that the most likely transition involves electron transfer from the hole-bearing oxygen to other neighboring oxygen ions. This reinforces the experimental conclusions obtained for defective alpha-quartz. Notably, the two lowest, most dominant excitation energies calculated by BB1K-TDDFT (1.99 and 3.03 eV) show excellent agreement with experiment (1.96 and 2.85 eV [B. K. Meyer, J.M. Spaeth, and J.A. Weil, J. Phys. C: Solid State Phys. 17, L31 (1987)]) clearly outperforming the CIS method and other DFT calculations available in the literature.     

Accelerating the convergence of the total energy evaluation in density functional theory calculations

A special feature of the Strutinsky shell correction method (SCM) [D. Ullmo et al., Phys. Rev. B 63, 125339 (2001)] and the recently proposed orbital-corrected orbital-free density functional theory (OO-DFT) [B. Zhou and Y. A. Wang, J. Chem. Phys. 124, 081107 (2006)] is that the second-order corrections are incorporated in the total energy evaluation. In the SCM, the series expansion of the total electronic energy is essentially the Harris functional with its second-order correction. Unfortunately, a serious technical problem for the SCM is the lack of the exact Kohn-Sham (KS) density rho KS(r) required for the evaluation of the second-order correction. To overcome this obstacle, we design a scheme that utilizes the optimal density from a high-quality density mixing scheme to approximate rho KS(r). Recently, we proposed two total energy density functionals, i.e., the Zhou-Wang-lambda (ZW lambda) and the Wang-Zhou-alpha (WZ alpha) functionals, for use in the OO-DFT method. If the two interpolation parameters, lambda and alpha, are chosen to allow the second-order errors of the ZW lambda and the WZ alpha functionals to vanish, these two functionals reduce to the Hohenberg-Kohn-Sham functional with its second-order correction. Again, the optimal density from a high-quality density mixing scheme is used to approximate rho KS(r) in the evaluation of lambda and alpha. This approach is tested in iterative KS-DFT calculations on systems with different chemical environments and can also be generalized for use in other iterative first-principles quantum chemistry methods.     

Hole transfer energetics in structurally distorted DNA: the nucleosome core particle

The dynamics of long-range hole transport (HT) through DNA are critically dependent on the relative energies of guanine radical cation states. Electrostatic contacts with protein fragments and changes in the secondary structure of the DNA helix are expected to directly influence the stability of a guanine radical cation. This expectation is especially relevant when considering DNA HT in the eukaryotic nucleus, where DNA is condensed into nucleosome core particles (NCPs), the fundamental building blocks of chromatin. Using quantum-chemical calculations, we consider how the electrostatic interactions between the DNA nucleobases and the surrounding protein and water atoms and the structural changes in DNA arising from compaction into a NCP affect the energetics of hole transfer between guanine sites. We find that structural distortions of DNA can have dramatic consequences for the stability of a guanine radical cation, and therefore, these effects must be taken into account during the modeling of in vivo DNA HT and in the interpretation of experimental findings. When the electrostatic potential arising from the water and basic histone proteins is included we find that DNA-histone contacts, particularly between arginine residues and the DNA minor groove, destabilize the hole state on specific guanine residues. Therefore, contacts between the DNA nucleobases and basic amino acids have the potential to perturb the sites of preferred hole stability in DNA.     

Implementation of empirical dispersion corrections to density functional theory for periodic systems

A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT‐D3) is implemented in the plane‐wave program package VASP. The DFT‐D3 implementation is compared with an implementation of the earlier DFT‐D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT‐D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three‐center terms as implemented in the original molecular DFT‐D3 version is investigated. The empirical parameters are taken from the original DFT‐D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid‐state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT‐D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106). © 2012 Wiley Periodicals, Inc.     

Valence bond descriptions of benzene and cyclobutadiene and their counterparts with localized bonds

Geometry optimizations have been performed for benzene and cyclobutadiene and for the corresponding moieties with nonresonating double bonds, viz. 1,3,5‐cyclohexatriene and 1,3‐cyclobutadiene. The calculations were done using the valence bond self‐consistent field method including orbital optimization. Both strictly local and delocalized p‐like orbitals were used for the π system, which influences the strengths of the π bonds. The calculations result in geometries and resonance and stabilization energies for benzene and cyclobutadiene, which are compared with theoretical models of aromaticity. The importance of resonance is discussed. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Electronic structures of permethyloligosilane radical cations at the ground and low-lying excited states

The electronic structures at the ground and low-lying excited states of permethyloligosilane radical cations, Si_n(CH₃)_2n+2⁺ (n = 4-7), have been investigated using DFT and ab initio calculations. The calculations showed that positive charge (hole) is delocalized along the Si-Si main chain at the ground and first excited states. On the other hand, the hole is transferred to the methyl side-chain at the second and higher excited states. From these results, it was concluded that hole can move along the Si-Si main chain at thermal conditions. Also, it was predicted that intermolecular hole hopping takes place by photo-irradiation to the permethyloligosilane radical cation. The mechanism of hole transfer was discussed on the basis of the results.     

Molecular mechanics-valence bond method for planar conjugated hydrocarbon cations

We present an extension of the molecular mechanics-valence bond (MMVB) hybrid method to study ground and excited states of planar conjugated hydrocarbon cations. Currently, accurate excited state calculations on these systems are limited to expensive ab initio studies of smaller systems: up to 15 active electrons in 16 pi orbitals with complete active space self-consistent field (CASSCF) theory using high symmetry. The new MMVB extension provides a faster, cheaper treatment to investigate larger cation systems with more than 24 active orbitals. Extension requires both new matrix elements and new parameters: In this paper we present both, for the limited planar case. The scheme is tested for the planar radical cations of benzene, naphthalene, anthracene, and phenanthrene. Calculated MMVB relative energies are in good agreement with CASSCF results for equilibrium geometries on the ground and first excited states, and conical intersections.