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...

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design.

The conductor-like screening model (COSMO) was used to investigate the solvent influence on electronic g-values of organic radicals. The previously studied diphenyl nitric oxide and di-tert-butyl nitric oxide radicals were taken as test cases. The calculations employed spin-unrestricted density functional theory and the BP and B3LYP density functionals. The g-tensors were calculated as mixed second derivative properties with respect to the external magnetic field and the electron magnetic moment. The first-order response of the Kohn-Sham orbitals with respect to the external magnetic field was determined through the coupled-perturbed DFT approach. The spin-orbit coupling operator was treated using an accurate multicenter spin-orbit mean-field (SOMF) approach. Provided that important hydrogen bonds are explicitly modeled by a supermolecule approach and that the basis set is sufficiently saturated, the COSMO calculations lead to accurate predictions of isotropic g-shifts with deviations of not more than 100 ppm relative to experiment. Very accurate results were obtained by employing a recently developed self-consistent modification of the COSMO method to real solvents (COSMO-RS), which we briefly introduce in this paper as direct COSMO-RS (D-COSMO-RS). This model gives isotropic g-shifts of similar high accuracy for water without using the supermolecule approach. This is an important result because it solves many of the problems associated with the supermolecule approach such as local minima and the choice of a suitable model system. Thus, the self-consistent D-COSMO-RS incorporates some specific solvation effects into continuum models, in particular it appears to successfully model the effects of hydrogen bonding. Although not yet widely validated, this opens a novel approach for the calculation of properties which so far only could be calculated by the inclusion of explicit solvent molecules in continuum solvation methods.

The COSMO-RS method, a combination of the quantum chemical dielectric continuum solvation model COSMO with a statistical thermodynamics treatment for more realistic solvation (RS) simulations, has been used for the direct prediction of pKa constants of a large variety of 64 organic and inorganic acids. A highly significant correlation of r(2) = 0.984 with a standard deviation of only 0.49 between the calculated values of the free energies of dissociation and the experimental pKa values was found, without any special adjustment of the method. Thus, we have a theoretical a priori prediction method for pKa, which has the regression constant and the slope as only adjusted parameters. Such a method can be of great value in many areas of physical chemistry, especially in pharmaceutical and agrochemical industry. To our surprise, the slope of pKa vs ΔGdiss is only 58% of the theoretically expected value of 1/RTln(10). A careful analysis with respect to different contributions as well as a comparison with the work of other authors excludes the possibility that the discrepancy is due to weaknesses of the calculation method. Hence, we must conclude that the experimental pKa scale depends differently on the free energy of dissociation than generally assumed.