Full spectral decomposition of ring currents

Within the ipsocentric method for calculation of molecular magnetic response, projection of perturbed orbitals onto the virtual orbital space allows partition of induced current density into contributions from individual virtual excitations between occupied and unoccupied orbitals, enabling detailed assignment of the origin of currents in, e.g., benzene, cyclooctatetraene, borazine, coronene, and corannulene. Whereas delocalized currents in benzene and planar cyclooctatetraene are described by transitions within the valence space, localized currents in the borazine pi system involve excitations outside the valence space.     

Asymptotic correction of the exchange-correlation kernel of time-dependent density functional theory for long-range charge-transfer excitations

Time-dependent density functional theory (TDDFT) calculations of charge-transfer excitation energies omegaCT are significantly in error when the adiabatic local density approximation (ALDA) is employed for the exchange-correlation kernel fxc. We relate the error to the physical meaning of the orbital energy of the Kohn-Sham lowest unoccupied molecular orbital (LUMO). The LUMO orbital energy in Kohn-Sham DFT--in contrast to the Hartree-Fock model--approximates an excited electron, which is correct for excitations in compact molecules. In CT transitions the energy of the LUMO of the acceptor molecule should instead describe an added electron, i.e., approximate the electron affinity. To obtain a contribution that compensates for the difference, a specific divergence of fxc is required in rigorous TDDFT, and a suitable asymptotically correct form of the kernel fxc(asymp) is proposed. The importance of the asymptotic correction of fxc is demonstrated with the calculation of omegaCT(R) for the prototype diatomic system HeBe at various separations R(He-Be). The TDDFT-ALDA curve omegaCT(R) roughly resembles the benchmark ab initio curve omegaCT CISD(R) of a configuration interaction calculation with single and double excitations in the region R=1-1.5 A, where a sizable He-Be interaction exists, but exhibits the wrong behavior omegaCT(R)相似文献   

Away from generalized gradient approximation: orbital-dependent exchange-correlation functionals

The local-density approximation of density functional theory (DFT) is remarkably accurate, for instance, for geometries and frequencies, and the generalized gradient approximations have also made bond energies quite reliable. Sometimes, however, one meets with failure in individual cases. One of the possible routes towards better functionals would be the incorporation of orbital dependence (which is an implicit density dependency) in the functionals. We discuss this approach both for energies and for response properties. One possibility is the use of the Hartree-Fock-type exchange energy expression as orbital-dependent functional. We will argue that in spite of the increasing popularity of this approach, it does not offer any advantage over Hartree-Fock for energies. We will advocate not to apply the separation of exchange and correlation, which is so ingrained in quantum chemistry, but to model both simultaneously. For response properties the energies and shapes of the virtual orbitals are crucial. We will discuss the benefits that Kohn-Sham potentials can offer which are derived from either an orbital-dependent energy functional, including the exact-exchange functional, or which can be obtained directly as orbital-dependent functional. We highlight the similarity of the Hartree-Fock and Kohn-Sham occupied orbitals and orbital energies, and the essentially different meanings the virtual orbitals and orbital energies have in these two models. We will show that these differences are beneficial for DFT in the case of localized excitations (in a small molecule or in a fragment), but are detrimental for charge-transfer excitations. Again, orbital dependency, in this case in the exchange-correlation kernel, offers a solution.     

Optimized effective potential method: is it possible to obtain an accurate representation of the response function for finite orbital basis sets?

The optimized effective potential (OEP) equations are solved in a matrix representation using the orbital products of occupied and virtual orbitals for the representation of both the local potential and the response function. This results in a direct relationship between the matrix elements of local and nonlocal operators for the exchange-correlation potential. The effect of the truncation of the number of such products in the case of finite orbital basis sets on the OEP orbital and total energies and on the spectrum of eigenvalues of the response function is examined. Test calculations for Ar and Ne show that rather large AO basis sets are needed to obtain an accurate representation of the response function.     

Response calculations based on an independent particle system with the exact one-particle density matrix: excitation energies

Adiabatic response time-dependent density functional theory (TDDFT) suffers from the restriction to basically an occupied → virtual single excitation formulation. Adiabatic time-dependent density matrix functional theory allows to break away from this restriction. Problematic excitations for TDDFT, viz. bonding-antibonding, double, charge transfer, and higher excitations, are calculated along the bond-dissociation coordinate of the prototype molecules H(2) and HeH(+) using the recently developed adiabatic linear response phase-including (PI) natural orbital theory (PINO). The possibility to systematically increase the scope of the calculation from excitations out of (strongly) occupied into weakly occupied ("virtual") natural orbitals to larger ranges of excitations is explored. The quality of the PINO response calculations is already much improved over TDDFT even when the severest restriction is made, to virtually the size of the TDDFT diagonalization problem (only single excitation out of occupied orbitals plus all diagonal doubles). Further marked improvement is obtained with moderate extension to allow for excitation out of the lumo and lumo+1, which become fractionally occupied in particular at longer distances due to left-right correlation effects. In the second place the interpretation of density matrix response calculations is elucidated. The one-particle reduced density matrix response for an excitation is related to the transition density matrix to the corresponding excited state. The interpretation of the transition density matrix in terms of the familiar excitation character (single excitations, double excitations of various types, etc.) is detailed. The adiabatic PINO theory is shown to successfully resolve the problematic cases of adiabatic TDDFT when it uses a proper PI orbital functional such as the PILS functional.     

