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.

We have formulated and implemented the multireference Mukherjee's coupled cluster method with connected singles, doubles, and perturbative triples [MR MkCCSD(T)] in the ACES II program package. Assessment of the new method has been performed on the first three electronic states of the oxygen molecule and on the automerization barrier of cyclobutadiene, where a comparison with other multireference CC treatments and with experimental data where available. The MR MkCCSD(T) method seems to be a promising candidate for an accurate, yet computationally tractable, treatment of systems where the static correlation plays an important role.

In this paper, we report on the development of a parallel implementation of the coupled-cluster (CC) Green function formulation (GFCC) employing single and double excitations in the cluster operator (GFCCSD). A key aspect of this work is the determination of the frequency dependent self-energy, Σ(ω). The detailed description of the underlying algorithm is provided, including approximations used that preserve the pole structure of the full GFCCSD method, thereby reducing the computational costs while maintaining an accurate character of methodology. Furthermore, for systems with strong local correlation, our formulation reveals a diagonally dominate block structure where as the non-local correlation increases, the block size increases proportionally. To demonstrate the accuracy of our approach, several examples including calculations of ionization potentials for benchmark systems are presented and compared against experiment.

In this paper we discuss a new formalism for producing an analytic coupled-cluster (CC) Green's function for an N-electron system by shifting the poles of similarity transformed Hamiltonians represented in N - 1 and N + 1 electron Hilbert spaces. Simple criteria are derived for the states in N - 1 and N + 1 electron spaces that are then corrected in the spectral resolution of the corresponding matrix representations of the similarity transformed Hamiltonian. The accurate description of excited state processes within a Green's function formalism would be of significant importance to a number of scientific communities ranging from physics and chemistry to engineering and the biological sciences. This is because the Green's function methodology provides a direct path for not only calculating properties whose underlying origins come from coupled many-body interactions but also provides a straightforward path for calculating electron transport, response, and correlation functions that allows for a direct link with experiment. As a special case of this general formulation, we discuss the application of this technique for Green's function defined by the CC with singles and doubles representation of the ground-state wave function.