Subspace Methods for the Simulation of Molecular Response Properties on a Quantum Computer.

We explore Davidson methods for obtaining excitation energies and other linear response properties within the recently developed quantum self-consistent linear response (q-sc-LR) method. Davidson-type methods allow for obtaining only a few selected excitation energies without explicitly constructing the electronic Hessian since they only require the ability to perform Hessian-vector multiplications. We apply the Davidson method to calculate the excitation energies of hydrogen chains (up to H10) and analyze aspects of statistical noise for computing excitation energies on quantum simulators. Additionally, we apply Davidson methods for computing linear response properties such as static polarizabilities for H2, LiH, H2O, OH-, and NH3, and show that unitary coupled cluster outperforms classical projected coupled cluster for molecular systems with strong correlation. Finally, we formulate the Davidson method for damped (complex) linear response, with application to the nitrogen K-edge X-ray absorption of ammonia, and the C6 coefficients of H2, LiH, H2O, OH-, and NH3.

Quantum Equation of Motion with Orbital Optimization for Computing Molecular Properties in Near-Term Quantum Computing.

Determining the properties of molecules and materials is one of the premier applications of quantum computing. A major question in the field is how to use imperfect near-term quantum computers to solve problems of practical value. Inspired by the recently developed variants of the quantum counterpart of the equation-of-motion (qEOM) approach and the orbital-optimized variational quantum eigensolver (oo-VQE), we present a quantum algorithm (oo-VQE-qEOM) for the calculation of molecular properties by computing expectation values on a quantum computer. We perform noise-free quantum simulations of BeH2 in the series of STO-3G/6-31G/6-31G* basis sets and of H4 and H2O in 6-31G using an active space of four electrons and four spatial orbitals (8 qubits) to evaluate excitation energies, electronic absorption, and, for twisted H4, circular dichroism spectra. We demonstrate that the proposed algorithm can reproduce the results of conventional classical CASSCF calculations for these molecular systems.

Quantum Mechanical Versus Polarizable Embedding Schemes: A Study of the Xray Absorption Spectra of Aqueous Ammonia and Ammonium.

The X-ray absorption spectra of aqueous ammonia and ammonium are computed using a combination of coupled cluster singles and doubles (CCSD) with different quantum mechanical and molecular mechanical embedding schemes. Specifically, we compare frozen Hartree-Fock (HF) density embedding, polarizable embedding (PE), and polarizable density embedding (PDE). Integrating CCSD with frozen HF density embedding is possible within the CC-in-HF framework, which circumvents the conventional system-size limitations of standard coupled cluster methods. We reveal similarities between PDE and frozen HF density descriptions, while PE spectra differ significantly. By including approximate triple excitations, we also investigate the effect of improving the electronic structure theory. The spectra computed using this approach show an improved intensity ratio compared to CCSD-in-HF. Charge transfer analysis of the excitations shows the local character of the pre-edge and main-edge, while the post-edge is formed by excitations delocalized over the first solvation shell and beyond.

Understanding X-ray absorption in liquid water using triple excitations in multilevel coupled cluster theory.

X-ray absorption (XA) spectroscopy is an essential experimental tool to investigate the local structure of liquid water. Interpretation of the experiment poses a significant challenge and requires a quantitative theoretical description. High-quality theoretical XA spectra require reliable molecular dynamics simulations and accurate electronic structure calculations. Here, we present the first successful application of coupled cluster theory to model the XA spectrum of liquid water. We overcome the computational limitations on system size by employing a multilevel coupled cluster framework for large molecular systems. Excellent agreement with the experimental spectrum is achieved by including triple excitations in the wave function and using molecular structures from state-of-the-art path-integral molecular dynamics. We demonstrate that an accurate description of the electronic structure within the first solvation shell is sufficient to successfully model the XA spectrum of liquid water within the multilevel framework. Furthermore, we present a rigorous charge transfer analysis of the XA spectrum, which is reliable due to the accuracy and robustness of the electronic structure methodology. This analysis aligns with previous studies regarding the character of the prominent features of the XA spectrum of liquid water.

Which Options Exist for NISQ-Friendly Linear Response Formulations?

