TURBOMOLE: Today and Tomorrow.

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Pushing the limits: Efficient wavefunction methods for excited states in complex systems using frozen-density embedding.

Frozen density embedding (FDE) is an embedding method for complex environments that is simple for users to set up. It reduces the computation time by dividing the total system into small subsystems and approximating the interaction by a functional of their densities. Its combination with wavefunction methods is, however, limited to small- or medium-sized molecules because of the steep scaling in computation time of these methods. To mitigate this limitation, we present a combination of the FDE approach with pair natural orbitals (PNOs) in the TURBOMOLE software package. It combines the uncoupled FDE (FDEu) approach for excitation energy calculations with efficient implementations of second-order correlation methods in the ricc2 and pnoccsd programs. The performance of this combination is tested for tetraazaperopyrene (TAPP) molecular crystals. It is shown that the PNO truncation error on environment-induced shifts is significantly smaller than the shifts themselves and, thus, that the local approximations of PNO-based wavefunction methods can without the loss of relevant digits be combined with the FDE method. Computational wall times are presented for two TAPP systems. The scaling of the wall times is compared to conventional supermolecular calculations and demonstrates large computational savings for the combination of FDE- and PNO-based methods. Additionally, the behavior of excitation energies with the system size is investigated. It is found that the excitation energies converge quickly with the size of the embedding environment for the TAPPs investigated in the current study.

Employing Pseudopotentials to Tackle Excited-State Electron Spill-Out in Frozen Density Embedding Calculations.

In frozen density embedding (FDE), the properties of a target molecule are computed in the presence of an effective embedding potential, which accounts for the attractive and repulsive contributions of the environment. The formally exact embedding potential, however, is in practice calculated using explicit kinetic-energy functionals for which the resulting potentials are in many cases not repulsive enough to account fully for Pauli repulsion by the electrons of the environment and to compensate thereby the strong electron-nuclear attraction. For the excited states on the target molecule, this leads to charge spill-out when diffuse basis functions are included, which allow that valence electrons are excited to those regions of the environment where the strong nuclear attraction is not sufficiently compensated by repulsive contributions. To reduce this insufficiency, we propose in the present work the inclusion of atomic all-electron pseudopotentials for all environment atoms on top of the conventional embedding potential. In the current work, the pseudopotentials are applied for computing vertical excitation energies of local excited states in complex systems employing the second-order algebraic diagrammatic construction (ADC(2)) scheme. The proposed approach leads to significantly reduced charge spill-out and an improved agreement of FDE and supermolecular calculations in the frozen solvent approximation. In particular, when diffuse functions are employed, the mean absolute deviation (MAD) is reduced from 0.27 to 0.05 eV for the investigated cases.