Accurately solving the electronic Schrodinger equation of small atoms and molecules using explicitly correlated (r12-)MR-CI. VIII. Valence excited states of methylene (CH2)

We compute the adiabatic transition energies of methylene (CH(2)) from the ground state to the lowest electronically excited valence states using the r(12)-MR-ACPF-2 method with a large basis set and an extended reference space. We recall that this method aims at reaching the basis-set and full configuration interaction (CI) limits simultaneously. Our best excitation energies, T(e) (T(0)), are 9.22 (8.87) (a (1)A(1), corrected for relativistic and adiabatic effects), 31.98 (31.86) (b (1)B(1)), and 57.62 (57.18) kcal mol(-1) (c (1)A(1)) (both uncorrected). We are able to reach the respective basis-set limits that closely that the remaining errors of our (uncorrected) calculations are clearly due to the MR-ACPF-2 method. While we are unable to assess the error of the latter method in a systematic way, we still believe that it is rather unlikely that the errors of our excitation energies exceed +/-0.10 kcal mol(-1). We finally observe that our (corrected) a state values deviate by only -0.10 (-0.10) kcal mol(-1) from the results of Csaszar et al. [J. Chem. Phys. 118, 10631 (2003)]--who did careful extrapolations to the valence full-CI and basis-set limits and added a correction for the core correlation--and that the deviation from experiment is only -0.13 (-0.13) kcal mol(-1). From these excellent agreements we conclude that our excitation energies to the b and c states are similarly accurate.     

Registering the Amica electronic structure code in the extensible computational chemistry environment

We describe the integration and use of the Amica software package ("Atoms & Molecules In Chemical Accuracy") within the Extensible Computational Chemistry Environment (Ecce). Amica is capable of accurately solving the electronic Schrodinger equation of small atoms and molecules using terms that are linear in the interelectronic distances, r(12), on multireference level of theory, but it requires expert knowledge to configure and execute its algorithms. Ecce is a comprehensive suite of tools that support the computational chemistry research processes of computation setup, execution, and analysis through a convenient graphical user interface. Additionally, Ecce was architected with mechanisms to integrate alternative electronic structure codes. The successful integration of Amica within Ecce validates the architecture of the latter and brings the high-accuracy capabilities of Amica to a wider audience.     

Ab Initio prediction of possible crystal structures for general organic molecules

The performance of a new crystal packing procedure for the ab initio prediction of possible molecular crystal structures is presented. The method is based upon only molecular information, i.e., no unit cell parameters are assumed to be known. The search for the global crystal energy minimum and all local minima inside an energy window is derived from Monte Carlo simulated annealing methods and has been applied to various organic molecules containing heteroatoms and polar groups. A systematic evaluation of the search method and of the quality of the potential energy function has been established. It is demonstrated that the packing of general organic molecules is possible even with standard force fields like CHARMM provided that the charges defining the electrostatic interactions are based upon physical models rather than transferable empirical parameters. Concepts of crystal packing that were based till now upon assumptions and speculations could be proved or disproved by solving directly the extended global optimization problem related to crystal packing. Crystal structures of molecules as complex as those treated in this article have not been, till now, predicted by a computational approach. In one case, a disagreement between the predicted and experimental structure was evident and, based upon the computations, we suspect that the published structure is the wrong one. © 1992 by John Wiley & Sons, Inc.     

Transition-state energy and geometry, exothermicity, and van der Waals wells on the F + H2 --> FH + H ground-state surface calculated at the r12-ACPF-2 level

Explicitly correlated averaged coupled-pair functional methods have been used to compute the ground-state Born-Oppenheimer potential energy surface for the F + HH' --> FH + H' reaction at the F + HH' and FH + H' asymptotes, the F...HH', and FH...H van der Waals wells, the reaction transition state, and at points along the intrinsic reaction coordinate connecting all of these stationary points. To these energies, corrections for spin-orbit coupling and scalar relativistic effects were added to produce total electronic energies whose accuracy is demonstrated to be very high (e.g., 0.1 kcal mol(-1)). The final data are used to refine the two-body parts of the currently best three-dimensional potential energy surface for this reaction, to predict several spectroscopic parameters of the species involved, and to offer accurate estimates of the title reaction's exothermicity (32.0 kcal mol(-1)) and activation barrier (1.8 kcal mol(-1)) as well as the geometry of the transition state.