Linear response (LR) theory is a powerful tool in classic quantum chemistry crucial to understanding photoinduced processes in chemistry and biology. However, performing simulations for large systems and in the case of strong electron correlation remains challenging. Quantum computers are poised to facilitate the simulation of such systems, and recently, a quantum linear response formulation (qLR) was introduced [Kumar et al., J. Chem. Theory Comput. 2023, 19, 9136-9150]. To apply qLR to near-term quantum computers beyond a minimal basis set, we here introduce a resource-efficient qLR theory, using a truncated active-space version of the multiconfigurational self-consistent field LR ansatz. Therein, we investigate eight different near-term qLR formalisms that utilize novel operator transformations that allow the qLR equations to be performed on near-term hardware. Simulating excited state potential energy curves and absorption spectra for various test cases, we identify two promising candidates, dubbed "proj LRSD" and "all-proj LRSD".

Genome-Driven Discovery of Hygrocins in Streptomyces rapamycinicus.

Ansamycins, represented by the antituberculosis drug rifamycin, are an important family of natural products. To obtain new ansamycins, Streptomyces rapamycinicus IMET 43975 harboring an ansamycin biosynthetic gene cluster was fermented in a 50 L scale, and subsequent purification work led to the isolation of five known and four new analogues, where hygrocin W (2) belongs to benzoquinonoid ansamycins, and the other three hygrocins, hygrocins X-Z (6-8), are new seco-hygrocins. The structures of ansamycins (1-8) were determined by the analysis of spectroscopic (1D/2D NMR and ECD) and MS spectrometric data. The Baeyer-Villiger enzyme which catalyzed the ester formation in the ansa-ring was confirmed through in vivo CRISPR base editing. The discovery of these compounds further enriches the structural diversity of ansamycins.

The variational quantum eigensolver self-consistent field method within a polarizable embedded framework.

We formulate and implement the Variational Quantum Eigensolver Self Consistent Field (VQE-SCF) algorithm in combination with polarizable embedding (PE), thereby extending PE to the regime of quantum computing. We test the resulting algorithm, PE-VQE-SCF, on quantum simulators and demonstrate that the computational stress on the quantum device is only slightly increased in terms of gate counts compared to regular VQE-SCF. On the other hand, no increase in shot noise was observed. We illustrate how PE-VQE-SCF may lead to the modeling of real chemical systems using a simulation of the reaction barrier of the Diels-Alder reaction between furan and ethene as an example.

Eremane, viscidane and isozizaene diterpenoids from the leaves of Eremophila rigida and their absolute configurations.

Previously undescribed eremane, viscidane, and isozizaene diterpenoids, eremorigidanes A-F, along with six known O-methylated flavonoids and three known triterpenoids were isolated and identified from the leaves of Eremophila rigida Chinnock by combined use of high-resolution PTP1B inhibition profiling, semipreparative- and analytical-scale HPLC separations, HPLC-PDA-HRMS analysis, and NMR spectroscopy. The absolute configuration of the unreported diterpenoids were determined by comparison of their experimental and calculated ECD spectra as well as by biosynthetic arguments. All isolates were evaluated for their PTP1B inhibitory activities, which revealed the flavonoid penduletin (3) to show inhibition with an IC50 value of 18.3 µM, and the triterpenoids 3,4-seco-olean-12-ene-3,28-dioic acid (15), oleanolic acid (16), and 3-oxo-oleanolic acid (17) to show inhibition with IC50 values of 55.7, 9.9, and 6.3 µM, respectively. The preliminary structure-activity relationship (SAR) of isolated flavonoids and triterpenoids is discussed. Plausible biosynthetic steps involved in eremane and isozizaene metabolism are presented and discussed.

Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-Accelerated Computer Architectures.

We present the development and implementation of relativistic coupled cluster linear response theory (CC-LR), which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin for the magnetic properties) as well as taking into account the finite lifetime of excited states in the framework of damped response theory. We showcase our implementation, which is capable to offload the computationally intensive tensor contractions characteristic of coupled cluster theory onto graphical processing units, in the calculation of (a) frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a particular emphasis on the calculation of valence absorption cross sections for the I2 molecule; (b) indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H2Se-H2O dimer as a prototypical system containing hydrogen bonds; and (c) optical rotations at the sodium D line for hydrogen peroxide analogues (H2Y2, Y = O, S, Se, Te). Thanks to this implementation, we are able to show the similarities in performance, but often the significant discrepancies, between CC-LR and approximate methods such as density functional theory. Comparing standard CC response theory with the flavor based upon the equation of motion formalism, we find that for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region, such as spin-spin couplings, on the other hand, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table, it may become increasingly difficult to measure pure optical rotation at the sodium D line due to the appearance of absorbing states.

