Variational Pair-Density Functional Theory: Dealing with Strong Correlation at the Protein Scale.

Multiconfigurational pair-density functional theory (MC-PDFT) offers a promising solution to the challenges faced by traditional density functional theory (DFT) in addressing molecular systems containing transition metals, open-shells, or strong correlations in general. By utilizing both the density and on-top pair-density, MC-PDFT can make use of a more flexible multiconfigurational wave function to capture the necessary static correlation, while the pair-density functional also includes the effect of dynamic correlation. So far, MC-PDFT has been used after a multiconfigurational self-consistent field (MCSCF) step, using the orbitals and configuration interaction coefficients from the converged MCSCF wave function to compute PDFT energies and properties. Here, instead, we propose to perform a direct optimization of the wave function using the pair-density functionals, resulting in a variational formulation of MC-PDFT. We derive the expressions for the wave function gradient and illustrate their similarity to standard MCSCF equations. Furthermore, we illustrate the accuracy on a set of singlet-triplet gaps as well as dissociation curves. Our findings highlight one of MC-PDFT's standout features: a reduced dependency on the active space size compared to conventional multiconfigurational wave function methodologies. Additionally, we show that the computational cost of MC-PDFT is potentially lower than MCSCF and often on-par with standard Kohn-Sham DFT, which is demonstrated by performing a MC-PDFT calculation of the entire ferredoxin protein with 1447 atoms and nearly 12 000 basis functions.

Multiconfigurational Pair-Density Functional Theory Is More Complex than You May Think.

Multiconfigurational pair-density functional theory (MC-PDFT) is a promising way to describe both strong and dynamic correlations in an inexpensive way. The functionals in MC-PDFT are often "translated" from standard spin density functionals. However, these translated functionals can in principle lead to "translated spin densities" with a nonzero imaginary component. Current developments so far neglect this imaginary part by simply setting it to zero. In this work, we show how this imaginary component is actually needed to reproduce the correct physical behavior in a range of cases, especially low-spin open shells. We showcase the resulting formalism on both local density approximation and generalized gradient approximation functionals and illustrate the numerical behavior by benchmarking a number of singlet-triplet splittings (ST gaps) of organic diradicals and low-lying excited states of some common organic molecules. The results demonstrate that this scheme improves existing translated functionals and gives more accurate results, even with minimal active spaces.

Complex Linear Response Functions for a Multiconfigurational Self-Consistent Field Wave Function in a High Performance Computing Environment.

We present novel developments for the highly efficient evaluation of complex linear response functions of a multiconfigurational self-consistent field (MCSCF) wave function as implemented in MultiPsi. Specifically, expressions for the direct evaluation of linear response properties at given frequencies using the complex polarization propagator (CPP) approach have been implemented, within both the Tamm-Dancoff approximation (TDA) and the random phase approximation (RPA). Purely real algebra with symmetric and antisymmetric trial vectors in a shared subspace is used wherein the linear response equations are solved. Two bottlenecks of large scale MC-CPP calculations, namely, the memory footprint and computational time, are addressed. The former is addressed by limiting the size of the subspace of trial vectors by using singular value decomposition (SVD) on either orbital or CI subspaces. The latter is addressed using an efficient parallel implementation as well as the strategy of dynamically adding linear response equations at near-convergence to neighboring roots. Furthermore, a novel methodology for decomposing MC-CPP spectra in terms of intuitive orbital excitations in an approximate fashion is presented. The performance of the code is illustrated with several numerical examples, including the X-ray spectrum of a molecule with nearly one hundred atoms. Additionally, for X-ray spectroscopy, the effect of including or excluding the core orbital in the active space on small covalent metal complexes is discussed.

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Complete Active Space Methods for NISQ Devices: The Importance of Canonical Orbital Optimization for Accuracy and Noise Resilience.

To avoid the scaling of the number of qubits with the size of the basis set, one can divide the molecular space into active and inactive regions, which is also known as complete active space methods. However, selecting the active space alone is not enough to accurately describe quantum mechanical effects such as correlation. This study emphasizes the importance of optimizing the active space orbitals to describe correlation and improve the basis-dependent Hartree-Fock energies. We will explore classical and quantum computation methods for orbital optimization and compare the chemically inspired ansatz, UCCSD, with the classical full CI approach for describing the active space in both weakly and strongly correlated molecules. Finally, we will investigate the practical implementation of a quantum CASSCF, where hardware-efficient circuits must be used and noise can interfere with accuracy and convergence. Additionally, we will examine the impact of using canonical and noncanonical active orbitals on the convergence of the quantum CASSCF routine in the presence of noise.

Soft X-ray signatures of cationic manganese-oxo systems, including a high-spin manganese(V) complex.

