A critical study of basis set effects and the use of approximate natural orbitals in SCF-CI calculations of molecular geometries and heats of reaction

SCF-CI calculations have been performed on a number of chemical reactions between closed shell molecules in order to determine the heats of reaction. Contracted Gaussian type atomic basis sets of three different qualities were used and the CI calculations were performed in a truncated approximate natural orbital space. The conclusions to be drawn from these calculations are rather pessimistic. For heats of reaction, errors up to 6 kcal/mole are obtained on the SCF-level with a double zeta plus polarization atomic basis. A further improvement is only possible if extended basis sets are used. Correlation effects on heats of reaction are of the same size and CI calculations are therefore only meaningful with large atomic basis sets.For the CI calculations a one-electron space of approximate natural orbitals, obtained from second order RS perturbation theory, was used. Different truncations, using the occupation number as criterion, were tested. The general conclusion is that errors in energy differences obtained with a truncated basis set are of the same magnitude as the error in the total correlation energy. In practice this means that not more than 20–30% of the approximate natural orbitals can be deleted if the error is to be kept less than a few kcal/mole.Finally the truncation error in calculations of bond distances was tested for a few cases. Errors of around 10% of the total change due to correlation were found when 30% of the lowest occupied natural orbitals were deleted.     

Toward cocaine esterase therapeutics

Cocaine is among the most reinforcing of all drugs of abuse, yet no effective pharmacotherapy is available. Herein, we report the development and characterization of phage-displayed cocaine esterases with pharmacologically relevant kinetic parameters (kcat/Km approximately 104 M-1 s-1).     

Direct dynamics study of ultrafast vibrational energy relaxation in ice Ih

Car-Parrinello molecular dynamics (CPMD) and a previously developed wave packet model are used to study ultrafast relaxation in water clusters. Water clusters of 15 water molecules are used to represent ice Ih. The relaxation is studied by exciting a symmetric or an asymmetric stretch mode of the central water molecule. The CPMD results suggest that relaxation occurs within 100 fs. This is in agreement with experimental work by Woutersen and Bakker and the earlier wave packet calculations. The CPMD results further indicate that the excitation energy is transferred both intramolecularly and intermolecularly on roughly the same time scale. The intramolecular energy transfer occurs predominantly between the symmetric and asymmetric modes while the bend mode is largely left unexcited on the short time scale studied here.     

A CNDO/2 calculation on stilbene

The equilibrium conformations of trans- and cis-stilbene and the rotational barrier between the two isomers were calculated using the CNDO/2 method. The results for the equilibrium angles and distances are in good agreement with experimental data. However, according to the calculations the phenyl rings would have a nearly free rotation which is only limited by steric hindrance. This result obviously deviates from the current view. The results are discussed in relation to previous theoretical and experimental evidence.

---

Zusammenfassung Die Gleichgewichtskonformation von trans- und cis-Stilben und die Rotationsbarriere zwischen zwei Isomeren wurde nach der CNDO/2-Methode berechnet. Die Ergebnisse in bezug auf die Geometrie stimmen mit experimentellen Daten gut überein, die Rotationsbarriere für die Drehung der Phenylringe dagegen ist praktisch Null.

---

Résumé La méthode CNDO/2 a été employée pour calculer les conformations d'équilibre du trans- et du cis-stilbène ainsi que la barrière de rotation entre ces deux isomères. Angles d'équilibre et distances calculés sont en bon accord avec les données expérimentales. Cependant les calculs indiquent que les cycles phényles devraient posséder une rotation libre soumise seulement à l'empêchement stérique. Ce résultat s'écarte de l'opinion courante. Les résultats sont discutés en relation avec les données théoriques et expérimentales antérieures.

---

---

This work was supported by grant No 2741-4 from The Swedish Natural Science Research Council and grant No 69-927/U 697 from The Swedish Board for Technical Development.     

A multiconfigurational scf study of the transformation of methylene peroxide to dioxirane Intermediates in the ozonolysis of ethylene

The transformation of methylene peroxide to dioxirane, molecules of importance for understanding the mechanism of ozonolysis of ethylene, has been studied with ab initio MC SCF calculations. The results, an exothermicity of 106 kJ/mol and reaction barrier of 90 kJ/mol, support the Criegee mechanism for liquid-phase ozonolysis and the idea that dioxirane is important for gas-phase ozonolysis of ethylene.     

A new method for approximate estimates of the dispersion interaction between two molecules

A new method to evaluate the dispersion interaction between two weakly interacting closed shell molecules is presented. The method is based on a second-order sum-over-states perturbation method, where the nominator is approximated by products of one electron integrals.