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.

and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.• Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal ? Take down policyIf you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

In the present work, we describe a more accurate and efficient variant of the chain-of-spheres algorithm (COSX) for exchange matrix computations. Higher accuracy for the numerical integration is obtained with new grids that were developed using global optimization techniques. With our new default grids, the average absolute energy errors are much lower than 0.1 kcal/mol, which is desirable to achieve “chemical accuracy.” Although the size of the new grids is increased by roughly a factor of 2.5, the excellent efficiency of the original COSX implementation is still further improved in most cases. The evaluation of the analytic electrostatic potential integrals was significantly accelerated by a new implementation of rolled-out versions of the Dupuis–Rys–King and Head-Gordon–Pople algorithms. Compared to our earlier implementation, a twofold speedup is obtained for the frequently used triple- ζ basis sets, while up to a 16-fold speedup is observed for quadruple- ζ basis sets. These large gains are a consequence of both the more efficient integral evaluation and the intermediate exchange matrix computation in a partially contracted basis when generally contracted shells occur. With our new RIJCOSX implementation, we facilitate accurate self-consistent field (SCF) binding energy calculations on a large supra-molecular complex composed of 320 atoms. The binding-energy errors with respect to the fully analytic results are well below 0.1 kcal/mol for the cc-pV(T/Q)Z basis sets and even smaller than for RIJ with fully analytic exchange. At the same time, our RIJCOSX SCF calculation even with the cc-pVQZ basis and the finest grid is 21 times faster than the fully analytic calculation.

The complete active space self-consistent-field (CASSCF) linear response method for the simulation of ultraviolet-visible (UV/Vis) absorption and electronic circular dichroism (ECD) spectra of large open-shell molecules is presented. By using a one-index transformed Hamiltonian, the computation of the most time-consuming intermediates can be pursued in an integral-direct fashion, which allows us to employ the efficient resolution-of-the-identity and overlap-fitted chain-of-spheres approximation. For the iterative diagonalization, pairs of Hermitian and anti-Hermitian trial vectors are used which facilitate, on the one hand, an efficient solution of the pair-structured generalized eigenvalue problem in the reduced space, and on the other hand, make the full multiconfigurational random phase approximation as efficient as the corresponding Tamm-Dancoff approximation. Electronic transitions are analyzed and characterized in the particle-hole picture by natural transition orbitals that are introduced for CASSCF linear response theory. For a small organic radical, we can show that the accuracy of simulated UV/Vis absorption spectra with the CASSCF linear response approach is significantly improved compared to the popular state-averaged CASSCF method. To demonstrate the efficiency of the implementation, the 50 lowest roots of a large Ni triazole complex with 231 atoms are computed for the simulated UV/Vis and ECD spectra.

We present a new implementation of a trust-region augmented Hessian approach (TRAH-SCF) for restricted and unrestricted Hartree–Fock and Kohn–Sham methods. With TRAH-SCF, convergence can always be achieved even with tight convergence thresholds, which requires just a modest number of iterations. Our convergence benchmark study and our illustrative applications focus on open-shell molecules, also antiferromagnetically coupled systems, for which it is notoriously complicated to converge the Roothaan–Hall self-consistent field (SCF) equations. We compare the number of TRAH iterations to reach convergence with those of Pulay’s original and Kollmar’s (K) variants of the direct inversion of the iterative subspace (DIIS) method and also analyze the obtained SCF solutions. Often, TRAH-SCF finds a symmetry-broken solution with a lower energy than DIIS and KDIIS. For unrestricted calculations, this is accompanied by a larger spin contamination, i.e., larger deviation from the desired spin-restricted ⟨ S2⟩ expectation value. However, there are also rare cases in which DIIS finds a solution with a lower energy than KDIIS and TRAH. In rare cases, both TRAH-SCF and KDIIS may also converge to a non- aufbau solution. For those calculations, standard DIIS always diverges. For cases that converge smoothly with either method, TRAH usually needs more iterations to converge than DIIS and KDIIS because for every new set of orbitals, the level-shifted Newton–Raphson equations are solved approximately and iteratively. In such cases, the total runtime of TRAH-SCF is still competitive with the DIIS-based approaches even if extended basis sets are employed, which is illustrated for a large hemocyanin model complex.

