Physical analysis of the diatomic “chemical” energy components

 The recent “chemical energy component analysis” permits the total energy of a molecule to be presented approximately but to good accuracy as a sum of atomic and diatomic energy contributions. Here the diatomic energy components are further decomposed into terms of different physical origin: electrostatics (in point-charge approximation and the distributed charge corrections), exchange effects, diatomic overlap and atomic basis extension terms. This analysis may provide us with a deeper insight into the factors influencing both the chemical bonds and the nonbonded interatomic interactions. Received: 6 May 2002 / Accepted: 13 November 2002 / Published online: 19 March 2003 Acknowledgements. The authors are indebted to the Hungarian Scientific Research Fund for partial financial support (grant no. OTKA T29716). Correspondence to: I. Mayer e-mail: mayer@chemres.hu     

Second-order Møller–Plesset perturbation theory with terms linear in the interelectronic coordinates and exact evaluation of three-electron integrals

 The second-order correlation energy of M?ller–Plesset perturbation theory is computed for the neon atom using a wave function that depends explicitly on the interelectronic coordinates (MP2-R12). The resolution-of-identity (RI) approximation, which is invoked in the standard formulation of MP2-R12 theory, is largely avoided by rigorously computing the necessary three-electron integrals. The basis-set limit for the second-order correlation energy is reached to within 0.1 mE _h. A comparison with the conventional RI-based MP2-R12 method shows that only three-electron integrals over s and p orbitals need to be computed exactly, indicating that the RI approximation can be safely used for integrals involving orbitals of higher angular momentum. Received: 9 May 2001 / Accepted: 31 October 2001 / Published online: 9 January 2002     

Unitary perturbation theory: a generalization of two-by-two rotations

A generalization of the two-by-two rotation technique is proposed, permitting a whole row (column) of the matrix to be treated simultaneously. The method is based on the explicit analytical evaluation of the matrix exponent representing a symmetric combination of the individual rotations. Besides constructing the unitary transformation matrices, a new orthogonalization algorithm is also proposed. The resulting “unitary perturbation theory” and orthogonalization method may be useful in different areas. Received: 15 November 1999 / Accepted: 30 January 2000 / Published online: 19 April 2000     

Gradient of the ZAPT2 energy

 An explicit expression for the analytical first derivative of the Z-averaged perturbation theory taken to second order energy, due to Lee and Jayatilaka, is presented for application to high-spin systems described by a restricted open-shell Hartree–Fock wavefunction. The use of frozen core orbitals is incorporated into the derivation. Received: 23 April 2001 / Accepted: 31 August 2001 / Published online: 9 January 2002     

Near-degeneracy corrections for second-order perturbation theory: comparison of two approaches

 We compare two approximate perturbation schemes which were developed recently to deal with the (quasi)degeneracy problem in many-body perturbation theory. We conclude that although the two methods were introduced on quite different theoretical grounds, their performances are quite similar, and present an improvement over traditional perturbation theory. Both methods are cheap in computation time, but cannot compete in accuracy with more sophisticated schemes such as complete-active-space perturbation theory or dressed particle theories. Received: 1 August 2000 / Accepted: 2 August 2000 / Published online: 19 January 2001     

Evaluation of the Coulomb energy in relativistic self-consistent-field theory

 Computational schemes are presented with which to evaluate the electrostatic Coulomb energy in relativistic molecular electronic structure calculations using a basis of four-component Dirac spinor amplitudes. We demonstrate that algorithms may be constructed and implemented which differ only in minor details from those in common use in nonrelativistic quantum chemistry, and that the four-component formalism is neither as complicated nor as expensive as has been suggested recently in the literature. Spherically symmetrical atomic basis sets are presented which indicate that accurate representations of the Coulomb energy may be obtained using modest expansions of the electronic density in a scalar auxiliary basis set of spherical harmonic Gaussian-type functions. Received: 15 April 2002 / Accepted: 15 May 2002 / Published online: 29 July 2002     

An intermolecular perturbation approach to hydrogen bondings in systems in the ground and excited electronic states

The intermolecular interaction energy for binary systems in the ground and excited electronic states was partitioned into the Coulomb, exchange-repulsion, induction, dispersion and charge-transfer interaction terms by the perturbation expansion method. The various interaction terms were evaluated for the hydrogen bondings in (HF)₂, (H₂O)₂, (CH₃OH)₂, (RCOOH)₂, and HF·H₂O in various geometrical configurations. It has been found that the Coulombic interaction plays a dominant role in the stability of these hydrogen bonded systems. The method was further applied to the HCOOH·H₂O codimer in both the ground and excited singlet electronic states. The results were in accord with the well-known water solvent effects on the shifts of absorption spectral bands.     