Insights into localization, energy ordering, and substituent effect in excited states of azobenzenes from coupled cluster calculations of nuclear spin-induced circular dichroism.

Nuclear spin-induced circular dichroism (NSCD) is a molecular effect of differential absorption of left- and right-circularly polarized light due to nuclear spins in the molecule. In this work, new tools for its calculation are presented. Specifically, analytic expressions for the computation of the K term of NSCD have been derived and implemented for the second-order coupled cluster singles and doubles (CC2) model. NSCD results obtained thereby for three derivatives of azobenzenes have been compared with results from time-dependent density functional theory (TD-DFT). The complementary information that could be obtained from NSCD measurements compared to NMR for these three species is discussed.

Auger-Meitner and X-ray Absorption Spectra of Ethylene Cation: Insight into Conical Intersection Dynamics.

We present a theoretical investigation of the near-edge X-ray absorption fine structure and the Auger-Meitner decay spectra of ethylene and its cation. Herein, we demonstrate that our method, coupled with the nuclear ensemble approach, successfully reproduces the natural bandwidth structure of the experimental resonant Auger-Meitner decay spectra of ethylene, which is not very well reproduced within the Franck-Condon approximation. Furthermore, we analyze the Auger-Meitner decay spectra of the ethylene cation in light of minimum energy conical intersection structures involving the two lowest cationic states (D1 and D0), providing valuable insights into the ultrafast D1/D0 relaxation dynamics. Our results suggest that Auger-Meitner electron spectroscopy can help elucidate the mechanism behind the initial 20 fs of the relaxation dynamics.

Electronic Characterization of Glycolaldehyde: Experimental and Theoretical Insights from the Core- and Valence-Level Spectroscopy.

The core-level electron excitation and ionization spectra of glycolaldehyde have been investigated by photoabsorption and photoemission spectroscopy at both carbon and oxygen K-edges; the valence ionization spectra were also recorded by photoelectron spectroscopy in the UV-vis region. The spectra are interpreted by means of ab initio calculations based on the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) and coupled cluster singles, doubles, and perturbative are in good agreement with the experimental results, and many of the observed features are assigned. The photoabsorption spectra are not only dominated by transitions from core-level orbitals to unoccupied π and σ orbitals but also show structures due to Rydberg transitions.

Orthogonal Reversed-Phase C18 and Pentafluorophenyl HPLC Separation for Phytochemical Profiling of Serrulatanes in Eremophila denticulata.

Serrulatanes constitute a class of unique diterpenoids derived from all-Z nerylneryl diphosphate rather than the conventional all-E diterpenoid precursor geranylgeranyl diphosphate and thus provide an intriguing expansion of the chemical space of plant specialized metabolites. Plants of the Australian Eremophila genus are rich sources of structurally diverse serrulatanes. Here, we report the identification of 15 hitherto undescribed serrulatanes (eremoculatanes A-N), together with 16 previously reported compounds, from the EtOAc extract of Eremophila denticulata leaves. Isolation was performed by a combined use of systematic HPLC-PDA-HRMS-based phytochemical profiling and orthogonal reversed-phase C18 and pentafluorophenyl separations. Among the new compounds isolated, eremoculatane A contains a C12 backbone, for which the configuration was established by comparison of experimentally measured and theoretically calculated ECD spectra. The antihyperglycemic and antibacterial activities of the E. denticulata extract were investigated by high-resolution inhibition profiling, and they indicated that major constituents, mainly serrulatanes and flavonoids, contributed to the observed activity of the extract. One flavonoid, eupafolin (4), displayed moderate α-glucosidase inhibitory activity with an IC50 value of 41.3 µM, and four serrulatanes (8, 9, 19g, and 19j) showed more than 50% PTP1B inhibition at 200 µM.

Efficient Protocol for Computing MCD Spectra in a Broad Frequency Range Combining Resonant and Damped CC2 Quadratic Response Theory.