The performance and relationship among range-separated schemes for density functional theory

The performance and relationship among different range-separated (RS) hybrid functional schemes are examined using the Coulomb-attenuating method (CAM) with different values for the fractions of exact Hartree-Fock (HF) exchange (α), long-range HF (β), and a range-separation parameter (μ), where the cases of α + β = 1 and α + β = 0 were designated as CA and CA0, respectively. Attenuated PBE exchange-correlation functionals with α = 0.20 and μ = 0.20 (CA-PBE) and α = 0.25 and μ = 0.11 (CA0-PBE) are closely related to the LRC-ωPBEh and HSE functionals, respectively. Time-dependent density functional theory calculations were carried out for a number of classes of molecules with varying degrees of charge-transfer (CT) character to provide an assessment of the accuracy of excitation energies from the CA functionals and a number of other functionals with different exchange hole models. Functionals that provided reasonable estimates for local and short-range CT transitions were found to give large errors for long-range CT excitations. In contrast, functionals that afforded accurate long-range CT excitation energies significantly overestimated energies for short-range CT and local transitions. The effects of exchange hole models and parameters developed for RS functionals for CT excitations were analyzed in detail. The comparative analysis across compound classes provides a useful benchmark for CT excitations.     

Long-range excitations in time-dependent density functional theory

Adiabatic time-dependent density functional theory fails for excitations of a heteroatomic molecule composed of two open-shell fragments at large separation. Strong frequency dependence of the exchange-correlation kernel is necessary for both local and charge-transfer excitations. The root of this is the static correlation created by the step in the exact Kohn-Sham ground-state potential between the two fragments. An approximate nonempirical kernel is derived for excited molecular dissociation curves at large separation. Our result is also relevant when the usual local and semilocal approximations are used for the ground-state potential, as static correlation there arises from the coalescence of the highest occupied and lowest unoccupied orbital energies as the molecule dissociates.     

Local orbitals by minimizing powers of the orbital variance

It is demonstrated that a set of local orthonormal Hartree-Fock (HF) molecular orbitals can be obtained for both the occupied and virtual orbital spaces by minimizing powers of the orbital variance using the trust-region algorithm. For a power exponent equal to one, the Boys localization function is obtained. For increasing power exponents, the penalty for delocalized orbitals is increased and smaller maximum orbital spreads are encountered. Calculations on superbenzene, C(60), and a fragment of the titin protein show that for a power exponent equal to one, delocalized outlier orbitals may be encountered. These disappear when the exponent is larger than one. For a small penalty, the occupied orbitals are more local than the virtual ones. When the penalty is increased, the locality of the occupied and virtual orbitals becomes similar. In fact, when increasing the cardinal number for Dunning's correlation consistent basis sets, it is seen that for larger penalties, the virtual orbitals become more local than the occupied ones. We also show that the local virtual HF orbitals are significantly more local than the redundant projected atomic orbitals, which often have been used to span the virtual orbital space in local correlated wave function calculations. Our local molecular orbitals thus appear to be a good candidate for local correlation methods.     

Calculations of the optical spectra of hydrocarbon radical cations based on Koopmans' theorem

The first few bands in the optical spectra of radical cations can often be interpreted in terms of A-type transitions that involve electron promotions from doubly occupied to the singly occupied molecular orbital (SOMO) and/or B-type transition which involve electron promotion from the SOMO to virtual molecular orbitals. We had previously demonstrated that, by making use of Koopmans' theorem, the energies of A-type transitions can be related to orbital energy differences between lower occupied MOs and the highest occupied MO (HOMO) in the neutral molecule, calculated at the geometry of the radical cation. We now propose that the energies of B-type transitions can be related similarly to energy differences between the lowest unoccupied MO (LUMO) and higher virtual MOs in the dication, also calculated at the geometry of the radical cation, by way of an extension of Koopmans' theorem to virtual MOs similar to that used sometimes to model resonances in electron scattering experiments. The optical spectra of the radical cations of several polyenes and aromatic compounds, the matrix spectra of which are known (or presented here for the first time), and for which CASSCF/CASPT2 calculations are available, are discussed in terms of these Koopmans-based models. Then the spectra of five poly(bicycloalkyl)-protected systems and that of hexabenzocoronene, compounds not amenable to higher level calculations, are examined and it is found that the Koopmans-type calculations allow a satisfactory interpretation of most of the features in these spectra. These simple calculations therefore provide a computationally inexpensive yet effective way to assign optical transitions in radical ions. Limitations of the model are discussed.     