Multireference perturbation configuration interaction V. Third-order energy contributions in the Møller–Plesset and Epstein–Nesbet partitions

 A simple implementation of third-order perturbation theory applied to a multireference zero-order wavefunction is presented. Two different partitions of the Hamiltonian (M?ller–Plesset baricentric and Epstein–Nesbet) are considered. Two test cases, CH₂ and N₂, are examined. The third-order results are shown to be in good agreement in either partition and are generally an improvement with respect to the second-order results. The phenomenon of intruder states, absent in Epstein–Nesbet, appears to be magnified in the M?ller–Plesset partition. Received: 22 November 2001 / Accepted: 22 February 2002 / Published online: 3 May 2002     

Density functional theory study of the mechanism of the proline-catalyzed intermolecular aldol reaction

Transition structures associated with the C-C bond-formation step of the proline-catalyzed intermolecular aldol reaction between acetone and isobutyraldehyde have been studies using density functional theory methods at the B3LYP/6-31G** computational level. A continuum model has been selected to represent solvent effects. For this step, which is the stereocontrolling and rate-determining step, four reactive channels corresponding to the syn and anti arrangement of the active methylene of the enamine relative to the carboxylic acid group of l-proline and the re and si attack modes to both faces of the aldehyde carbonyl group have been analyzed. The B3LYP/6-31G** energies are in good agreement with experiment, allowing us to explain the origin of the catalysis and stereoselectivity for these proline-catalyzed aldol reactions. Received: 2 April 2002 / Accepted: 18 July 2002 / Published online: 11 October 2002 Acknowledgements. This work was supported by research funds provided by the Ministerio de Educación y Cultura of the Spanish Government by DGICYT (project PB98–1429). All the calculations were performed on a Cray–Silicon Graphics Origin 2000 of the Servicio de Informática de la Universidad de Valencia. We are most indebted to this center for providing us with computer capabilities. Correspondence to: L. R. Domingo e-mail: domingo@utopia.uv.es     

A multireference perturbation theory study on the vertical electronic spectrum of thiophene

The vertical electronic spectrum of the thiophene molecule is investigated by means of second and third order multireference perturbation theory (NEVPT). Single-state and quasi-degenerate NEVPT calculations of more than 25 singlet excited states have been performed. The study is addressed to the theoretical characterization of the four lowest-energy valence states, as well as the 3s, 3p and 3d Rydberg states. In addition, the excitation energies of two and valence states are also reported. For almost all the excited states, coupled cluster calculations (CCSD and CCSDR(3)) have been also carried out, using the same geometry and basis set used for the NEVPT ones, in order to make the comparison between the results of the two methods meaningful. A remarkable accordance between the NEVPT and CC excitation energies is found. The present results, over all, confirm the experimental assignments but, above all, represent an important contribution to the assignments of some low-energy states, valence and Rydberg, for which a firm interpretation is not available in the literature.     

A constant denominator perturbation theory for molecular energy

A constant denominator perturbation theory is developed based on a zeroth order Hamiltonian characterized by degenerate subsets of orbitals. Such a formulation allows for a decoupling of the numerators of the perturbation sequence, allowing for much more rapid evaluation of the resultant sums. For example, the full fourth order theory can be evaluated as an N ⁶ step rather than N ⁷, where N is proportional to the basis set.Although the theory is general, a constant denominator is chosen for this study as the difference between the average occupied and average virtual orbital energies scaled so that the first order wavefunction yields the lowest possible variational bound. The third order correction then appears naturally as a scaled Langhoff-Davidson correction. The full fourth order with this partitioning is developed. Results are presented within the localized bond model utilizing both the Pariser-Parr-Pople and CNDO/2 model Hamiltonians. The second order theory presents a useful bound, usually containing a good deal of the basis set correlation. In all cases examined the fourth order theory shows remarkable stability, even in those cases in which the Nesbet-Epstein partitioning seems poorly convergent, and the Moller-Plesset theory uncertain.     

Energy decomposition in the topological theory of atoms in molecules and in the linear combination of atomic orbitals formalism: a note

 It is shown that the molecular energy calculated at the self-consistent-field level can be strictly expressed as a sum of one- and two-atom energy components in the framework of Bader's topological theory of atoms in molecules (AIM). The expressions of our recent “chemical energy component analysis” can be obtained from the AIM ones as some linear combination of atomic orbitals mappings of the integrations over the atomic basins. Received: 15 June 2000 / Accepted: 4 October 2000 / Published online: 19 January 2001     

