Theoretical study of excited states of DNA base dimers and tetramers using optimally tuned range‐separated density functional theory

Excited states of various DNA base dimers and tetramers including Watson‐Crick H‐bonding and stacking interactions have been investigated by time‐dependent density functional theory using nonempirically tuned range‐separated exchange (RSE) functionals. Significant improvements are found in the prediction of excitation energies and oscillator strengths, with results comparable to those of high‐level coupled‐cluster (CC) models (RI‐CC2 and EOM‐CCSD(T)). The optimally‐tuned RSE functional significantly outperforms its non‐tuned (default) version and widely‐used B3LYP functional. Compared to those high‐level CC benchmarks, the large mean absolute deviations of conventional functionals can be attributed to their inappropriate amount of exact exchange and large delocalization errors which can be greatly eliminated by tuning approach. Furthermore, the impacts of H‐bonding and π‐stacking interactions in various DNA dimers and tetramers are analyzed through peak shift of simulated absorption spectra as well as corresponding change of absorption intensity. The result indicates the stacking interaction in DNA tetramers mainly contributes to the hypochromicity effect. The present work provides an efficient theoretical tool for accurate prediction of optical properties and excited states of nucleobase and other biological systems. © 2015 Wiley Periodicals, Inc.     

Assessment of range‐separated exchange functionals and nonempirical functional tuning for calculating the static second hyperpolarizabilities of streptocyanines

DFT method can severely overestimate the response properties for π‐conjugation systems. The range‐separated exchange and recently developed optimal IP‐tuning process are evaluated on the prediction of static second hyperpolarizabilities of streptocyanines of increasing molecular length. The finite field results have shown that the exact exchange at midium and long distance can relieve only a part of the overshooting but still beyond satisfaction. The exact exchange at short distance has the oppsite effects showing the failure of converntional hGGA. The optimal tuned range‐separated exchange functionals show little improvements performing worse than the default ones. Importantly, the electronic structure–property relationship, bond order alternation‐γ , is not well established with DFT method. © 2017 Wiley Periodicals, Inc.     

Performance of range‐separated hybrid exchange–correlation functionals for the calculation of magnetic exchange coupling constants of organic diradicals

The prediction of magnetic behavior is important for the design of new magnetic materials. Kohn–Sham density functional theory is popular for this purpose, although one should be careful about choosing the right exchange–correlation functional. Here, we perform a statistical analysis to test different range‐separated hybrid density functionals for the calculation of magnetic exchange coupling constants J of fourteen organic diradicals. Our analysis suggests that in absolute terms the MN12SX functional performs best among the series of twelve functionals studied here (including the popular B3LYP), followed by N12SX functionals along with Scuseria's HSE series of functionals. LC‐ PBE was found to be the least accurate, which is in contrast with its good performance for calculating J for transition metal complexes. The HSE family of functionals and B3LYP are the only functionals to reproduce the qualitative trends of the coupling constants correctly for the ferromagnetically coupled diradicals under study. © 2017 Wiley Periodicals, Inc.     

Implementation of nuclear gradients of range‐separated hybrid density functionals and benchmarking on rotational constants for organic molecules

We have implemented the nuclear gradient for several range‐separated hybrid density functionals in the general quantum chemistry code ORCA. To benchmark the performance, we have used a recently published set of back‐corrected gas phase rotational constants, which we extended by three molecules. In our evaluation, CAM‐B3LYP‐D3 and ωB97X‐D3 show great accuracy, and are surpassed by B2PLYP‐D3 only. Lower‐cost alternatives to quadruple‐ζ basis set‐based calculations, among them a smaller basis set and the use of resolution‐of‐the‐identity approaches, are assessed and shown to yield acceptable deviations. In addition, the Hartree‐Fock‐based back‐correction method is compared to a density functional theory alternative, which largely shows consistency between the two. A new, well‐performing, spin‐component scaled MP2 variant is designed and discussed, as well. © 2014 Wiley Periodicals, Inc.     

Prediction of excited‐state properties of oligoacene crystals using polarizable continuum model‐tuned range‐separated hybrid functional approach

