Biomedical subjects

Oleg Gritsenko

Publications and source records attributed to Oleg Gritsenko.

3 recordsLinked to original sources

PubMed2006

A simple natural orbital mechanism of "pure" van der Waals interaction in the lowest excited triplet state of the hydrogen molecule.

A treatment of van der Waals (vdW) interaction by density-matrix functional theory requires a description of this interaction in terms of natural orbitals (NOs) and their occupation numbers. From an analysis of the configuration-interaction (CI) wave function of the 3Sigmau + state of H2 and the exact NO expansion of the two-electron triplet wave function, we demonstrate that the construction of such a functional is straightforward in this case. A quantitative description of the vdW interaction is already obtained with, in addition to the standard part arising from the Hartree-Fock determinant /1sigmag(r1)1sigmau(r2)/, only two additional terms in the two-electron density, one from the first "excited" determinant /2sigmag(r1)2sigmau(r2)/ and one from the state of 3Sigmau + symmetry belonging to the (1pig)1(1piu)1 configuration. The potential-energy curve of the 3Sigmau + state calculated around the vdW minimum with the exact density-matrix functional employing only these eight NOs and NO occupations is in excellent agreement with the full CI one and reproduces well the benchmark potential curve of Kolos and Wolniewicz [J. Chem. Phys. 43, 2429 (1965)]. The corresponding terms in the two-electron density rho2(r1,r2), containing specific products of NOs combined with prefactors that depend on the occupation numbers, can be shown to produce exchange-correlation holes that correspond precisely to the well-known intuitive picture of the dispersion interaction as an instantaneous dipole-induced dipole (higher multipole) effect. Indeed, (induced) higher multipoles account for almost 50% of the total vdW bond energy. These results serve as a basis for both a density-matrix functional theory of van der Waals bonding and for the construction of orbital-dependent functionals in density-functional theory that could be used for this type of bonding.

Journal Article↗

PubMed2005

An improved density matrix functional by physically motivated repulsive corrections.

An improved density matrix functional [correction to Buijse and Baerends functional (BBC)] is proposed, in which a hierarchy of physically motivated repulsive corrections is employed to the strongly overbinding functional of Buijse and Baerends (BB). The first correction C1 restores the repulsive exchange-correlation (xc) interaction between electrons in weakly occupied natural orbitals (NOs) as it appears in the exact electron pair density rho(2) for the limiting two-electron case. The second correction C2 reduces the xc interaction of the BB functional between electrons in strongly occupied NOs to an exchange-type interaction. The third correction C3 employs a similar reduction for the interaction of the antibonding orbital of a dissociating molecular bond. In addition, C3 applies a selective cancellation of diagonal terms in the Coulomb and xc energies (not for the frontier orbitals). With these corrections, BBC still retains a correct description of strong nondynamical correlation for the dissociating electron pair bond. BBC greatly improves the quality of the BB potential energy curves for the prototype few-electron molecules and in several cases BBC reproduces very well the benchmark ab initio potential curves. The average error of the self-consistent correlation energies obtained with BBC3 for prototype atomic systems and molecular systems at the equilibrium geometry is only ca. 6%.

Journal Article↗

PubMed2004

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)<<omegaCT CISD(R) at large R. The TDDFT curve obtained with fxc (asymp) however approaches omegaCT CISD(R) closely in the region R=3-10 A. Then, the adequate rigorous TDDFT approach should interpolate between the LDA/GGA ALDA xc kernel for excitations in compact systems and fxc(asymp) for weakly interacting fragments and suitable interpolation expressions are considered.

Journal Article↗