A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and openshell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr 2 dimer, exploring zeolitecatalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.Keywords quantum chemistry, software, electronic structure theory, density functional theory, electron correlation, computational modelling, Q-Chem Disciplines Chemistry CommentsThis article is from Molecular Physics: An International Journal at the Interface Between Chemistry and Physics 113 (2015): 184, doi:10.1080/00268976.2014. RightsWorks produced by employees of the U.S. Government as part of their official duties are not copyrighted within the U.S. The content of this document is not copyrighted. Authors 185A summary of the technical advances that are incorporated in the fourth major release of the Q-CHEM quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly corre...

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

A Laplace-transformed second-order Moller-Plesset perturbation theory (MP2) method is presented, which allows to achieve linear scaling of the computational effort with molecular size for electronically local structures. Also for systems with a delocalized electronic structure, a cubic or even quadratic scaling behavior is achieved. Numerically significant contributions to the atomic orbital (AO)-MP2 energy are preselected using the so-called multipole-based integral estimates (MBIE) introduced earlier by us [J. Chem. Phys. 123, 184102 (2005)]. Since MBIE provides rigorous upper bounds, numerical accuracy is fully controlled and the exact MP2 result is attained. While the choice of thresholds for a specific accuracy is only weakly dependent upon the molecular system, our AO-MP2 scheme offers the possibility for incremental thresholding: for only little additional computational expense, the numerical accuracy can be systematically converged. We illustrate this dependence upon numerical thresholds for the calculation of intermolecular interaction energies for the S22 test set. The efficiency and accuracy of our AO-MP2 method is demonstrated for linear alkanes, stacked DNA base pairs, and carbon nanotubes: e.g., for DNA systems the crossover toward conventional MP2 schemes occurs between one and two base pairs. In this way, it is for the first time possible to compute wave function-based correlation energies for systems containing more than 1000 atoms with 10 000 basis functions as illustrated for a 16 base pair DNA system on a single-core computer, where no empirical restrictions are introduced and numerical accuracy is fully preserved.

The importance of NMR spectroscopy for modern chemistry and biochemistry cannot be overestimated. Starting with the classical NMR experiments in 1946 by Purcell and Bloch, [1] contributions by numerous scientists have propelled NMR spectroscopy to be an extremely powerful tool in the investigation of structure and dynamics of molecular systems both in solution and in the solid state (e.g., reference [2]). Despite this progress, the understanding and reliable assignment of observed experimental spectra often remains a highly difficult task. Thus, theoretical methods can be extremely useful, which is the focus of this work.The most reliable way to predict NMR spectra for a specific molecular system is to calculate the NMR chemical shieldings by using quantum-chemical methods. For this purpose an entire hierarchy of methods (and basis sets) exists, which allows the exact result to be systematically approached. The only drawback is, however, the dramatic growth of the computational effort in approaching the exact solution and in increasing the number of atoms in a molecular system. Nevertheless, the hierarchy of ab initio methods allows approximate solutions to be selected and validated, so that error bars can be estimated and the simplest, reliable approximation for studying a specific class of molecular systems can be found.In recent years there has been much progress with respect to the size of molecular systems that can nowadays be treated by using Hartree-Fock (HF) and density-functional methods, due to a reduction of the scaling of the computational effort to linear. Progress has been made mostly for the calculation of energetics (e.g., see references [3][4][5][6][7][8]), the optimization of structures by using analytic gradients (e.g., references [9,10]), and for the calculation of molecular properties (e.g., references [6,11]). However, the linear-scaling computation of[*] Prof.

We present a simple but accurate preselection method based on Schwarz integral estimates to determine the significant elements of the exact exchange matrix before its evaluation, thus providing an asymptotical linear-scaling behavior for non-metallic systems. Our screening procedure proves to be highly suitable for exchange matrix calculations on massively parallel computing architectures, such as graphical processing units, for which we present a first linear-scaling exchange matrix evaluation algorithm.