Linear-scaling formation of Kohn-Sham Hamiltonian: application to the calculation of excitation energies and polarizabilities of large molecular systems

We present calculations of excitation energies and polarizabilities in large molecular systems at the local-density and generalized-gradient approximation levels of density-functional theory (DFT). Our results are obtained using a linear-scaling DFT implementation in the program system DALTON for the formation of the Kohn-Sham Hamiltonian. For the Coulomb contribution, we introduce a modification of the fast multipole method to calculations over Gaussian charge distributions. It affords a simpler implementation than the original continuous fast multipole method by partitioning the electrostatic Coulomb interactions into "classical" and "nonclassical" terms which are explicitly evaluated by linear-scaling multipole techniques and a modified two-electron integral code, respectively. As an illustration of the code, we have studied the singlet and triplet excitation energies as well as the static and dynamic polarizabilities of polyethylenes, polyenes, polyynes, and graphite sheets with an emphasis on the trends observed with system size.     

Excited state properties,fluorescence energies,and lifetimes of a poly(fluorene‐phenylene), based on TD‐DFT investigation

The structural and electronic properties of fluorene‐phenylene copolymer (FP)_n, n = 1–4 were studied by means of quantum chemical calculations based on density functional theory (DFT) and time dependent density functional theory (TD‐DFT) using B3LYP functional. Geometry optimizations of these oligomers were performed for the ground state and the lowest singlet excited state. It was found that (FP)_n is nonplanar in its ground state while the electronic excitations lead to planarity in its S₁ state. Absorption and fluorescence energies were calculated using TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods. Vertical excitation energies and fluorescence energies were obtained by extrapolating these values to infinite chain length, resulting in extrapolated values for vertical excitation energy of 2.89 and 2.87 eV, respectively. The S₁ ← S₀ electronic excitation is characterized as a highest occupied molecular orbital to lowest unoccupied molecular orbital transition and is distinguishing in terms of oscillator strength. Fluorescence energies of (FP)_n calculated from TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods are 2.27 and 2.26 eV, respectively. Radiative lifetimes are predicted to be 0.55 and 0.51 ns for TD‐B3LYP/SVP and TD‐B3LYP/SVP+ calculations, respectively. These fundamental information are valuable data in designing and making of promising materials for LED materials. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Adsorption‐induced changes of intramolecular optical transitions: PTCDA/NaCl and PTCDA/KCl

Structural and optical properties of isolated perylene‐3,4,9,10‐tetracarboxylic acid dianhydride molecules adsorbed on (100) oriented NaCl and KCl surfaces were studied theoretically to analyze the recently observed red‐shift of the optical excitation spectrum after adsorption (Müller et al., Phys. Rev. B, 2011, 83, 241203; Paulheim et al. Phys. Chem. Chem. Phys., 2013, 15, 4906). The ground‐state structures were obtained by periodic dispersion‐corrected density functional theory (DFT) calculations. For the excited‐state calculations, nonperiodic time‐dependent DFT methods were applied for a cluster model embedded in point charges. The range‐separated hybrid functional CAM‐B3LYP was used. Correlation‐consistent basis sets were used and the calculated excitation energies were extrapolated to the complete basis set limit. The shift of the first optical excitation energy was analyzed in terms of electronic and geometric contributions. It was found that both the distortion of the molecule due to the interaction with the surface and the electrostatic potential of the surface play an important role. © 2015 Wiley Periodicals, Inc.     

Highly efficient implementation of pseudospectral time‐dependent density‐functional theory for the calculation of excitation energies of large molecules