A methodology combining the polarizable continuum model and optimally‐tuned range‐separated (RS) hybrid functional was proposed for the quantitative characterization of the excited‐state properties in oligoacene (from anthracene to hexacene) crystals. We show that it provides lowest vertical singlet and triplet excitation energies, singlet‐triplet gap, and exciton binding energies in very good agreement with the available experimental data. We further find that it significantly outperforms its non‐tuned RS counterpart and the widely used B3LYP functional, and even many‐body perturbation theory within the GW approximation (based on a PBE starting point). Hence, this approach provides an easily applicable and computationally efficient tool to study the excited‐state properties of organic solids of complexity. © 2017 Wiley Periodicals, Inc.     

The self‐consistent charge density functional tight‐binding theory study of carbon adatoms using tuned Hubbard U parameters

The self‐consistent charge density functional tight‐binding (DFTB) theory is a useful tool for realizing the electronic structures of large molecular complex systems. In this study, the electronic structure of C₆₁ formed by fullerene C₆₀ with a carbon adatom is analyzed, using the fully localized limit and pseudo self‐interaction correction methods of DFTB to adjust the Hubbard U parameter (DFTB + U). The results show that both the methods used to adjust U can significantly reduce the molecular orbital energy of occupied states localized on the defect carbon atom and improve the gap between highest occupied molecular orbital(HOMO) and lowest unoccupied molecular orbital(LUMO) of C₆₁. This work will provide a methodological reference point for future DFTB calculations of the electronic structures of carbon materials.     

Reaction energetics on long‐range corrected density functional theory: Diels–Alder reactions

The possibility of quantitative reaction analysis on the orbital energies of long‐range corrected density functional theory (LC‐DFT) is presented. First, we calculated the Diels–Alder reaction enthalpies that have been poorly given by conventional functionals including B3LYP functional. As a result, it is found that the long‐range correction drastically improves the reaction enthalpies. The barrier height energies were also computed for these reactions. Consequently, we found that dispersion correlation correction is also crucial to give accurate barrier height energies. It is, therefore, concluded that both long‐range exchange interactions and dispersion correlations are essentially required in conventional functionals to investigate Diels–Alder reactions quantitatively. After confirming that LC‐DFT accurately reproduces the orbital energies of the reactant and product molecules of the Diels–Alder reactions, the global hardness responses, the halves of highest occupied molecular orbital (HOMO)‐lowest unoccupied molecular orbital (LUMO) energy gaps, along the intrinsic reaction coordinates of two Diels–Alder reactions were computed. We noticed that LC‐DFT results satisfy the maximum hardness rule for overall reaction paths while conventional functionals violate this rule on the reaction pathways. Furthermore, our results also show that the HOMO‐LUMO gap variations are close to the reaction enthalpies for these Diels–Alder reactions. Based on these results, we foresee quantitative reaction analysis on the orbital energies. © 2012 Wiley Periodicals, Inc.     

A simple approximation for the Pauli potential yielding self‐consistent electron densities exhibiting proper atomic shell structure

A simple approximation for the Pauli potential for the groundstate of atomic systems is given, which in connection with Hohenberg–Kohn variational procedure yields self‐consistent electron densities exhibiting proper atomic shell structure. © 2015 Wiley Periodicals, Inc.     

Algebraic modifications to second quantization for non‐Hermitian complex scaled hamiltonians with application to a quadratically convergent multiconfigurational self‐consistent field method

The algebraic structure for creation and annihilation operators defined on orthogonal orbitals is generalized to permit easy development of bound‐state techniques involving the use of non‐Hermitian Hamiltonians arising from the use of complex‐scaling or complex‐absorbing potentials in the treatment of electron scattering resonances. These extensions are made possible by an orthogonal transformation of complex biorthogonal orbitals and states as opposed to the customary unitary transformation of real orthogonal orbitals and states and preserve all other formal and numerical simplicities of existing bound‐state methods. The ease of application is demonstrated by deriving the modified equations for implementation of a quadratically convergent multiconfigurational self‐consistent field (MCSCF) method for complex‐scaled Hamiltonians but the generalizations are equally applicable for the extension of other techniques such as single and multireference coupled cluster (CC) and many‐body perturbation theory (MBPT) methods for their use in the treatment of resonances. This extends the domain of applicability of MCSCF, CC, MBPT, and methods based on MCSCF states to an accurate treatment of resonances while still using L² real basis sets. Modification of all other bound‐state methods and codes should be similarly straightforward. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

About the compatibility between ansatzes and constraints for a local formulation of orbital‐free density functional theory

