Direct computation of general chemical energy differences: Application to ionization potentials, excitation, and bond energies
Abstract
Chemists are mainly interested in energy differences. In contrast, most quantum chemical methods yield the total energy which is a large number compared to the difference and has therefore to be computed to a higher relative precision than would be necessary for the difference alone. Hence, it is desirable to compute energy differences directly, thereby avoiding the precision problem. Whenever it is possible to find a parameter which transforms smoothly from an initial to a final state, the energy difference can be obtained by integrating the energy derivative with respect to that parameter (c.f., thermodynamic integration or adiabatic connection methods). If the dependence on the parameter is predominantly linear, accurate results can be obtained by single-point integration. In density functional theory (DFT) and Hartree-Fock, we applied the formalism to ionization potentials, excitation energies, and chemical bond breaking. Example calculations for ionization potentials and excitation energies showed that accurate results could be obtained with a linear estimate. For breaking bonds, we introduce a non-geometrical parameter which gradually turns the interaction between two fragments of a molecule on. The interaction changes the potentials used to determine the orbitals as well as constraining the orbitals to be orthogonal.
- Authors:
- ORNL
- Publication Date:
- Research Org.:
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States); Center for Computational Sciences
- Sponsoring Org.:
- USDOE Office of Science (SC)
- OSTI Identifier:
- 978169
- DOE Contract Number:
- DE-AC05-00OR22725
- Resource Type:
- Journal Article
- Resource Relation:
- Journal Name: Journal of Chemical Physics; Journal Volume: 125
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; ACCURACY; CHEMICAL BONDS; EXCITATION; FUNCTIONALS; IONIZATION POTENTIAL; THERMODYNAMICS
Citation Formats
Beste, Ariana, Harrison, Robert J, and Yanai, Takeshi. Direct computation of general chemical energy differences: Application to ionization potentials, excitation, and bond energies. United States: N. p., 2006.
Web. doi:10.1063/1.2244559.
Beste, Ariana, Harrison, Robert J, & Yanai, Takeshi. Direct computation of general chemical energy differences: Application to ionization potentials, excitation, and bond energies. United States. doi:10.1063/1.2244559.
Beste, Ariana, Harrison, Robert J, and Yanai, Takeshi. Sun .
"Direct computation of general chemical energy differences: Application to ionization potentials, excitation, and bond energies". United States.
doi:10.1063/1.2244559.
@article{osti_978169,
title = {Direct computation of general chemical energy differences: Application to ionization potentials, excitation, and bond energies},
author = {Beste, Ariana and Harrison, Robert J and Yanai, Takeshi},
abstractNote = {Chemists are mainly interested in energy differences. In contrast, most quantum chemical methods yield the total energy which is a large number compared to the difference and has therefore to be computed to a higher relative precision than would be necessary for the difference alone. Hence, it is desirable to compute energy differences directly, thereby avoiding the precision problem. Whenever it is possible to find a parameter which transforms smoothly from an initial to a final state, the energy difference can be obtained by integrating the energy derivative with respect to that parameter (c.f., thermodynamic integration or adiabatic connection methods). If the dependence on the parameter is predominantly linear, accurate results can be obtained by single-point integration. In density functional theory (DFT) and Hartree-Fock, we applied the formalism to ionization potentials, excitation energies, and chemical bond breaking. Example calculations for ionization potentials and excitation energies showed that accurate results could be obtained with a linear estimate. For breaking bonds, we introduce a non-geometrical parameter which gradually turns the interaction between two fragments of a molecule on. The interaction changes the potentials used to determine the orbitals as well as constraining the orbitals to be orthogonal.},
doi = {10.1063/1.2244559},
journal = {Journal of Chemical Physics},
number = ,
volume = 125,
place = {United States},
year = {Sun Jan 01 00:00:00 EST 2006},
month = {Sun Jan 01 00:00:00 EST 2006}
}
-
A direct method (D-Delta-MBPT(2)) to calculate second-order ionization potentials (IPs), electron affinities (EAs), and excitation energies is developed. The Delta-MBPT(2) method is defined as the correlated extension of the Delta-HF method. Energy differences are obtained by integrating the energy derivative with respect to occupation numbers over the appropriate parameter range. This is made possible by writing the second-order energy as a function of the occupation numbers. Relaxation effects are fully included at the SCF level. This is in contrast to linear response theory, which makes the D-Delta-MBPT(2) applicable not only to single excited but also higher excited states. We showmore »
-
Computation of the temperature of the superconducting transition for chemical compounds with A15 structure in terms of the Debye temperature and the effective and reduced ionization potentials
A formula analogous to the McMillan formula is developd is derived without use of the electron-phonon interaction parameter or the Coulomb pseudo-Potential.(AIP) -
C sub 60 and its ions: Electronic structure, ionization potentials, and excitation energies
Ab initio restricted Hartree-Fock calculations were performed on the ground states and excited states of C{sub 60}, C{sub 60}{sup +}, C{sub 60}{sup {minus}}, C{sub 60}{sup 2{minus}}, and C{sub 60}{sup 3{minus}}. Icosahedral symmetry was used throughout the calculations. The photoelectron spectrum of C{sub 60} was computed by the {Delta}SCF method. Electron affinities and excitation energies were obtained. C{sub 60}{sup {minus}} was computed to be more stable than C{sub 60}. C{sub 60}{sup 2{minus}} was estimated to be more stable than C{sub 60} while C{sub 60}{sup 3{minus}} was estimated to be higher in energy than C{sub 60}. -
Estimated core electron binding energies and their application in the interpretation of valence ionization potentials
Empirical group shifts for estimate of the core electron binding energies of compounds of 14 elements are evaluated. It is shown, for elements of the first row of the periodic table, that simple one-parameter groups shifts are inadequate for estimating core binding energies of highly positively charged atoms with polarizable ligands. In the case of carbon 1s binding energies, a method involving two parameters - an electronegativity parameter and a polarizability parameter - is devised which significantly improves the accuracy of the predictions. The nitrogen 1s binding energies of the aminoboranes, NH/sub 2/BH/sub 2/, N(CH/sub 3/)/sub 2/BF/sub 2/, NH(CH/sub 3/)BF/submore » -
TABULATED VALUES OF BOND DISSOCIATION ENERGIES, IONIZATION POTENTIALS, AND ELECTRON AFFINITIES FOR SOME MOLECULES FOUND IN HIGH-TEMPERATURE CHEMICAL REACTIONS
Values of bond dissociation energies, ionization potentials, and electron affinities are tabulated for some monatomic, diatomic, and polyatomic molecules, which are found in many high-temperature chemical reactions including combustion reactions. 113 references are given to journals and books published since 1950. (P.C.H.)