We have developed and implemented pseudospectral time‐dependent density‐functional theory (TDDFT) in the quantum mechanics package Jaguar to calculate restricted singlet and restricted triplet, as well as unrestricted excitation energies with either full linear response (FLR) or the Tamm–Dancoff approximation (TDA) with the pseudospectral length scales, pseudospectral atomic corrections, and pseudospectral multigrid strategy included in the implementations to improve the chemical accuracy and to speed the pseudospectral calculations. The calculations based on pseudospectral time‐dependent density‐functional theory with full linear response (PS‐FLR‐TDDFT) and within the Tamm–Dancoff approximation (PS‐TDA‐TDDFT) for G2 set molecules using B3LYP/6‐31G*^* show mean and maximum absolute deviations of 0.0015 eV and 0.0081 eV, 0.0007 eV and 0.0064 eV, 0.0004 eV and 0.0022 eV for restricted singlet excitation energies, restricted triplet excitation energies, and unrestricted excitation energies, respectively; compared with the results calculated from the conventional spectral method. The application of PS‐FLR‐TDDFT to OLED molecules and organic dyes, as well as the comparisons for results calculated from PS‐FLR‐TDDFT and best estimations demonstrate that the accuracy of both PS‐FLR‐TDDFT and PS‐TDA‐TDDFT. Calculations for a set of medium‐sized molecules, including C_n fullerenes and nanotubes, using the B3LYP functional and 6‐31G^** basis set show PS‐TDA‐TDDFT provides 19‐ to 34‐fold speedups for C_n fullerenes with 450–1470 basis functions, 11‐ to 32‐fold speedups for nanotubes with 660–3180 basis functions, and 9‐ to 16‐fold speedups for organic molecules with 540–1340 basis functions compared to fully analytic calculations without sacrificing chemical accuracy. The calculations on a set of larger molecules, including the antibiotic drug Ramoplanin, the 46‐residue crambin protein, fullerenes up to C₅₄₀ and nanotubes up to 14×(6,6), using the B3LYP functional and 6‐31G^** basis set with up to 8100 basis functions show that PS‐FLR‐TDDFT CPU time scales as N^2.05 with the number of basis functions. © 2016 Wiley Periodicals, Inc.     

New relativistic ANO basis sets for transition metal atoms

New basis sets of the atomic natural orbital (ANO) type have been developed for the first, second, and third row transition metal atoms. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive and negative ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies, electron affinities, and excitation energies for all atoms and polarizabilities for spherically symmetric atoms. These calculations include spin-orbit coupling using a variation-perturbation approach. Computed ionization energies have an accuracy better than 0.2 eV in most cases. The accuracy of computed electron affinities is the same except in cases where the experimental values are smaller than 0.5 eV. Accurate results are obtained for the polarizabilities of atoms with spherical symmetry. Multiplet levels are presented for some of the third row transition metals.     

Spectral/structural data of helium atoms with exponential‐cosine‐screened coulomb potentials

Atomic systems with exponential‐cosine‐screened Coulomb (ECSC) or static screened Coulomb (SSC) potentials have drawn considerable attention recently due to the possible applications to atoms in plasma environments. In this work, we develop a computing scheme to deal with the electron–electron correlation terms with the screened Coulomb interactions instead of the conventional expansion method using the Gegenbauer's addition theorem. Based on this approach, we investigate the helium atom with the ECSC potentials. Bypassing the complex expansion functions for electron–electron interactions, the proposed approach simplifies the calculations greatly and provides an advantage on programming. The results are found to be in good agreement with the existing data. Bound‐state energies, oscillator strengths, and multipole polarizabilities varying with the screening parameters are presented. Comparisons of the ECSC potentials are made with the SSC potentials. The influence of screening effect on the energies, oscillator strengths, and polarizabilities is discussed. © 2015 Wiley Periodicals, Inc.     

Substituent Effects in the Absorption Spectra of Phenol Radical Species: Origin of the Redshift Caused by 3,5‐Dimethoxyl Substitution

The ground‐state equilibrium geometries, electronic structures and vertical excitation energies of methyl‐ and methoxyl‐substituted phenol radical cations and phenoxyl radicals have been investigated using time‐dependent density‐functional theory (namely TD‐B3LYP) and complete‐active‐space second‐order perturbation theory (CASPT2). The “anomalous” large redshifts of the absorption maxima of the phenol radical species observed in the ultraviolet–visible spectral region upon di‐meta‐methoxyl substitution are reproduced by the calculations. Furthermore, these “anomalous” shifts which were unexplained to date can be rationalized on the basis of a qualitative molecular orbital perturbation analysis.     

Density functional self-consistent quantum mechanics/molecular mechanics theory for linear and nonlinear molecular properties: Applications to solvated water and formaldehyde

A combined quantum mechanics/molecular mechanics (QM/MM) method is described, where the polarization between the solvent and solute is accounted for using a self-consistent scheme linear in the solvent polarization. The QM/MM method is implemented for calculation of energies and molecular response properties including the calculation of linear and quadratic response functions using the density-functional theory (DFT) and the Hartree-Fock (HF) theory. Sample calculations presented for ground-state energies, first-order ground-state properties, excitation energies, first-order excited state properties, polarizabilities, first-hyperpolarizabilities, and two-photon absorptions strengths of formaldehyde suggests that DFT may in some cases be a sufficiently reliable alternative to high-level theory, such as coupled-cluster (CC) theory, in modeling solvent shifts, whereas results obtained with the HF wave function deviate significantly from the CC results. Calculations carried out on water gives results that also are comparable with CC calculations in accuracy for ground-state and first-order properties. However, to obtain such accuracy an exchange-correlation functional capable of describing the diffuse Rydberg states must be chosen.     

