3,132 research outputs found

    Density functional scheme for calculating the pair density

    Full text link
    The density functional scheme for calculating the pair density is presented by means of the constrained-search technique. The resultant single-particle equation takes the form of the modified Hartree-Fock equation which contains the kinetic contribution of the exchange-correlation energy functional as the correlation potential. The practical form of the kinetic contribution is also proposed with the aid of the scaling relations of the kinetic energy functionals.Comment: 5 page

    Robust methods for predicting the transition states of chemical reactions: new approaches that focus on key coordinates

    Get PDF
    A new method for optimizing transition state and minima structures using redundant internal coordinates is presented. The new method is innovative because it allows the user to select a few key reduced coordinates, whose Hessian components will be accurately computed by finite differencing; the remaining elements of the Hessian are approximated with a quasi-Newton method. Usually the reduced coordinates are the coordinates that are involved in bond breaking/forming. In order to develop this method, several other innovations were made, including ways to (a) select the key reduced coordinates automatically, (b) guess the transition state quickly and efficiently, (c) choose dihedrals so that the “linear angle problem” is avoided, (d) robustly convert redundant internal coordinates to Cartesian coordinates. These, and other technical developments (e.g., new quasi-Newton Hessians, new trust-radius updates), were validated using a database of 7000 initial transition-state guesses for a diverse set of 140 chemical reactions

    Subsystem constraints in variational second order density matrix optimization: curing the dissociative behavior

    Get PDF
    A previous study of diatomic molecules revealed that variational second-order density matrix theory has serious problems in the dissociation limit when the N-representability is imposed at the level of the usual two-index (P, Q, G) or even three-index (T1, T2) conditions [H. van Aggelen et al., Phys. Chem. Chem. Phys. 11, 5558 (2009)]. Heteronuclear molecules tend to dissociate into fractionally charged atoms. In this paper we introduce a general class of N-representability conditions, called subsystem constraints, and show that they cure the dissociation problem at little additional computational cost. As a numerical example the singlet potential energy surface of BeB+ is studied. The extension to polyatomic molecules, where more subsystem choices can be identified, is also discussed.Comment: published version;added reference
    • …
    corecore