Polarizable Embedded RI-CC2 Method for Two-Photon Absorption Calculations

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.

Multiscale Modeling of Two-Photon Probes for Parkinson's Diagnostics Based on Monoamine Oxidase B Biomarker

Monoamine oxidase B (MAO-B) is a potential biomarker for Parkinson's disease (PD), a neurodegenerative disease associated with the loss of motor activities in human subjects. The disease state is associated with dopamine deprival, and so the inhibitors of MAO-B can serve as therapeutic drugs for PD. Since the expression level of MAO-B directly correlates to the disease progress, the distribution and population of this enzyme can be employed to monitor disease development. One of the approaches available for estimating the population is two-photon imaging. The ligands used for two-photon imaging should have high binding affinity and binding specificity toward MAO-B along with significant two-photon absorption cross sections when they are bound to the target. In this article, we study using a multiscale modeling approach, the binding affinity and spectroscopic properties (one- and two-photon absorption) of three (Flu1, Flu2, Flu3) of the currently available probes for monitoring the MAO-B level. We report that the binding affinity of the probes can be explained using the molecular size and binding cavity volume. The experimentally determined one-photon absorption spectrum is well reproduced by the employed QM/MM approaches, and the most accurate spectral shifts, on passing from one probe to another, are obtained at the coupled-cluster (CC2) level of theory. An important conclusion from this study is also the demonstration that intrinsic molecular two-photon absorption strengths (δ^2PA) increase in the order δ^2PA (Flu1) > δ^2PA (Flu2) > δ^2PA (Flu3). This is in contrast with experimental data, which predict similar values of two-photon absorption cross sections for Flu1 and Flu3. We demontrate, based on the results of electronic-structure calculations for Flu1 that this discrepancy cannot be explained by an explicit account for neighboring residues (which could lead to charge transfer between a probe and neighboring aromatic amino acids thus boosting δ^2PA). In summary, we show that the employed multiscale approach not only can optimize two-photon absorption properties and verify binding affinity, but it can also help in detailed analyses of experimental data.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

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.

Analytic Excited State Gradients for the QM/MM Polarizable Embedded Second-Order Algebraic Diagrammatic Construction for the Polarization Propagator PE-ADC(2)

An implementation of a QM/MM embedding in a polarizable environment is presented for second-order Møller-Plesset perturbation theory, MP2, for ground state energies and molecular gradients and for the second-order Algebraic Diagrammatic Construction, ADC(2), for excitation energies and excited state molecular gradients. In this implementation of PE-MP2 and PE-ADC(2), the polarizable embedded Hartree-Fock wave function is used as uncorrelated reference state. The polarization-correlation cross terms for the ground and excited states are included in this model via an approximate coupling density. A Lagrangian formulation is used to derive the relaxed electron densities and molecular gradients. The resolution-of-the-identity approximation speeds up the calculation of four-index electron repulsion integrals in the molecular orbital basis. As a first application, the method is used to study the photophysical properties of host-guest complexes where the accuracy and weaknesses of the model are also critically examined. It is demonstrated that the ground state geometries of the full quantum mechanical calculation for the supermolecule can be well reproduced. For excited state geometries, the deviations from the supermolecular calculation are slightly larger, but still the environment effects are captured qualitatively correctly, and energy gaps between the ground and excited states are obtained with sufficient accuracy.

COSMO-RI-ADC(2) excitation energies and excited state gradients

Evaluating vertical excitation energies and excited state analytic gradients in solution at COSMO-ADC(2).

Excited states in large molecular systems through polarizable embedding

Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.

Pyrrolo[3,2-b]pyrroles-From Unprecedented Solvatofluorochromism to Two-Photon Absorption

A combined experimental and theoretical study of the two-photon absorption (2PA) properties of a series of quadrupolar molecules possessing a highly electron-rich heterocyclic core, pyrrolo[3,2-b]pyrrole, is presented. In agreement with quantum-chemical calculations, large 2PA cross-section values, σ2PA ≈10(2) -10(3)  GM (1 GM=10(50)  cm(4)  s photon(-1) ), are observed at wavelengths of 650-700 nm, which correspond to the two-photon allowed but one-photon forbidden transitions. The calculations also predict that increased planarity of this molecule through removal of two N-substituents leads to further increase in the σ2PA values. Surprisingly, the most quadrupolar pyrrolo[3,2-b]pyrrole derivative, containing two 4-nitrophenyl substituents at positions 2 and 5, demonstrates a very strong solvatofluorochromic effect, with a fluorescence quantum yield as high as 0.96 in cyclohexane, whereas the fluorescence vanishes in DMSO.