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 review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

Self-assembled conjugated molecular wires containing thiophene up to 6 nm in length were grown layer-by-layer using click chemistry. Reflection-absorption infrared spectroscopy, ellipsometry and X-ray photoelectron spectroscopy were used to follow the stepwise growth. The electronic structure of the conjugated wires was studied with cyclic voltammetry and UV-vis spectroscopy as well as computationally with density functional theory (DFT). The current-voltage curves (±1 V) of the conjugated molecular wires were measured with conducting probe atomic force microscopy (CP-AFM) in which the molecular wire film bound to a gold substrate was contacted with a conductive AFM probe. By systematically measuring the low bias junction resistance as a function of length for molecules 1-4 nm long, we extracted the structure dependent tunneling attenuation factor (β) of 3.4 nm(-1) and a contact resistance of 220 kΩ. The crossover from tunneling to hopping transport was observed at a molecular length of 4-5 nm with an activation energy of 0.35 eV extracted from Arrhenius plots of resistance versus temperature. DFT calculations revealed localizations of spin densities (polarons) on molecular wire radical cations. The calculations were employed to gauge transition state energies for hopping of polarons along wire segments. Individual estimated transition state energies were 0.2-0.4 eV, in good agreement with the experimental activation energy. The transition states correspond to flattening of dihedral angles about specific imine bonds. These results open up possibilities to further explore the influence of molecular architecture on hopping transport in molecular junctions, and highlight the utility of DFT to understand charge localization and associated hopping-based transport.

A correct description of electronically excited states is critical to the interpretation of visible-ultraviolet spectra, photochemical reactions, and excited-state charge-transfer processes in chemical systems. We have recently proposed a theory called multiconfiguration pair-density functional theory (MC-PDFT), which is based on a combination of multiconfiguration wave function theory and a new kind of density functional called an on-top density functional. Here, we show that MC-PDFT with a first-generation on-top density functional performs as well as CASPT2 for an organic chemistry database including valence, Rydberg, and charge-transfer excitations. The results are very encouraging for practical applications.

The ability of Kohn−Sham density functional theory (KS-DFT) to accurately predict various types of electronic excitation energies with (necessarily approximate) exchangecorrelation functionals faces several challenges. Chief among these is that valence excitations are usually inherently multiconfigurational and therefore best treated by functionals with local exchange, whereas Rydberg and charge-transfer excitations are often better treated with nonlocal exchange. The question arises regarding whether one can optimize a functional such that all three kinds of excitations (valence, Rydberg, and charge transfer, including long-range charge transfer) are treated in a balanced and accurate way. The goal of the present work is to try to answer that question and then to optimize a functional with the best possible balanced behavior. Of the variety of functional types available, we choose to use a range-separated hybrid meta functional for the following reasons: (i) Range separation allows the percentage of Hartree−Fock (HF) exchange to change with interelectronic separation, and therefore, one can have 100% HF exchange at large interelectronic separations, which gives good performance for long-range charge-transfer excitations, while the range separation allows one to simultaneously have smaller values of HF exchange at small and intermediate interelectronic separations, giving good performance for valence and Rydberg excitations. (ii) Meta functionals allow one to obtain better accuracy with high HF exchange than is possible with functionals whose local part depends only on spin densities and their gradients. This work starts with the range-separated hybrid meta functional M11 and reoptimizes it (with strong smoothness restraints) against electronic excitation energies and ground-state properties to obtain a new functional called revM11 that gives good performance for all three types of electronic excitations and at the same time gives very good predictions across the board for ground-state properties.