Manganese-oxo species catalyze key reactions, including C-H bond activation or dioxygen formation in natural photosynthesis. To better understand relevant reaction intermediates, we characterize electronic states and geometric structures of [MnOn]+ manganese-oxo complexes that represent a wide range of manganese oxidation states. To this end, we apply soft X-ray spectroscopy in a cryogenic ion trap, combined with multiconfigurational wavefunction calculations. We identify [MnO2]+ as a rare high-spin manganese(V) oxo complex with key similarities to six-coordinated manganese(V) oxo systems that are proposed as reaction intermediates in catalytic dioxygen bond formation.

Unravelling the mechanism of pH-regulation in dinoflagellate luciferase.

Dinoflagellates are the dominant source of bioluminescence in coastal waters. The luminescence reaction involves the oxidation of luciferin by a luciferase enzyme, which only takes place at low pH. The pH-dependence has previously been linked to four conserved histidines. It has been suggested that their protonation might induce a conformational change in the enzyme, thereby allowing substrate access to the binding pocket. Yet, the precise mechanism of luciferase activation has remained elusive. Here, we use computational tools to predict the open structure of the luciferase in Lingulodinium polyedra and to decipher the nature of the opening mechanism. Through accelerated molecular dynamics simulations, we demonstrate that the closed-open conformational change likely takes place via a tilt of the pH-regulatory helix-loop-helix domain. Moreover, we propose that the molecular basis for the transition is electrostatic repulsion between histidine-cation pairs, which destabilizes the closed conformation at low pH. Finally, by simulating truncated mutants, we show that eliminating the C-terminus alters the shape of the active site, effectively inactivating the luciferase.

Exact semi-classical light-matter interaction operator applied to two-photon processes with strong relativistic effects.

X-ray processes involve interactions with high-energy photons. For these short wavelengths, the perturbing field cannot be treated as constant, and there is a need to go beyond the electric-dipole approximation. The exact semi-classical light-matter interaction operator offers several advantages compared to the multipole expansion such as improved stability and ease of implementation. Here, the exact operator is used to model x-ray scattering in metal K pre-edges. This is a relativistic two-photon process where absorption is dominated by electric-dipole forbidden transitions. With the restricted active space state-interaction approach, spectra can be calculated even for the multiconfigurational wavefunctions including second-order perturbation. However, as the operator itself depends on the transition energy, the cost for evaluating integrals for hundreds of thousands unique transitions becomes a bottleneck. Here, this is solved by calculating the integrals in a molecular-orbital basis that only runs over the active space, combined with a grouping scheme where the operator is the same for close-lying transitions. This speeds up the calculations of single-photon processes and is critical for the modeling of two-photon scattering processes. The new scheme is used to model Kα resonant inelastic x-ray scattering of iron-porphyrin complexes with relevance to studies of heme enzymes, for which the total computational time is reduced by several orders of magnitude with an effect on transition intensities of 0.1% or less.

Origin of core-to-core x-ray emission spectroscopy sensitivity to structural dynamics.

Recently, coherent structural dynamics in the excited state of an iron photosensitizer was observed through oscillations in the intensity of Kα x-ray emission spectroscopy (XES). Understanding the origin of the unexpected sensitivity of core-to-core transitions to structural dynamics is important for further development of femtosecond time-resolved XES methods and, we believe, generally necessary for interpretation of XES signals from highly non-equilibrium structures that are ubiquitous in photophysics and photochemistry. Here, we use multiconfigurational wavefunction calculations combined with atomic theory to analyze the emission process in detail. The sensitivity of core-to-core transitions to structural dynamics is due to a shift of the minimum energy metal-ligand bond distance between 1s and 2p core-hole states. A key effect is the additional contraction of the non-bonding 3s and 3p orbitals in 1s core-hole states, which decreases electron-electron repulsion and increases overlap in the metal-ligand bonds. The effect is believed to be general and especially pronounced for systems with strong bonds. The important role of 3s and 3p orbitals is consistent with the analysis of radial charge and spin densities and can be connected to the negative chemical shift observed for many transition metal complexes. The XES sensitivity to structural dynamics can be optimized by tuning the emission energy spectrometer, with oscillations up to ±4% of the maximum intensity for the current system. The theoretical predictions can be used to design experiments that separate electronic and nuclear degrees of freedom in ultrafast excited state dynamics.

Toward a Computational Ecotoxicity Assay.