Coupled cluster response theory offers a path to high-accuracy calculations of spectroscopic properties, such as magnetic circular dichroism (MCD). However, divergence or slow convergence issues are often encountered for electronic transitions in high-energy regions with a high density of states. This is here addressed for MCD by an implementation of damped quadratic response theory for resolution-of-identity coupled cluster singles-and-approximate-doubles (RI-CC2), along with an implementation of the MCD A term from resonant response theory. Combined, damped and resonant response theory calculations provide an efficient strategy to obtain MCD spectra over a broad frequency range and for systems that include highly symmetric molecules with degenerate excited states. The protocol is illustrated by application to zinc tetrabenzoporphyrin in the energy region of 2-8 eV and comparison to experimental data. Timings are reported for the resonant and damped approaches, showing that a greater part of the calculation time is consumed by the construction of the building blocks for the final MCD ellipticity. A recommendation on how to use the procedure is outlined.

Description of the KLL Auger-Meitner decay spectra of argon following primary and satellite core-ionized states.

The K-edge photoelectron and KLL Auger-Meitner decay spectra of Argon have been investigated computationally at the restricted active space perturbation theory to the second order level using biorthonormally transformed orbital sets. Binding energies were computed for the Ar 1s primary ionization, as well as for satellite states originated from shake-up and shake-off processes. Based on our calculations, the contributions of shake-up and shake-off states to the KLL Auger-Meitner spectra of Argon have been completely elucidated. Our results are compared with recent state-of-the-art experimental measurements on Argon.

Identification of new PTP1B-inhibiting decipiene diterpenoid esters from Eremophila clarkei by high-resolution PTP1B inhibition profiling, enzyme kinetics analysis, and molecular docking.

In this study, an extract of the leaves of Eremophila clarkei Oldfield & F.Muell. showed protein tyrosine phosphatase 1B (PTP1B) inhibitory activity with an IC50 value of 33.0 µg/mL. The extract was therefore investigated by high-resolution PTP1B inhibition profiling to pinpoint the constituents responsible for the activity. Subsequent isolation and purification using analytical-scale HPLC led to identification of eight previously undescribed decipiene diterpenoids, eremoclarkanes A-H, as well as eremoclarkic acid, a biogenetically related new phenolic acid. In addition, one known decipiene diterpenoid and ten known O-methylated flavonoids were isolated. The structures of the isolated compounds were elucidated by extensive analysis of their HRMS and 1D and 2D NMR spectra. The absolute configuration of decipiene diterpenoids was determined by comparison of experimental and calculated ECD spectra. The flavonoid hispidulin (2b) and the four decipiene diterpenoids 13a, 13b, 13f, and 14b exhibited PTP1B inhibitory activity with IC50 values ranging from 22.8 to 33.6 µM. This is the first report of PTP1B inhibitory activity of decipienes, and enzyme kinetics revealed that 13a and 13b are competitive inhibitors of PTP1B, whereas 13f and 14b displayed mixed-type-mode inhibition of PTP1B. Finally, molecular docking indicated that 13a, 13b, 13f, and 14b showed comparable binding affinity towards the active and/or allosteric site of PTP1B enzyme. Structure-activity relationship (SAR) of the identified O-methylated flavonoids and decipiene diterpenoids towards PTP1B is discussed. Plausible enzymatic and photochemically driven routes for the formation of the decipienes and conversion products thereof are presented and discussed.

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.

A theoretical study of the time-resolved x-ray absorption spectrum of the photoionized BT-1T cation.

The time-resolved x-ray absorption spectrum of the BT-1T cation (BT-1T+) is theoretically simulated in order to investigate the charge transfer reaction of the system. We employ both trajectory surface hopping and quantum dynamics to simulate the structural evolution over time and the changes in the state populations. To compute the static x-ray absorption spectra (XAS) of the ground and excited states, we apply both the time-dependent density functional theory and the coupled cluster singles and doubles method. The results obtained are in good agreement between the methods. It is, furthermore, found that the small structural changes that occur during the reaction have little effect on the static XAS. Hence, the tr-XAS can be computed based on the state populations determined from a nuclear dynamics simulation and one set of static XAS calculations, utilizing the ground state optimized geometry. This approach can save considerable computational resources, as the static spectra need not to be calculated for all geometries. As BT-1T is a relatively rigid molecule, the outlined approach should only be considered when investigating non-radiative decay processes in the vicinity of the Franck-Condon point.