Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

Psi4 is an ab initio electronic structure program providing methods such as HartreeFock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion

A theoretical model chemistry designed to achieve high accuracy for enthalpies of formation of atoms and small molecules is described. This approach is entirely independent of experimental data and contains no empirical scaling factors, and includes a treatment of electron correlation up to the full coupled-cluster singles, doubles, triples and quadruples approach. Energies are further augmented by anharmonic zero-point vibrational energies, a scalar relativistic correction, first-order spin-orbit coupling, and the diagonal Born-Oppenheimer correction. The accuracy of the approach is assessed by several means. Enthalpies of formation (at 0 K) calculated for a test suite of 31 atoms and molecules via direct calculation of the corresponding elemental formation reactions are within 1 kJ mol(-1) to experiment in all cases. Given the quite different bonding environments in the product and reactant sides of these reactions, the results strongly indicate that even greater accuracy may be expected in reactions that preserve (either exactly or approximately) the number and types of chemical bonds.

Theoretical investigations of charge transport in organic materials are generally based on the "energy splitting in dimer" method and routinely assume that the transport parameters (site energies and transfer integrals) determined from monomer and dimer calculations can be reliably used to describe extended systems. Here, we demonstrate that this transferability can fail even in molecular crystals with weak van der Waals intermolecular interactions, due to the substantial (but often ignored) impact of polarization effects, particularly on the site energies. We show that the neglect of electronic polarization leads to qualitatively incorrect values and trends for the transfer integrals computed with the energy splitting method, even in simple prototypes such as ethylene or pentacene dimers. The polarization effect in these systems is largely electrostatic in nature and can change dramatically upon transition from a dimer to an extended system. For example, the difference in site energy for a prototypical "face-to-edge" one-dimensional stack of pentacene molecules is calculated to be 30% greater than that in the "face-to-edge" dimer, whereas the site energy difference in the pentacene crystal is vanishingly small. Importantly, when computed directly in the framework of localized monomer orbitals, the transfer integral values for dimer and extended systems are very similar.