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.

A convenient and inherently more secure communication channel for encoding messages via specifically designed molecular keys is introduced by combining advanced encryption standard cryptography with molecular steganography. The necessary molecular keys require large structural diversity, thus suggesting the application of multicomponent reactions. Herein, the Ugi four-component reaction of perfluorinated acids is utilized to establish an exemplary database consisting of 130 commercially available components. Considering all permutations, this combinatorial approach can unambiguously provide 500,000 molecular keys in only one synthetic procedure per key. The molecular keys are transferred nondigitally and concealed by either adsorption onto paper, coffee, tea or sugar as well as by dissolution in a perfume or in blood. Re-isolation and purification from these disguises is simplified by the perfluorinated sidechains of the molecular keys. High resolution tandem mass spectrometry can unequivocally determine the molecular structure and thus the identity of the key for a subsequent decryption of an encoded message.

Nuclear shielding calculations for local hybrid (LH) functionals with position-dependent exact-exchange admixtures within a coupled-perturbed Kohn−Sham (CPKS) framework have been implemented into the TURBOMOLE code using efficient seminumerical integration techniques to deal with two-electron integrals. When using gauge-including atomic orbitals, LHs generate additional terms within the "pre-loop" section of the CPKS scheme compared to global hybrid (GH) functionals, related to perturbed electron-repulsion integrals. These terms have been implemented and tested in detail, together with dependencies on grid sizes and integral screening procedures. Even with relatively small grids, a seminumerical treatment of GHs reproduces analytical GH results with high accuracy while improving scaling with system and basis-set sizes significantly. The extra terms generated by LHs in the pre-loop part increase the scaling of that contribution slightly, but the advantages compared to the analytical scheme are largely retained, in particular for the typically large basis sets used in NMR shift calculations, allowing for a very efficient computational scheme. An initial comparison of four first-generation LHs based on LDA exchange for a shielding test set of 15 small main-group molecules against high-level CCSD(T) benchmark data indicates a substantial reduction of the systematically underestimated shieldings compared to semilocal functionals or GHs for non-hydrogen nuclei when a so-called t-LMF is used to control the position dependence of the exact-exchange admixture. In contrast, proton shieldings are underestimated with this LMF, while an LH with a so-called s-LMF performs much better. These results are discussed in the context of experience for other properties, and they suggest directions for further improvements of LHs regarding nuclear shieldings.

We present a highly efficient implementation for density functional calculations of chemical shielding constants. It employs the multipole-accelerated resolution of the identity for the calculation of the Coulomb part, which complements the usage of low order scaling routines for the evaluation of the exchange-correlation part and the Hartree-Fock exchange part. Introduced errors for shifts of chemical shielding constants of H, C, F, and P are evaluated for respective test sets of molecules and are related to the accuracy of shifts obtained with hybrid and nonhybrid functionals of the generalized gradient approximation type as well as for meta-GGA functionals themselves. Efficiency is demonstrated for α-d-glucose chains with more than 2500 atoms on a single CPU as well as with an OpenMP parallelized version.

The parent phosphinidene SIMesPH (1) (SIMes = 1,3-bis(2,4,6-trimethylphenyl)imidazolin-2-ylidene) was treated with benzyl potassium to form the metalated species SIMesPK (2), which was used as a synthon for the new group 13/15 cycles and cage compounds [SIMesPGatBu] (3), [SIMesP(GatBu)Cl] (4) and [K(SIMesP)AltBu] (5).