Computational linear dependence in molecular electronic structure calculations using universal basis sets

Distributed universal even‐tempered basis sets have been developed over recent years that are capable of supporting Hartree–Fock energies to an accuracy approaching the sub‐μHartree level. These basis sets have also been exploited in correlation studies, in applications to polyatomic molecules, and in the calculation of electric properties, such as multipole moments, polarizabilities, and hyperpolarizabilities. Jorge and coworkers have also developed universal basis sets and have recently reported applications to diatomic molecular systems. In this article, we compare the molecular calculations reported by Jorge and coworkers with our previous studies. Particular attention is given to the degree of computational linear dependence associated with the various basis sets employed and the consequential effects of the accuracy of the calculated energies. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Analytical time-dependent density functional derivative methods within the RI-J approximation, an approach to excited states of large molecules

Time-dependent density functional theory (TDDFT) is now well established as an efficient method for molecular excited state treatments. In this work, we introduce the resolution of the identity approximation for the Coulomb energy (RI-J) to excited state gradient calculations. In combination with nonhybrid functionals, the RI-J approximation leads to speed ups in total timings of an order of magnitude compared to the conventional method; this is demonstrated for oligothiophenes with up to 40 monomeric units and adamantane clusters. We assess the accuracy of the computed adiabatic excitation energies, excited state structures, and vibrational frequencies on a set of 36 excited states. The error introduced by the RI-J approximation is found to be negligible compared to deficiencies of standard basis sets and functionals. Auxiliary basis sets optimized for ground states are suitable for excited state calculations with small modifications. In conclusion, the RI-J approximation significantly extends the scope of applications of analytical TDDFT derivative methods in photophysics and photochemistry.     

Grid‐based density functional calculations of many‐electron systems

Exploratory variational pseudopotential density functional calculations are performed for the electronic properties of many‐electron systems in the 3D cartesian coordinate grid (CCG). The atom‐centered localized gaussian basis set, electronic density, and the two‐body potentials are set up in the 3D cubic box. The classical Hartree potential is calculated accurately and efficiently through a Fourier convolution technique. As a first step, simple local density functionals of homogeneous electron gas are used for the exchange‐correlation potential, while Hay‐Wadt‐type effective core potentials are employed to eliminate the core electrons. No auxiliary basis set is invoked. Preliminary illustrative calculations on total energies, individual energy components, eigenvalues, potential energy curves, ionization energies, and atomization energies of a set of 12 molecules show excellent agreement with the corresponding reference values of atom‐centered grid as well as the grid‐free calculation. Results for three atoms are also given. Combination of CCG and the convolution procedure used for classical Coulomb potential can provide reasonably accurate and reliable results for many‐electron systems. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

DFT/TDDFT Studies on the Electronic Structures and Spectral Properties of Carbazole‐based Blue Light‐emitting Dendrimers

The purpose of this paper is to provide an in‐depth investigation of the electronic and optical properties of two series of carbazole‐based blue light‐emitting dendrimers, including 1 – 6 six oligomers. These materials show great potential for application in organic light‐emitting diodes as efficient blue‐light and red‐light emitting materials due to the tuning of the optical and electronic properties by the use of different electron donors (D) and electron acceptors (A). The geometric and electronic structures of these compounds in the ground state are calculated using density functional theory (DFT) and the ab initio HF, whereas the lowest singlet excited states were optimized by ab initio single excitation configuration interaction (CIS). All DFT calculations are performed using the B3LYP functional on 6‐31G* basis set. The outcomes show that the highest occupied molecular orbitals (HOMOs), lowest occupied molecular orbitals (LUMOs), energies gaps, ionization potentials, electron affinities and reorganization energies of each molecular are affected by different D and A moieties and different substitute positions.     

ERKALE—A flexible program package for X‐ray properties of atoms and molecules