Computational developments in generalized valence bond calculations

In this article a procedure for generating starting orbitals for generalized valence bond (GVB) calculations is presented. This is achieved by selecting orbitals which correspond to specific bonds or electron pairs. These orbitals can be identified from the localized molecular orbitals, for both occupied and virtual orbitals, which are obtained through a unitary transformation of the Hartree-Fock canonical molecular orbitals using the Boys's localization method. A scheme has also been implemented which achieves optimum convergence of the pairwise orbital optimization. An object-oriented GVB program is developed which automatically generates reliable initial GVB orbitals, leading to proper and fast convergence. © 1996 by John Wiley & Sons, Inc.     

Fast calculation of excited-state potentials for rare-gas diatomic molecules: Ne2 and Ar2

A fast method for obtaining excited-state potentials of rare-gas diatomic molecules is described. Two types of excited orbitals are used: molecular orbitals calculated in the field of a singly charged molecular ion, and atomic orbitals (properly symmetrized) obtained in a similar atomic system. The RPA equations are solved within the manifold of excitations from the highest occupied orbital in each symmetry to the lowest excited orbital of either type in each symmetry. A simple model for estimating the dynamic correlation correction to excitation (and ionization) energies is given. Applications to excited states of Ne₂(^1,3Σ⁺_{g, u}, ^1,3Π_{g, u}) and Ar₂(^1,3Σ⁺_{g, u}) are described. Two-electron integral transformations involve only three orbitals of each symmetry, and the RPA matrices are four-dimensional. The computational effort required for all excited-state potentials adds less than one-tenths (in terms of computer time) to the effort involved in the preliminary ground state Hartree—Fock calculations. The resulting potentials compare favorably with more elaborate CI calculations and give good agreement with spectroscopic and scattering data. Potential curves for the molecular ions are also given.     

A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.     

A local second-order Møller-Plesset method with localized orbitals: a parallelized efficient electron correlation method

Using orthogonal localized occupied orbitals we have developed and implemented a parallelized local second-order M?ller-Plesset (MP2) method based on the idea developed by Head-Gordon and co-workers. A subset of nonorthogonal correlation functions (the orbital domain) was assigned to each of the localized occupied orbitals using a distance criterion and excitations from localized occupied orbitals that were arranged into subsets. The correlation energy was estimated using a partial diagonalization and an iterative efficient method for solving large-scale linear equations. Some illustrative calculations are provided for molecules with up to 1484 Cartesian basis sets. The orbital domain sizes were found to be independent of the molecular size, and the present local MP2 method covered about 98%-99% of the correlation energy of the conventional canonical MP2 method.     

Theoretical study of the electronic spectroscopy of CO adsorbed on Pt(111)

The excited states of CO adsorbed on the Pt(111) surface are studied using a time-dependent density functional theory formalism. To reduce the computational cost, electronic excitations are computed within a reduced single excitation space. Using cluster models of the surface, excitation energies are computed for CO in the on-top, threefold, and bridge binding sites. On adsorption, there is a lowering of the 5sigma orbital energy. This leads to a large blueshift in the 5sigma- -> pi(CO*) excitation energy for all adsorption sites. The 1pi and 4sigma orbital energies are lowered to a lesser extent, and smaller shifts in the corresponding excitation energies are predicted. For the larger clusters, pi* excitations at lower energies are observed. These transitions correspond to excitations to virtual orbitals of pi* character which lie below the pi* orbitals of gas phase CO. These orbitals are associated predominantly with the metal atoms of the cluster. The excitation energies are also found to be sensitive to changes in the adsorption geometry. The electronic spectrum of CO on Pt(111) is simulated and the assignment of the bands observed in experimental electron energy loss spectroscopy discussed.     

An orbital-invariant and strictly size extensive post-Hartree-Fock correlation functional

A strictly size extensive post-Hartree-Fock correlation functional being invariant with respect to orbital transformations within the occupied and virtual subspaces is presented. While avoiding the necessity to solve additional Z vector equations for the calculation of properties and energy gradients, this functional reproduces almost exactly the results of coupled-cluster singles doubles (CCSD) calculations. In particular, it is demonstrated that the method is rigorous in the sense that it can be systematically improved by the perturbative inclusion of triple excitations in the same way as CCSD. As to the computational cost, the presented approach is somewhat more expensive than the CCSD if the energy is variationally optimized with respect to both the orbitals and the excitation amplitudes. Replacement of orbital optimization by the Brueckner condition reduces the computational cost by a factor of two, thus making the method less expensive than CCSD.     

