Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Over the past decade, about 50 potential energy surfaces (PESs) for polyatomics with 4-11 atoms and for clusters have been calculated using the permutationally invariant polynomial method. This is a general, mainly linear least-squares method for precise mathematical fitting of tens of thousands of electronic energies for reactive and nonreactive systems. A brief tutorial of the methodology is given, including several recent improvements. Recent applications to the formic acid dimer (the current record holder in size for a reactive system), the H-HO complex, and four protonated water clusters [H(HO)] are given. The last application also illustrates extension to large clusters using the many-body representation.

Many model potential energy surfaces (PESs) have been reported for water; however, none are strictly from “first-principles”. Here we report such a potential, based on a many-body representation at the CCSD­(T) level of theory up to the four-body interaction. The new PES is benchmarked for the isomers of the water hexamer for dissociation energies, harmonic frequencies, and unrestricted diffusion Monte Carlo (DMC) calculations of the zero-point energies of the Prism, Book, and Cage isomers. Dissociation energies of several isomers of the 20-mer agree well with recent benchmark energies. Exploratory DMC calculations on this cluster verify the robustness of the new PES for quantum simulations. The accuracy and speed of the new PES are demonstrated for standard condensed phase properties, i.e., the radial distribution function and the self-diffusion constant. Quantum effects are shown to be substantial for these observables and also needed to bring theory into excellent agreement with experiment.