ERKALE is a novel software program for computing X‐ray properties, such as ground‐state electron momentum densities, Compton profiles, and core and valence electron excitation spectra of atoms and molecules. The program operates at Hartree–Fock or density‐functional level of theory and supports Gaussian basis sets of arbitrary angular momentum and a wide variety of exchange‐correlation functionals. ERKALE includes modern convergence accelerators such as Broyden and ADIIS and it is suitable for general use, as calculations with thousands of basis functions can routinely be performed on desktop computers. Furthermore, ERKALE is written in an object oriented manner, making the code easy to understand and to extend to new properties while being ideal also for teaching purposes. © 2012 Wiley Periodicals, Inc.     

Ab initio calculation of the electronic absorption of functionalized octahedral silsesquioxanes via time-dependent density functional theory with range-separated hybrid functionals

Recent advances in the ability to functionalize octahedral silsesquioxanes with different photoactive ligands, and thereby tune their optical properties, suggest that these molecules may serve as potential building blocks of light-harvesting, photovoltaic, and photonic devices. In this paper we report extensive ab initio calculations of the excitation energies underlying the absorption spectra of these systems. The calculations are based on density functional theory for the ground electronic state and time-dependent density functional theory for the excited electronic states. The ability of the commonly used B3LYP functional to reproduce the experimentally observed absorption excitation energies is compared to that of recently developed range-separated hybrid functionals. The importance of pairing the range-separated hybrid functionals with basis sets that include diffuse and polarization basis functions is demonstrated in the case of vinyl-functionalized silsesquioxanes. Absorptive excitation energies are then calculated and compared with experiment for octahedral silsesquioxanes functionalized with larger ligands. The tunability of optical properties is demonstrated by considering the effect on the excitation energies of functionalizing the ligands with electron-donating or -withdrawing groups.     

Linear-scaling implementation of molecular response theory in self-consistent field electronic-structure theory

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field theories for the calculation of frequency-dependent molecular response properties and excitation energies is presented, based on a nonredundant exponential parametrization of the one-electron density matrix in the atomic-orbital basis, avoiding the use of canonical orbitals. The response equations are solved iteratively, by an atomic-orbital subspace method equivalent to that of molecular-orbital theory. Important features of the subspace method are the use of paired trial vectors (to preserve the algebraic structure of the response equations), a nondiagonal preconditioner (for rapid convergence), and the generation of good initial guesses (for robust solution). As a result, the performance of the iterative method is the same as in canonical molecular-orbital theory, with five to ten iterations needed for convergence. As in traditional direct Hartree-Fock and Kohn-Sham theories, the calculations are dominated by the construction of the effective Fock/Kohn-Sham matrix, once in each iteration. Linear complexity is achieved by using sparse-matrix algebra, as illustrated in calculations of excitation energies and frequency-dependent polarizabilities of polyalanine peptides containing up to 1400 atoms.     

Mechanism for large first hyperpolarizabilities of phosphonic acid stilbene derivatives

This paper presents calculations of dipole moments (mu), static polarizabilities (alpha), and first hyperpolarizabilities (beta) of phosphonic acid stilbene derivatives calculated in the framework of density functional theory. These calculations were performed using a finite field approach implemented in the density functional program ALLCHEM and were of an all-electron type using local exchange-correlation functional and specially designed basis sets. The molecular structures have been fully optimized using the semiempirical program MSINDO. Some of the investigated stilbenes have been synthesized very recently while others are described for the first time. Donor and acceptor groups of these analogues have been modified and the influence of these changes on the first hyperpolarizabilities has been investigated. This work demonstrates that the nonlinear optical response beta of these compounds increases dramatically when the acceptor moiety is displaced by analogues containing alkali metal groups. A general mechanism for the design of novel nonlinear optical materials with large first hyperpolarizabilities is described.     

Singlet excitation energies of thiophene derivatives of fluorene: TD‐DFT study

Fluorene‐thiophene (FT)‐based oligomers and polymers and their derivatives are good candidates for organic blue light‐emitting diodes. In this work, the intrinsic properties of the ground and excited states of FT monomer and its derivatives are studied. The ground‐state optimized structures and energies are obtained using molecular orbital theory and density functional theory (DFT). The ground‐state potential energy curves or surfaces of FT and its derivatives are also obtained. All derivatives are nonplanar in their electronic ground states. The character and energy of the first 20 singlet–singlet electronic transitions are investigated by applying the time‐dependent density functional theory (TD‐DFT) approximations to the correspondingly optimized ground‐state geometries. The lowest singlet state is studied with the configuration interaction (singles) approach (CIS). Excitation energies are red shifted when the FT unit or its derivatives are extended longitudinally. CIS results suggest geometry relaxation in the first singlet excited state. When available, a comparison is made with experimental results. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007