A singularity excluded approximate expansion scheme in relativistic density functional theory

 A singularity excluded approximate expansion (SEAX) scheme, which can be considered as one between Breit-Pauli expansion and RA expansion schemes, is proposed to expand the total energy of 4-component relativistic density functional theory. The one-electron equation can be derived variationally from the approximate total energy expression. The Hamiltonian of the one-electron equation is bounded from below and can be dealt with variationally, and the gauge dependency error in the ZORA method is essentially eliminated. It is easier to solve the SEAX equation than the IORA equation. The results related to the valence orbitals by solving the scalar SEAX equation agree very well with those by the scalar ZORA ESA method, and the results related to the inner-shell electrons of heavy elements by the two component SEAX calculations agree quite well with those by the 4-component relativistic density functional calculations. Received: 20 February 2002 / Accepted: 29 April 2002 / Published online: 8 July 2002     

Interaction energies for the water dimer by supermolecular methods and symmetry-adapted perturbation theory: the role of bond functions and convergence of basis subsets

 Using a systematic series of basis sets in supermolecular and symmetry-adapted intermolecular perturbation theory calculations it is examined how interaction energies of various water dimer structures change upon addition and shifting of bond functions. Their addition to augmented double- and triple-zeta basis sets brings the sum of the electron correlation contributions to the second-order interaction energy nearly to convergence, while accurate first-order electrostatic and exchange contributions require better than augmented quadruple-zeta quality. A scheme which combines the different perturbation energy contributions as computed in different basis subsets performs uniformly well for the various dimer structures. It yields a symmetry-adapted perturbation theory value of −21.08 kJ/mol for the energy of interaction of two vibrationally averaged water molecules compared to −21.29 kJ/mol when the full augmented triple-zeta basis set is used throughout. Received: 4 November 1999 / Accepted: 8 February 2000 / Published online: 12 May 2000     

A density functional theory study of the conformational properties of 1,2-ethanediamine: protonation and solvent effects

 The structures and the conformational energies of nonprotonated, monoprotonated and diprotonated 1,2-ethanediamine have been investigated through density functional theory. The relative performance of local and gradient-corrected functionals is discussed. The existence of hydrogen-bond formation has been determined with electron localisation function calculations. Proton affinities for nonprotonated and monoprotonated 1,2-ethanediamine have been calculated and are in agreement with experimental data. The influence of solvation has been accounted for through the self-consistent isodensity polarisable continuum model. The results for the nonprotonated conformers show that solvation stabilises those conformers which have the lone pair in an antiperiplanar conformation. Solvation of the monoprotonated conformer stabilises significantly the “anti” conformation, which is unstable in the gas phase. For the di-protonated species, solvation stabilises slightly the gauche conformer, which is unstable in the gas phase. Received: 28 September 1999 / Accepted: 2 May 2000 / Published online: 27 September 2000     

A perturbation theory without energy corrections

We employ logarithmic perturbation theory in a form that utilizes the force, instead of the potential, and yields equations for wavefunction‐correction terms without requiring any information of energy corrections at any stage. The knowledge of unperturbed eigenfunctions of other states is also unnecessary. The perturbed energy eigenvalue can be obtained as a series either by going back to the parent equation, or as an average value involving the perturbed state. The latter scheme applies to any other average property as well. Both ground and excited states of a few systems are chosen for demonstrative calculations. Influence of the nodal structure on the exponential function controlling spatial behavior of the probability density is discussed. Interrelations among specific correction terms are shown in the small‐ and large‐x regime; in the latter case, certain terms nicely sum up to yield the correct decay characteristics of the probability density as well. Relevance of the basic equation to an alternative, nonperturbative scheme of approximation is also outlined. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Convergence of multipole expanded intermolecular interaction energies for Gaussian-type-function and Slater-type-function basis sets

 It is shown that the multipole expansion of each order of the polarization series converges for large enough intermolecular distances when finite basis sets of Gaussian or Slater-type functions are used to approximate molecular response properties. Convergence of the multipole expansion for each order of the polarization series does not imply convergence of the polarization series itself. A corresponding convergence condition is extracted from the general perturbation theory in a finite-dimensional space and is applied to the H + H⁺ problem. Received: 29 September 1999 / Accepted: 22 May 2000 / Published online: 18 August 2000     

The electron affinities of transition metal atoms at the CCSD(T) and density functional levels of theory

 The electron affinities of Ti, V, Cr, Fe, Co, Ni, and Cu are computed using the density function theory and CCSD(T) approaches. Overall the CCSD(T) approach yields the best results. For this property, the B3LYP, BLYP, and BP86 functionals perform better than the BPW91, PBEPBE, and PBE1PBE ones. The accuracy of all the methods is higher if the number of 3delectrons is the same in the neutral atom and the anion. This is especially true for the density functional theory methods. Received: 23 January 2002 / Accepted: 1 April 2002 / Published online: 24 June 2002