A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree–Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized.

The 1,3-dipolar cycloadditions of ozone to ethyne and ethene provide extreme examples of multireference singlet-state chemistry, and they are examined here to test the applicability of several approaches to thermochemical kinetics of systems with large static correlation. Four different multireference diagnostics are applied to measure the multireference characters of the reactants, products, and transition states; all diagnostics indicate significant multireference character in the reactant portion of the potential energy surfaces. We make a more complete estimation of the effect of quadruple excitations than was previously available, and we use this with CCSDT/CBS estimation of Wheeler et al. (Wheeler, S. E.; Ess, D. H.; Houk, K. N. J. Phys. Chem. A 2008, 112, 1798 to make new best estimates of the van der Waals association energy, the barrier height, and the reaction energy to form the cycloadduct for both reactions. Comparing with these best estimates, we present comprehensive mean unsigned errors for a variety of coupled cluster, multilevel, and density functional methods. Several computational aspects of multireference reactions are considered: (i) the applicability of multilevel theory, (ii) the convergence of coupled cluster theory for reaction barrier heights, (iii) the applicability of completely renormalized coupled cluster methods to multireference systems, (iv) the treatment by density functional theory, (v) the multireference perturbation theory for multireference reactions, and (vi) the relative accuracy of scaling-type multilevel methods as compared with additive ones. It is found that scaling-type multilevel methods do not perform better than the additive-type multilevel methods. Among the 48 tested density functionals, only M05 reproduces the best estimates within their uncertainty. Multireference perturbation theory based on the complete-active-space reference wave functions constructed using a small number of reaction-specific active orbitals gives accurate forward barrier heights; however, it significantly underestimates reaction energies.

Shifts in the excitation energy of the organic chromophore, cis-7-hydroxyquinoline (cis-7HQ), corresponding to the π→π* transition in cis-7HQ and induced by the complexation with a variety of small hydrogen-bonded molecules, obtained with the frozen-density embedding theory (FDET), are compared with the results of the supermolecular equation-of-motion coupled-cluster (EOMCC) calculations with singles, doubles, and non-iterative triples, which provide the reference theoretical data, the supermolecular time-dependent density functional theory (TDDFT) calculations, and experiment. Unlike in the supermolecular EOMCC and TDDFT cases, where each complexation-induced spectral shift is evaluated by performing two separate calculations, one for the complex and another one for the isolated chromophore, the FDET shifts are evaluated as the differences of the excitation energies determined for the same many-electron system, representing the chromophore fragment with two different effective potentials. By considering eight complexes of cis-7HQ with up to three small hydrogen-bonded molecules, it is shown that the spectral shifts resulting from the FDET calculations employing non-relaxed environment densities and their EOMCC reference counterparts are in excellent agreement with one another, whereas the analogous shifts obtained with the supermolecular TDDFT method do not agree with the EOMCC reference data. The average absolute deviation between the complexation-induced shifts, which can be as large, in absolute value, as about 2000 cm-1, obtained using the non-relaxed FDET and supermolecular EOMCC approaches that represent two entirely different computational strategies, is only about 100 cm-1, i.e., on the same order as the accuracy of the EOMCC calculations. This should be contrasted with the supermolecular TDDFT calculations, which produce the excitation energy shifts that differ from those resulting from the reference EOMCC calculations by about 700 cm-1 on average. Among the discussed issues are the choice of the electronic density defining the environment with which the chromophore interacts, which is one of the key components of FDET considerations, the basis set dependence of the FDET, supermolecular TDDFT, and EOMCC results, the usefulness of the monomer vs supermolecular basis expansions in FDET considerations, and the role of approximations that are used to define the exchange-correlation potentials in FDET and supermolecular TDDFT calculations

The methylcobalamin cofactor (MeCbl), which is one of the biologically active forms of vitamin B12, has been the subject of many spectroscopic and theoretical investigations. Traditionally, the lowest-energy part of the photoabsorption spectrum of MeCbl (the so-called α/β band) has been interpreted as an S0→S1 electronic transition dominated by π→π* excitations associated with the C=C stretching of the corrin ring. However, a more quantitative band-shape analysis of the α/β spectral region, along with circular dichroism (CD), magnetic CD, and resonance Raman data, has revealed the presence of a second electronic transition that involves the Co-C(Me) bond weakening. Conversely, the lowest-energy excitations based on transient absorption spectroscopy measurements have been interpreted as metal-to-ligand charge transfer (MLCT) transitions. To resolve the existing controversy about the interpretation of the S1 state of MeCbl, calculations have been performed using two independent ab initio wavefunction-based methods. These include the modified variant of the second-order multiconfigurational quasi-degenerate perturbation theory (MC-XQDPT2), using complete active space self-consistent field orbitals, and the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) approach using restricted Hartree-Fock orbitals. It is shown that both ab initio methods provide a consistent description of the S1 state as having an MLCT character. In addition, the performance of different types of functionals, including hybrid (B3LYP, MPW1PW91, TPSSh), generalized-gradient-approximation-type (GGA-type) (BP86, BLYP, MPWPW91), meta-GGA (TPSS), and range-separated (CAM-B3LYP, LC-BLYP) approaches, has been examined and the results of the corresponding time-dependent density functional theory calculations have been benchmarked against the MC-XQDPT2 and EOM-CCSD data. The hybrid functionals support the interpretation in which the S1 state represents a π→π* transition localized on corrin, while pure GGA, meta-GGA, and LC-BLYP functionals produce results consistent with the MLCT assignment.