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.

We present an approach to accurately construct the few-state model Hamiltonians for singlet fission processes on the basis of an ab initio electronic structure method tailored to dimer wave functions, called an active space decomposition strategy. In this method, the electronic structure of molecular dimers is expressed in terms of a linear combination of products of monomer states. We apply this method to tetracene and pentacene, using monomer wave functions computed by the restricted active space (RAS) method. Near-exact wave functions are computed for π-electrons of dimers that contain up to 7 × 1012 electronic configurations. Our product ansatz preserves the diabatic picture of the minimal dimer model, allowing us to accurately identify model Hamiltonians. The wave functions obtained from the model Hamiltonians account for more than 99% of the total wave functions. The resulting model Hamiltonians are shown to be converged with respect to all the parameters in the model, and corroborate previously reported coupling strengths.

We have developed an active-space decomposition strategy for molecular dimers that allows for the efficient computation of the dimer's complete-active-space wavefunction while only constructing the monomers’ active-space wavefunctions. Dimer states are formed from linear combinations of direct products of localized orthogonal monomer states and Hamiltonian matrix elements are computed directly without explicitly constructing the product space. This decomposition is potentially exact in the limit where a full set of monomer states is included. The adiabatic states are then found by diagonalizing the dimer Hamiltonian matrix. We demonstrate the convergence of our method to a complete-active-space calculation of the full dimer with two test cases: the benzene and naphthalene dimers.

We report on the efficient turbomole implementation of quadratic response properties within the time-dependent density functional theory (TDDFT) context that includes the static and dynamic dipole hyperpolarizability, ground-to-excited-state two-photon absorption amplitudes (through a single residue) and state-to-state one-photon absorption amplitudes (through a double residue). Our implementation makes full use of arbitrary (including non-Abelian) point-group symmetry as well as permutational symmetry and enables the calculation of nonlinear properties with hybrid density functionals for molecules with hundreds of atoms and thousands of basis functions at a cost that is a fixed multiple of the cost of the corresponding linear properties. Using the PBE0 hybrid density functional, we show that excited-state absorption spectra computed within the pseudowavefunction approach contain the qualitative features of transient absorption spectra tracking excimer formation in perylene diimide dimers, two-photon absorption cross sections for a series of highly twisted fused porphyrin chains are semiquantitatively reproduced, and the computed dynamic hyperpolarizability of several calix[4]arene stereoisomers yield simulated hyper-Raleigh scattering signals consistent with experiment. In addition, we show that the incorrect pole structure of adiabatic TDDFT properties can cause incorrect excited-state absorption spectra and overly resonant hyperpolarizabilities, and discuss possible remedies.

γ-Graphyne is the most symmetric sp 2 /sp 1 allotrope of carbon, which can be viewed as graphene uniformly expanded through the insertion of two-carbon acetylenic units between all the aromatic rings. To date, synthesis of bulk γ-graphyne has remained a challenge. We here report the synthesis of multilayer γ-graphyne through crystallization-assisted irreversible cross-coupling polymerization. A comprehensive characterization of this new carbon phase is described, including synchrotron powder X-ray diffraction, electron diffraction, lateral force microscopy, Raman spectroscopy, infrared spectroscopy, and cyclic voltammetry. Experiments indicate that γ-graphyne is a 0.48 eV band gap semiconductor, with a hexagonal a-axis spacing of 6.88 Å and an interlayer spacing of 3.48 Å, which is consistent with theoretical predictions. The observed crystal structure has an aperiodic sheet stacking. The material is thermally stable up to 240 °C but undergoes transformation at higher temperatures. While conventional 2D polymerization and reticular chemistry rely on error correction through reversibility, we demonstrate that a periodic covalent lattice can be synthesized under purely kinetic control. The reported methodology is scalable and inspires extension to other allotropes of the graphyne family.

Perspective: Explicitly correlated electronic structure theory for complex systems

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.

Optimized and stringent chemical methods to profile nascent RNA expression are still in demand. Herein, we expand the toolkit for metabolic labeling of RNA through application of inverse electron demand Diels–Alder (IEDDA) chemistry. Structural examination of metabolic enzymes guided the design and synthesis of vinyl-modified nucleosides, which we systematically tested for their ability to be installed through cellular machinery. Further, we tested these nucleosides against a panel of tetrazines to identify those which are able to react with a terminal alkene, but are stable enough for selective conjugation. The selected pairings then facilitated RNA functionalization with biotin and fluorophores. We found that this chemistry not only is amenable to preserving RNA integrity but also endows the ability to both tag and image RNA in cells. These key findings represent a significant advancement in methods to profile the nascent transcriptome using chemical approaches.