Functional properties that are exact for the Hohenberg–Kohn functional may turn into mutually exclusive constraints at a given level of ansatz. This is exemplarily shown for the local density approximation. Nevertheless, it is possible to reach exactly the Kohn–Sham data from an orbital‐free density functional framework based on simple one‐point functionals by starting from the Levy–Perdew–Sahni formulation. The energy value is obtained from the density‐potential pair, and therefore does not refer to the functional dependence of the potential expression. Consequently, the potential expression can be obtained from any suitable model and is not required to follow proper scaling behavior.     

Toward assessment of density functionals for vibronic coupling in two‐photon absorption: A case study of 4‐nitroaniline

In this study, we predict vibronic two‐photon absorption (TPA) spectra for 4‐nitroaniline in vacuo. The simulations are performed using density functional theory and the approximate second‐order coupled‐cluster singles and doubles model CC2. Thereby we also demonstrate the possibility of simulations of vibronic TPA spectra with ab initio wavefunction methods that include electron correlation for medium‐sized systems. A special focus is put on the geometric derivatives of the second‐order transition moment and the dipole moment difference between the charge‐transfer excited state and the ground state. The results of CC2 calculations bring new insight into the vibronic coupling mechanism in TPA spectra of 4‐nitroniline and demonstrate that the mixed term is quite large and that it also exhibits a negative interference with the Franck‐Condon contribution. © 2015 Wiley Periodicals, Inc.     

New procedure to evaluate aromaticity at the density functional theory,Hartree–Fock,and post‐self‐consistent field levels

The spatial exchange interaction, arising from the exchange‐type two‐electron integrals ( ) between two different groups P and Q, is another driving force for the delocalization of π‐electrons besides orbital charge‐transfer and exchange interactions. We have developed a new combination program for restricted geometry optimization, in which all of the orbital and spatial interactions among isolated groups were excluded from the localized geometry of a conjugated molecule. This was achieved by deleting particular Fock elements and the 15 types of exchange‐type two‐electron integrals, ensuring that the corresponding π‐electrons are completely localized within their respective groups and the π‐orbitals are fully localized. The extra stabilization energy (ESE) of benzene is ?36.3 kcal/mol (B3LYP/6‐31G*), and the level of density functional theory, Hartree–Fock, and post‐self‐consistent field (Møller–Plesset 2, configuration interaction singles and doubles, and singles and doubles coupled‐cluster) and the basis sets have slight effect on the ESE. Based on the comparisons between our procedure, Morokuma's energy decomposition analysis and the block‐localized wave function method, it was confirmed that our program calculates reliable results. The nonaromaticity of acyclic polyenes and antiaromaticity of cyclobutadiene and planar cyclooctatetraene were also estimated. Comparison of the C? C single bond lengths in the ground state with its π‐localized geometries showed that shortening of the single bonds in acyclic polyenes and butadiyne should be attributed to different hybridization, demonstrating that the effect of π‐delocalization on single bonds is so small as to be negligible. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Density functional theory study of calix[4]arene‐N‐azacrown‐5, calix[4]arene‐N‐phenyl‐azacrown‐5, and their complexes with alkali‐metal cations: Na+, K+, and Rb+

Theoretical studies of 1,3‐alternate‐25,27‐bis(1‐methoxyethyl)calix[4]arene‐azacrown‐5 ( L₁ ), 1,3‐alternate‐25,27‐bis(1‐methoxyethyl)calix[4]arene‐N‐phenyl‐azacrown‐5 ( L₂ ), and the corresponding complexes M⁺/ L of L₁ and L₂ with the alkali‐metal cations: Na⁺, K⁺, and Rb⁺ have been performed using density functional theory (DFT) at B3LYP/6‐31G* level. The optimized geometric structures obtained from DFT calculations are used to perform natural bond orbital (NBO) analysis. The two main types of driving force metal–ligand and cation–π interactions are investigated. The results indicate that intermolecular electrostatic interactions are dominant and the electron‐donating oxygen offer lone pair electrons to the contacting RY* (1‐center Rydberg) or LP* (1‐center valence antibond lone pair) orbitals of M⁺ (Na⁺, K⁺, and Rb⁺). What's more, the cation–π interactions between the metal ion and π‐orbitals of the two rotated benzene rings play a minor role. For all the structures, the most pronounced changes in geometric parameters upon interaction are observed in the calix[4]arene molecule. In addition, an extra pendant phenyl group attached to nitrogen can promote metal complexation by 3D encapsulation greatly. In addition, the enthalpies of complexation reaction and hydrated cation exchange reaction had been studied by the calculated thermodynamic data. The calculated results of hydrated cation exchange reaction are in a good agreement with the experimental data for the complexes. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Electron smearing in DFT calculations: A case study of doxorubicin interaction with single‐walled carbon nanotubes