Very-low-lying electronic states result from n --> pi excitations in open-shell annulenes, annelated with alpha-dicarbonyl groups

UB3LYP/6-31G* calculations find that alpha-dicarbonyl-annelated cyclopentadienyl radical 1 has a sigma ground state, which is formed by excitation of an electron from the in-phase combination of carbonyl lone-pair orbitals into the singly occupied pi orbital. Similarly, tetrakis-annelated cyclooctatetraene 3 is calculated to have very-low-lying singlet and triplet excited states, which result from excitations of electrons from the b1g combination of lone pair orbitals into the empty pi nonbonding MO of the COT ring.     

Linear scaling local correlation approach for solving the coupled cluster equations of large systems

A linear scaling local correlation approach is proposed for approximately solving the coupled cluster doubles (CCD) equations of large systems in a basis of orthogonal localized molecular orbitals (LMOs). By restricting double excitations from spatially close occupied LMOs into their associated virtual LMOs, the number of significant excitation amplitudes scales only linearly with molecular size in large molecules. Significant amplitudes are obtained to a very good approximation by solving the CCD equations of various subsystems, each of which is made up of a cluster associated with the orbital indices of a subset of significant amplitudes and the local environmental domain of the cluster. The combined effect of these two approximations leads to a linear scaling algorithm for large systems. By using typical thresholds, which are designed to target an energy accuracy, our numerical calculations for a wide range of molecules using the 6-31G or 6-31G* basis set demonstrate that the present local correlation approach recovers more than 98.5% of the conventional CCD correlation energy.     

Large-scale RPA calculations of chiroptical properties of organic molecules: Program RPAC

The computational considerations involved in calculating ordinary and rotatory intensities and electronic excitation energies in the random phase approximation (RPA ) are examined. We employ a localized orbital formulation in order to analyze the results in terms of local and charge-transfer excitations. Occupied orbitals are localized by the Foster–Boys procedure. The virtual space is transformed into a localized “valence” set that maximizes dipole strengths with the occupied counterparts, and a delocalized remainder. The two-electron integral transformation is performed with an efficient algorithm, based on Diercksen's, that generates only the particle–hole-type integrals required in the RPA . The lowest solutions of the RPA equations are obtained iteratively using a modification of the Davidson-Liu simultaneous vector expansion method. This allows the inclusion of the entire set of particle–hole states supported by a basis set of up to 102 orbitals. Calculations at this level give better excitation energies and intensities than SDCI methods, at substantial savings in computational effort. Comparative timings, computed results and analysis in terms of localized orbitals are given for planar and distorted ethylene using extended atomic orbital bases including diffuse functions. The results for planar ethylene are in excellent agreement with experiment.     

Computational study of iron(II) systems containing ligands with nitrogen heterocyclic groups

Density functional theory (DFT) calculations show the higher energy HOMO (highest occupied molecular orbital) orbitals of four iron(II) diimine complexes are metal centered and the lower energy LUMO (lowest unoccupied molecular orbitals) are ligand centered. The energy of the orbitals correlates with electrochemical redox potentials of the complexes. Time-dependent density functional theory (TDDFT) calculations reveal ligand centered (LC) and metal-to-ligand charge transfer (MLCT) at higher energy than experimentally observed. TDDFT calculations also reveal the presence of d-d transitions which are buried under the MLCT and LC transitions. The difference in chemical and photophysical behavior of the iron complexes compared to that of their ruthenium analogues is also addressed.     

A DFT/TDDFT interpretation of the ground and excited states of porphyrin and porphyrazine complexes

A comprehensive treatment is given of the electronic excitation spectra of Mg, Zn and Ni complexes of porphyrin and porphyrazine using time-dependent density functional theory (TDDFT). It is emphasized that the Kohn–Sham (KS) molecular orbital (MO) method, which is the basis for the TDDFT calculations, affords a MO interpretation of the ground state electronic structure and of the nature of the excitations. This implies that a direct connection can be made to many previous MO treatments of the title compounds. We discuss in particular, how the original explanations of the intensity distribution over the lowest excitations (the Q and B bands) in terms of a cyclic polyene model, or even a free-electron model, can be reconciled with the actual molecular and electronic structure of these compounds being much more complicated than these simple models. A fragment approach is used, building the porphyrin ring from pyrrole rings and CH or N bridges. This leads directly to a simple interpretation of the orbitals of Gouterman's four-orbital model, which are responsible for the Q and B bands. It also leads to additional occupied π-orbitals which are absent in the cyclic polyene model and which need to be invoked to understand other features of the electronic spectra such as the origin of the N, L and M bands. Considerable attention is given to the intensities of the various transitions, explaining why the transitions within the so-called four-orbital model of Gouterman have large transition dipoles, why transitions from additional occupied π-orbitals have relatively small transition dipoles.