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 energy decomposition analysis (EDA) by Kitaura and Morokuma was redeveloped in the framework of the fragment molecular orbital method (FMO). The proposed pair interaction energy decomposition analysis (PIEDA) can treat large molecular clusters and the systems in which fragments are connected by covalent bonds, such as proteins. The interaction energy in PIEDA is divided into the same contributions as in EDA: the electrostatic, exchange‐repulsion, and charge transfer energies, to which the correlation (dispersion) term was added. The careful comparison to the ab initio EDA interaction energies for water clusters with 2–16 molecules revealed that PIEDA has the error of at most 1.2 kcal/mol (or about 1%). The analysis was applied to (H₂O)₁₀₂₄, the α helix, β turn, and β strand of polyalanine (ALA)₁₀, as well as to the synthetic protein (PDB code 1L2Y) with 20 residues. The comparative aspects of the polypeptide isomer stability are discussed in detail. © 2006 Wiley Periodicals, Inc.J Comput Chem 2007

The fragment molecular orbital (FMO) method makes possible nearly linear scaling calculations of large molecular systems, such as water clusters, proteins and DNA. In particular, FMO has been widely used in biochemical applications involving protein–ligand binding and drug design. The method has been efficiently parallelized suitable for petascale computing. Many commonly used wave functions and solvent models have been interfaced with FMO. We review the historical background of FMO, and summarize its method development and applications.

The fragment molecular-orbital (FMO) method was combined with the single-reference coupled-cluster (CC) theory. The developed method (FMO-CC) was applied at the CCSD and CCSD(T) levels of theory, for the cc-pVnZ family of basis sets (n=D,T,Q) to water clusters and glycine oligomers (up to 32 molecules/residues using as large basis sets as possible for the given system). The two- and three-body FMO-CC results are discussed at length, with emphasis on the basis-set dependence and three-body effects. Two- and three-body approximations based on interfragment distances were developed and the values appropriate for their accurate application carefully determined. The error in recovering the correlation energy was several millihartree for the two-body FMO-CC method and in the submillihartree range for the three-body FMO-CC method. In the largest calculations, we were able to perform the CCSD(T) calculations of (H2O)32 with the cc-pVQZ basis set (3680 basis functions) and (GLY)32 with the cc-VDZ basis set (712 correlated electrons). FMO-CC was parallelized using the upper level of the two-layer parallelization scheme. The computational scaling of the two-body FMO-CC method was demonstrated to be nearly linear. As an example of timings, CCSD(T) calculations of (H2O)32 with cc-pVDZ took 13min on an eight node 3.2-GHz Pentium4 cluster.