Thousands of anthropogenic chemicals are released into the environment each year, posing potential hazards to human and environmental health. Toxic chemicals may cause a variety of adverse health effects, triggering immediate symptoms or delayed effects over longer periods of time. It is thus crucial to develop methods that can rapidly screen and predict the toxicity of chemicals to limit the potential harmful impacts of chemical pollutants. Computational methods are being increasingly used in toxicity predictions. Here, the method of molecular docking is assessed for screening potential toxicity of a variety of xenobiotic compounds, including pesticides, pharmaceuticals, pollutants, and toxins derived from the chemical industry. The method predicts the binding energy of pollutants to a set of carefully selected receptors under the assumption that toxicity in many cases is related to interference with biochemical pathways. The strength of the applied method lies in its rapid generation of interaction maps between potential toxins and the targeted enzymes, which could quickly yield molecular-level information and insight into potential perturbation pathways, aiding in the prioritization of chemicals for further tests. Two scoring functions are compared: Autodock Vina and the machine-learning scoring function RF-Score-VS. The results are promising, although hampered by the accuracy of the scoring functions. The strengths and weaknesses of the docking protocol are discussed, as well as future directions for improving the accuracy for the purpose of toxicity predictions.

Modern quantum chemistry with [Open]Molcas.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Simulations of valence excited states in coordination complexes reached through hard X-ray scattering.

Hard X-ray spectroscopy selectively probes metal sites in complex environments. Resonant inelastic X-ray scattering (RIXS) makes it is possible to directly study metal-ligand interactions through local valence excitations. Here multiconfigurational wavefunction simulations are used to model valence K pre-edge RIXS for three metal-hexacyanide complexes by coupling the electric dipole-forbidden excitations with dipole-allowed valence-to-core emission. Comparisons between experimental and simulated spectra makes it possible to evaluate the simulation accuracy and establish a best-modeling practice. The calculations give correct descriptions of all LMCT excitations in the spectra, although energies and intensities are sensitive to the description of dynamical electron correlation. The consistent treatment of all complexes shows that simulations can rationalize spectral features. The dispersion in the manganese(iii) spectrum comes from unresolved multiple resonances rather than fluorescence, and the splitting is mainly caused by differences in spatial orientation between holes and electrons. The simulations predict spectral features that cannot be resolved in current experimental data sets and the potential for observing d-d excitations is also explored. The latter can be of relevance for non-centrosymmetric systems with more intense K pre-edges. These ab initio simulations can be used to both design and interpret high-resolution X-ray scattering experiments.

Vibrational wavepacket dynamics in Fe carbene photosensitizer determined with femtosecond X-ray emission and scattering.

The non-equilibrium dynamics of electrons and nuclei govern the function of photoactive materials. Disentangling these dynamics remains a critical goal for understanding photoactive materials. Here we investigate the photoinduced dynamics of the [Fe(bmip)2]2+ photosensitizer, where bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-pyridine, with simultaneous femtosecond-resolution Fe Kα and Kß X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS). This measurement shows temporal oscillations in the XES and XSS difference signals with the same 278 fs period oscillation. These oscillations originate from an Fe-ligand stretching vibrational wavepacket on a triplet metal-centered (3MC) excited state surface. This 3MC state is populated with a 110 fs time constant by 40% of the excited molecules while the rest relax to a 3MLCT excited state. The sensitivity of the Kα XES to molecular structure results from a 0.7% average Fe-ligand bond length shift between the 1 s and 2p core-ionized states surfaces.

Core-level nonlinear spectroscopy triggered by stochastic X-ray pulses.

Stochastic processes are highly relevant in research fields as different as neuroscience, economy, ecology, chemistry, and fundamental physics. However, due to their intrinsic unpredictability, stochastic mechanisms are very challenging for any kind of investigations and practical applications. Here we report the deliberate use of stochastic X-ray pulses in two-dimensional spectroscopy to the simultaneous mapping of unoccupied and occupied electronic states of atoms in a regime where the opacity and transparency properties of matter are subject to the incident intensity and photon energy. A readily transferable matrix formalism is presented to extract the electronic states from a dataset measured with the monitored input from a stochastic excitation source. The presented formalism enables investigations of the response of the electronic structure to irradiation with intense X-ray pulses while the time structure of the incident pulses is preserved.

OpenMolcas: From Source Code to Insight.

In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

Efficient calculations of a large number of highly excited states for multiconfigurational wavefunctions.

Electronically excited states play important roles in many chemical reactions and spectroscopic techniques. In quantum chemistry, a common technique to solve excited states is the multiroot Davidson algorithm, but it is not designed for processes like X-ray spectroscopy that involves hundreds of highly excited states. We show how the use of a restricted active space wavefunction together with a projection operator to remove low-lying electronic states offers an efficient way to reach single and double-core-hole states. Additionally, several improvements to the stability and efficiency of the configuration interaction (CI) algorithm for a large number of states are suggested. When applied to a series of transition metal complexes the new CI algorithm does not only resolve divergence issues but also leads to typical reduction in computational time by 70%, with the largest savings for small molecules and large active spaces. Together, the projection operator and the improved CI algorithm now make it possible to simulate a wide range of single- and two-photon spectroscopies. © 2019 Wiley Periodicals, Inc.

Tracing the 267 nm-Induced Radical Formation in Dimethyl Disulfide Using Time-Resolved X-ray Absorption Spectroscopy.

