Search

An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units. We demonstrate that the newmethod correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena. [ABSTRACT FROM AUTHOR]

Computational investigation of the longest wavelength excitations in a series of cyanines and linear n-acenes is undertaken with the use of standard spin-conserving linear response time-dependent density functional theory (TD-DFT) as well as its spin-flip variant and a ΔSCF method based on the ensemble DFT. The spin-conserving linear response TD-DFT fails to accurately reproduce the lowest excitation energy in these π-conjugated systems by strongly overestimating the excitation energies of cyanines and underestimating the excitation energies of n-acenes. The spin-flip TD-DFT is capable of correcting the underestimation of excitation energies of n-acenes by bringing in the non-dynamic electron correlation into the ground state; however, it does not fully correct for the overestimation of the excitation energies of cyanines, for which the non-dynamic correlation does not seem to play a role. The ensemble DFT method employed in this work is capable of correcting for the effect of missing non-dynamic correlation in the ground state of n-acenes and for the deficient description of differential correlation effects between the ground and excited states of cyanines and yields the excitation energies of both types of extended π-conjugated systems with the accuracy matching high-level ab initio multireference calculations. [ABSTRACT FROM AUTHOR]

The analytical energy gradient of the normalized elimination of the small component (NESC) method is derived for the first time and implemented for the routine calculation of NESC geometries and other first order molecular properties. Essential for the derivation is the correct calculation of the transformation matrix U relating the small component to the pseudolarge component of the wavefunction. The exact form of ∂U/∂λ is derived and its contribution to the analytical energy gradient is investigated. The influence of a finite nucleus model and that of the picture change is determined. Different ways of speeding up the calculation of the NESC gradient are tested. It is shown that first order properties can routinely be calculated in combination with Hartree-Fock, density functional theory (DFT), coupled cluster theory, or any electron correlation corrected quantum chemical method, provided the NESC Hamiltonian is determined in an efficient, but nevertheless accurate way. The general applicability of the analytical NESC gradient is demonstrated by benchmark calculations for NESC/CCSD (coupled cluster with all single and double excitation) and NESC/DFT involving up to 800 basis functions. [ABSTRACT FROM AUTHOR]

It is demonstrated that the LYP correlation functional is not suited to be used for the calculation of electron spin resonance hyperfine structure (HFS) constants, nuclear magnetic resonance spin-spin coupling constants, magnetic, shieldings and other properties that require a balanced account of opposite- and equal-spin correlation, especially in the core region. In the case of the HFS constants of alkali atoms, LYP exaggerates opposite-spin correlation effects thus invoking too strong in-out correlation effects, an exaggerated spin-polarization pattern in the core shells of the atoms, and, consequently, too large HFS constants. Any correlation functional that provides a balanced account of opposite- and equal-spin correlation leads to improved HFS constants, which is proven by comparing results obtained with the LYP and the PW91 correlation functional. It is suggested that specific response properties are calculated with the PW91 rather than the LYP correlation functional. [ABSTRACT FROM AUTHOR]

The regular approximation to the normalized elimination of the small component (NESC) in the modified Dirac equation has been developed and presented in matrix form. The matrix form of the infinite-order regular approximation (IORA) expressions, obtained in [Filatov and Cremer, J. Chem. Phys. 118, 6741 (2003)] using the resolution of the identity, is the exact matrix representation and corresponds to the zeroth-order regular approximation to NESC (NESC-ZORA). Because IORA (=NESC-ZORA) is a variationally stable method, it was used as a suitable starting point for the development of the second-order regular approximation to NESC (NESC-SORA). As shown for hydrogenlike ions, NESC-SORA energies are closer to the exact Dirac energies than the energies from the fifth-order Douglas–Kroll approximation, which is much more computationally demanding than NESC-SORA. For the application of IORA (=NESC-ZORA) and NESC-SORA to many-electron systems, the number of the two-electron integrals that need to be evaluated (identical to the number of the two-electron integrals of a full Dirac–Hartree–Fock calculation) was drastically reduced by using the resolution of the identity technique. An approximation was derived, which requires only the two-electron integrals of a nonrelativistic calculation. The accuracy of this approach was demonstrated for heliumlike ions. The total energy based on the approximate integrals deviates from the energy calculated with the exact integrals by less than 5×10-9 hartree units. NESC-ZORA and NESC-SORA can easily be implemented in any nonrelativistic quantum chemical program. Their application is comparable in cost with that of nonrelativistic methods. The methods can be run with density functional theory and any wave function method. NESC-SORA has the advantage that it does not imply a picture change. [ABSTRACT FROM AUTHOR]

A new method for calculating the indirect nuclear spin–spin coupling constant within the regular approximation to the exact relativistic Hamiltonian is presented. The method is completely analytic in the sense that it does not employ numeric integration for the evaluation of relativistic corrections to the molecular Hamiltonian. It can be applied at the level of conventional wave function theory or density functional theory. In the latter case, both pure and hybrid density functionals can be used for the calculation of the quasirelativistic spin–spin coupling constants. The new method is used in connection with the infinite-order regular approximation with modified metric (IORAmm) to calculate the spin–spin coupling constants for molecules containing heavy elements. The importance of including exact exchange into the density functional calculations is demonstrated. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A new method for relativistically corrected nuclear magnetic resonance (NMR) chemical shifts is developed by combining the individual gauge for the localized orbital approach for density functional theory with the normalized elimination of a small component using an effective potential. The new method is used for the calculation of the NMR chemical shifts of [SUP95]Mo and [SUP183]W in various molybdenum and tungsten compounds. It is shown that quasirelativistic corrections lead to an average improvement of calculated NMR chemical shift values by 300 and 120 ppm in the case of [SUP95]Mo and [SUP183]W, respectively, which is mainly due to improvements in the paramagnetic contributions. The relationship between electronic structure of a molecule and the relativistic paramagnetic corrections is discussed. Relativistic effects for the diamagnetic part of the magnetic shielding caused by a relativistic contraction of the s,p orbitals in the core region concern only the shielding values, however, have little consequence for the shift values because of the large independence from electronic structure and a cancellation of these effects in the shift values. It is shown that the relativistic corrections can be improved by level shift operators and a B3LYP hybrid functional, for which Hartree-Fock exchange is reduced to 15%. [ABSTRACT FROM AUTHOR]