In this work, two approaches for simulating X-ray absorption (XA) spectra with the complete active space self-consistent field (CASSCF) linear response (LR) method are introduced. The first approach employs the well-known core-valence separation (CVS) approximation, which is predominantly used by many other electronic structure methods for simulating X-ray spectra. The second ansatz uses the harmonic Davidson algorithm for finding interior eigenvalues that lie close to a target excitation energy shift and virtually solves a shifted-and-inverted (S&I) generalized eigenvalue problem. LR-CASSCF K-edge transition energies are systematically blueshifted though have consistently smaller errors than those of the CAS or restricted active space (RAS) configuration interaction (CI) methods. For simple molecules at which the core hole can only be created at a single site, the state-specific RASSCF or n-electron valence second-order perturbation theory/RASCI gave more accurate principal K-edge excitation energies. If the core hole can be created at multiple sites, the LR-CASSCF approaches perform much better than RASSCF. Moreover, we observed that the LR-CASSCF variants were the only MR methods discussed here that predicted correctly the order of O K-edge features in the ozone molecule and the permanganate ion. The peak separation of edge features in ozone was as accurate as with equation-ofmotion coupled cluster singles and doubles. The error of the CVS approximation turned out to be very system dependent and in some cases amounted up to 1.0 eV for the K-edge excitation energies. Those CVS errors are still acceptable if one considers the observed deviation from the experimental reference by 5-11 eV. The deviations made in the XAS intensities were even more pronounced. CVS and the full S&I oscillator strengths could differ even by a factor of 2.8. Since the S&I approach is at least as efficient as the LR-CASSCF method for valence excitations, future endeavors to improve the accuracy by accounting for dynamic correlation could be built on top of the full S&I approach.

A perturbation theory-based algorithm for the iterative orbital update in complete active space self-consistent-field (CASSCF) calculations is presented. Following Angeli et al. (J. Chem. Phys. 2002, 117, 10525), the first-order contribution of singly excited configurations to the CASSCF wave function is evaluated using the Dyall Hamiltonian for the determination of a zeroth-order Hamiltonian. These authors employ an iterative diagonalization of the first-order density matrix including the first-order correction arising from single excitations, whereas the present approach uses the single-excitation amplitudes directly for the construction of the exponential of an anti-Hermitian matrix resulting in a unitary matrix which can be used for the orbital update. At convergence, the single-excitation amplitudes vanish as a consequence of the generalized Brillouin's theorem. It is shown that this approach in combination with direct inversion of the iterative subspace (DIIS) leads to very rapid convergence of the CASSCF iteration procedure.

The accuracy of three different complete active space (CAS) self-consistent field (CASSCF) methods is investigated for the electronically excited-state benchmark set of SchreiberM.SchreiberM.18397056J. Chem. Phys.2008128134110. Comparison of the CASSCF linear response (LR) methods MC-RPA and MC-TDA and the state-averaged (SA) CASSCF method is made for 122 singlet excitation energies and 69 oscillator strengths. Of all CASSCF methods, when considering the complete test set, MC-RPA performs best for both excitation energies and oscillator strengths with a mean absolute error (MAE) of 0.74 eV and 51%, respectively. MC-TDA and SA-CASSCF show a similar accuracy for the excitation energies with a MAE of ∼1 eV with respect to more accurate coupled cluster (CC3) excitation energies. The opposite trend is observed for the subset of n → π* excitation energies for which SA-CASSCF exhibits the least deviations (MAE 0.65 eV). By looking at s-tetrazine in more detail, we conclude that better performance for the n → π* SA-CASSCF excitation energies can be attributed to a fortunate error compensation. For oscillator strengths, SA-CASSCF performs worst for the complete test set (MAE 100%) as well as for the subsets of n → π* (MAE 192%) and π → π* excitations (MAE 84.9%). In general, CASSCF gives the worst performance for excitation energies of all excited-state ab initio methods considered so far due to lacking the major part of dynamic electron correlation, though MC-RPA and TD-DFT (BP86) show similar performance. Among all LR-type methods, LR-CASSCF oscillator strengths are the ones with the least accuracy for the same reason. As state-specific orbital relaxation effects are accounted for in LR-CASSCF, oscillator strengths are significantly more accurate than those of MS-CASPT2. Our findings should encourage further developments of response theory-based multireference methods with higher accuracy and feasibility.