Robust, Accurate, and Efficient: Quantum Embedding Using the Huzinaga Level-Shift Projection Operator for Complex Systems

Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces

The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for approaches where hybrid DFT, second-order Møller-Plesset perturbation theory, or double hybrid DFT are embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms.

Theoretical and Numerical Comparison of Quantum- and Classical Embedding Models for Optical Spectra

Quantum-mechanical (QM) and classical embedding models approximate a supermolecular quantum-chemical calculation. This is particularly useful when the supermolecular calculation has a size that is out of reach for present QM models. Although QM and classical embedding methods share the same goal, they approach this goal from different starting points. In this study, we compare the polarizable embedding (PE) and frozen-density embedding (FDE) models. The former is a classical embedding model, whereas the latter is a density-based QM embedding model. Our comparison focuses on solvent effects on optical spectra of solutes. This is a typical scenario where super-system calculations including the solvent environment become prohibitively large. We formulate a common theoretical framework for PE and FDE models and systematically investigate how PE and FDE approximate solvent effects. Generally, differences are found to be small, except in cases where electron spill-out becomes problematic in the classical frameworks. In these cases, however, atomic pseudopotentials can reduce the electron-spill-out issue.

On the accuracy of orbital based multi-level approaches for closed-shell transition metal chemistry

In this work, we investigate the accuracy of the local molecular orbital molecular orbital (LMOMO) scheme and projection-based wave function-in-density functional theory (WF-in-DFT) embedding for the prediction of reaction energies and barriers of typical reactions involving transition metals. To analyze the dependence of the accuracy on the system partitioning, we apply a manual orbital selection for LMOMO as well as the so-called direct orbital selection (DOS) for both approaches. We benchmark these methods on 30 closed shell reactions involving 16 different transition metals. This allows us to devise guidelines for the manual selection as well as settings for the DOS that provide accurate results within an error of 2 kcal mol^-1 compared to local coupled cluster. To reach this accuracy, on average 55% of the occupied orbitals have to be correlated with coupled cluster for the current test set. Furthermore, we find that LMOMO gives more reliable relative energies for small embedded regions than WF-in-DFT embedding.

Environment Effects on X-Ray Absorption Spectra With Quantum Embedded Real-Time Time-Dependent Density Functional Theory Approaches

In this work we implement the real-time time-dependent block-orthogonalized Manby-Miller embedding (rt-BOMME) approach alongside our previously developed real-time frozen density embedding time-dependent density functional theory (rt-TDDFT-in-DFT FDE) code, and investigate these methods’ performance in reproducing X-ray absorption spectra (XAS) obtained with standard rt-TDDFT simulations, for model systems comprised of solvated fluoride and chloride ions ([X@(H2O)8−, X = F, Cl). We observe that for ground-state quantities such as core orbital energies, the BOMME approach shows significantly better agreement with supermolecular results than FDE for the strongly interacting fluoride system, while for chloride the two embedding approaches show more similar results. For the excited states, we see that while FDE (constrained not to have the environment densities relaxed in the ground state) is in good agreement with the reference calculations for the region around the K and L1 edges, and is capable of reproducing the splitting of the 1s1 (n + 1)p1 final states (n + 1 being the lowest virtual p orbital of the halides), it by and large fails to properly reproduce the 1s1 (n + 2)p1 states and misses the electronic states arising from excitation to orbitals with important contributions from the solvent. The BOMME results, on the other hand, provide a faithful qualitative representation of the spectra in all energy regions considered, though its intrinsic approximation of employing a lower-accuracy exchange-correlation functional for the environment induces non-negligible shifts in peak positions for the excitations from the halide to the environment. Our results thus confirm that QM/QM embedding approaches are viable alternatives to standard real-time simulations of X-ray absorption spectra of species in complex or confined environments.

Second-Order Active-Space Embedding Theory

Quantum embedding schemes are a promising way to extend multireference computations to large molecules with strong correlation effects localized on a small number of atoms. This work introduces a second-order active-space embedding theory [ASET(2)] which improves upon mean-field frozen embedding by treating fragment-environment interactions via an approximate canonical transformation. The canonical transformation employed in ASET(2) is formulated using the driven similarity renormalization group. The ASET(2) scheme is benchmarked on the N═N bond dissociation in pentyldiazene, the S₀ to S₁ excitation in 1-octene, and the interaction energy of the O₂-benzene complex. The ASET(2) explicit treatment of fragment-environment interactions beyond the mean-field level generally improves the accuracy of embedded computations, and it becomes necessary to achieve an accurate description of excitation energies of 1-octene and the singlet-triplet gap of the O₂-benzene complex.

The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory

The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.

Projector-Based Quantum Embedding for Molecular Systems: An Investigation of Three Partitioning Approaches

Projector-based embedding is a relatively recent addition to the collection of methods that seek to utilize chemical locality to provide improved computational efficiency. This work considers the interactions between the different proposed procedures for this method and their effects on the accuracy of the results. The interplay between the embedded background, projector type, partitioning scheme, and level of atomic orbital (AO) truncation are investigated on a selection of reactions from the literature. The Huzinaga projection approach proves to be more reliable than the level-shift projection when paired with other procedural options. Active subsystem partitioning from the subsystem projected AO decomposition (SPADE) procedure proves slightly better than the combination of Pipek-Mezey localization and Mulliken population screening (PMM). Along with these two options, a new partitioning criteria is proposed based on subsystem von Neumann entropy and the related subsystem orbital occupancy. This new method overlaps with the previous PMM method, but the screening process is computationally simpler. Finally, AO truncation proves to be a robust option for the tested systems when paired with the Huzinaga projection, with satisfactory results being acquired at even the most severe truncation level.

From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape

Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.

Density functional theory based embedding approaches for transition-metal complexes

Transition metal species are commonly discussed by considering the metal atom embedded in a ligand environment. This apparently makes them interesting targets for modern embedding strategies based on Kohn-Sham density functional theory (DFT), which aim at modelling accurate predictions for large systems by combining different quantum chemical methods. In this perspective, we will focus on subsystem density functional theory and projection-based embedding. We review the developments in the field for transition metal species, demonstrate benefits, drawbacks and analyse error sources of the different strategies using the example of chromium hexacarbonyle, before giving a perspective where the field is currently heading.