Disulfide bonds are pivotal for the structure, function, and stability of proteins, and understanding ultraviolet (UV)-induced S-S bond cleavage is highly relevant for elucidating the fundamental mechanisms underlying protein photochemistry. Here, the near-UV photodecomposition mechanisms in gas-phase dimethyl disulfide, a prototype system with a S-S bond, are probed by ultrafast transient X-ray absorption spectroscopy. The evolving electronic structure during and after the dissociation is simultaneously monitored at the sulfur L1,2,3-edges and the carbon K-edge with 100 fs (FWHM) temporal resolution using the broadband soft X-ray spectrum from a femtosecond high-order harmonics light source. Dissociation products are identified with the help of ADC and RASPT2 electronic-structure calculations. Rapid dissociation into two CH3S radicals within 120 ± 30 fs is identified as the major relaxation pathway after excitation with 267 nm radiation. Additionally, a 30 ± 10% contribution from asymmetric CH3S2 + CH3 dissociation is indicated by the appearance of CH3 radicals, which is, however, at least partly the result of multiphoton excitation.

Viewing the Valence Electronic Structure of Ferric and Ferrous Hexacyanide in Solution from the Fe and Cyanide Perspectives.

The valence-excited states of ferric and ferrous hexacyanide ions in aqueous solution were mapped by resonant inelastic X-ray scattering (RIXS) at the Fe L2,3 and N K edges. Probing of both the central Fe and the ligand N atoms enabled identification of the metal- and ligand-centered excited states, as well as ligand-to-metal and metal-to-ligand charge-transfer excited states. Ab initio calculations utilizing the RASPT2 method were used to simulate the Fe L2,3-edge RIXS spectra and enabled quantification of the covalencies of both occupied and empty orbitals of π and σ symmetry. We found that π back-donation in the ferric complex is smaller than that in the ferrous complex. This is evidenced by the relative amounts of Fe 3d character in the nominally 2π CN(-) molecular orbital of 7% and 9% in ferric and ferrous hexacyanide, respectively. Utilizing the direct sensitivity of Fe L3-edge RIXS to the Fe 3d character in the occupied molecular orbitals, we also found that the donation interactions are dominated by σ bonding. The latter was found to be stronger in the ferric complex, with an Fe 3d contribution to the nominally 5σ CN(-) molecular orbitals of 29% compared to 20% in the ferrous complex. These results are consistent with the notion that a higher charge at the central metal atom increases donation and decreases back-donation.

Molecular Orbital Simulations of Metal 1s2p Resonant Inelastic X-ray Scattering.

For first-row transition metals, high-resolution 3d electronic structure information can be obtained using resonant inelastic X-ray scattering (RIXS). In the hard X-ray region, a K pre-edge (1s→3d) excitation can be followed by monitoring the dipole-allowed Kα (2p→1s) or Kß (3p→1s) emission, processes labeled 1s2p or 1s3p RIXS. Here the restricted active space (RAS) approach, which is a molecular orbital method, is used for the first time to study hard X-ray RIXS processes. This is achieved by including the two sets of core orbitals in different partitions of the active space. Transition intensities are calculated using both first- and second-order expansions of the wave vector, including, but not limited to, electric dipoles and quadrupoles. The accuracy of the approach is tested for 1s2p RIXS of iron hexacyanides [Fe(CN)6](n-) in ferrous and ferric oxidation states. RAS simulations accurately describe the multiplet structures and the role of 2p and 3d spin-orbit coupling on energies and selection rules. Compared to experiment, relative energies of the two [Fe(CN)6](3-) resonances deviate by 0.2 eV in both incident energy and energy transfer directions, and multiplet splittings in [Fe(CN)6](4-) are reproduced within 0.1 eV. These values are similar to what can be expected for valence excitations. The development opens the modeling of hard X-ray scattering processes for both solution catalysts and enzymatic systems.

Analytical State-Average Complete-Active-Space Self-Consistent Field Nonadiabatic Coupling Vectors: Implementation with Density-Fitted Two-Electron Integrals and Application to Conical Intersections.

Analytical state-average complete-active-space self-consistent field derivative (nonadiabatic) coupling vectors are implemented. Existing formulations are modified such that the implementation is compatible with Cholesky-based density fitting of two-electron integrals, which results in efficient calculations especially with large basis sets. Using analytical nonadiabatic coupling vectors, the optimization of conical intersections is implemented within the projected constrained optimization method. The standard description and characterization of conical intersections is reviewed and clarified, and a practical and unambiguous system for their classification and interpretation is put forward. These new tools are subsequently tested and benchmarked for 19 different conical intersections. The accuracy of the derivative coupling vectors is validated, and the information that can be drawn from the proposed characterization is discussed, demonstrating its usefulness.