To address the choice of an appropriate value of electron smearing to facilitate self‐consistent field (SCF) convergence, we studied the interaction of doxorubicin with short armchair and zigzag single‐walled carbon nanotube models with closed caps, at the PWC/DNP level of density functional theory. By gradually reducing the electron smearing value from a large and most commonly used one of 0.005 Ha to zero (Fermi occupation), we monitored the changes in close contacts between the interacting species, total energy of the molecular system, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy and isosurfaces, HOMO‐LUMO gap energy, and plots of electrostatic potential. It became evident that the commonly used smearing values of ≥0.001 Ha can alter the results significantly (for example, by one order of magnitude for HOMO–LUMO gap energy). We suggest the setting of electron smearing value at 0.0001 Ha, which does not imply too high computation cost and can guarantee the results close to the ones obtained with Fermi occupation. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

A complete set of self‐consistent charge density‐functional tight‐binding parametrization of zinc chalcogenides (ZnX; X=O,S, Se,and Te)

We have developed a complete set of self‐consistent charge density‐functional tight‐binding parameters for Zn? X (X = Zn, O, S, Se, Te, Cd, H, C, and N). The transferability of the derived parameters has been tested against Pseudo Potential‐Perdew, Burke and Ernzerhof (PP‐PBE) calculations and experimental values (whenever available) for corresponding bulk systems (e.g., hexagonal close packing, zinc‐blende, and wurtzite(wz)), various kinds of nanostructures (such as nanowires, surfaces, and nanoclusters), and also some small molecular systems. Our results show that the derived parameters reproduce the structural and energetic properties of the above‐mentioned systems very well. With the derived parameter set, one can study zinc‐chalcogenide nanostructures of relatively large size which was otherwise prohibited by other methods. The Zn‐Cd parametrization developed in this article will help in studying large semiconductor hetero‐nanostructures of Zn and Cd chalcogenides such as ZnX/CdX core/shell nanoparticles, nanotubes, nanowires, and nanoalloys. © 2012 Wiley Periodicals, Inc.     

Electronic Structure of Thiolate‐bridged Diiron Complexes and a Single‐electron Oxidation Reaction: A Combination of Experimental and Computational Studies

Single‐electron oxidation of a diiron‐sulfur complex [Cp*Fe(μ‐bdt)FeCp*] ( 1 , Cp*=η⁵‐C₅Me₅; bdt=benzene‐1,2‐dithiolate) to [Cp*Fe(μ‐bdt)FeCp*]⁺ ( 2 ) has been experimentally conducted. The bdt ligand with redox‐active character has been computationally proposed to be a dianion (bdt^2?) rather than previously proposed monoanion (bdt^·?) radical in 1 though it has un‐equidistant aromatic C? C bond lengths. The ground state of 1 is predicted to be two low‐spin ferrous ions (S_Fe=0) and 2 has a medium‐spin ferric ion (S_Fe=1/2) and a low‐spin ferrous center (S_Fe=0), and the oxidation of 1 to 2 is calculated to be a single‐metal‐based process. Both complexes have no significant antiferromagnetic coupling character.     

Kinetics of radical‐molecule reactions in aqueous solution: A benchmark study of the performance of density functional methods

The performance of 18 density functional approximations has been tested for a very challenging task, the calculations of rate constants for radical‐molecule reactions in aqueous solution. Despite of the many difficulties involved in such an enterprise, six of them provide high quality results, and are recommended to that purpose. They are LC‐ωPBE, M06‐2X, BMK, B2PLYP, M05‐2X, and MN12SX, in that order. This trend was obtained using experimental data as reference. The other relevant aspects used in this benchmark are: (i) the SMD model for mimicking the solvent; (ii) the conventional transition state, the zero‐curvature tunneling correction, and the limit imposed by diffusion for the calculation of the rate constants. Even though changing any of these aspects might alter the trend in performance, at least, when using them, the aforementioned functionals can be successfully used to obtain high quality kinetic data for the kind of reactions investigated in this work. © 2014 Wiley Periodicals, Inc.