We discuss the partitioning of a quantum system through subsystem separation by unitary block-diagonalization (SSUB) applied to a Fock operator. For a one-particle Hilbert space, this separation can be formulated in a very general way. Therefore, it can be applied to very different partitionings ranging from those driven by features in the molecular structure (such as a solute surrounded by solvent molecules or an active site in an enzyme) to those that aim at an orbital separation (such as core-valence separation). Our framework embraces recent developments of Manby and Miller as well as the older ones of Huzinaga and Cantu. Projector-based embedding is simplified and accelerated by SSUB. Moreover, it directly relates to decoupling approaches for relativistic four-component many-electron theory. For a Fock operator based on the Dirac one-electron Hamiltonian, one would like to separate the so-called positronic (negative-energy) states from the electronic bound and continuum states. The exact two-component (X2C) approach developed for this purpose becomes a special case of the general SSUB framework and may therefore be viewed as a system-environment decoupling approach. Moreover, for SSUB, there exists no restriction with respect to the number of subsystems that are generated-in the limit, decoupling of all single-particle states is recovered, which represents exact diagonalization of the problem. The fact that a Fock operator depends on its eigenvectors poses challenges to all system-environment decoupling approaches and is discussed in terms of the SSUB framework. Apart from improved conceptual understanding, these relations bring about technical advances as developments in different fields can immediately cross-fertilize one another. As an important example, we discuss the atomic decomposition of the unitary block-diagonalization matrix in X2C-type approaches that can inspire approaches for the efficient partitioning of large total systems based on SSUB.
We discuss the partitioning of a quantum system through subsystem separation by unitary block-diagonalization (SSUB) applied to a Fock operator. For a one-particle Hilbert space, this separation can be formulated in a very general way. Therefore, it can be applied to very different partitionings ranging from those driven by features in the molecular structure (such as a solute surrounded by solvent molecules or an active site in an enzyme) to those that aim at an orbital separation (such as core-valence separation). Our framework embraces recent developments of Manby and Miller as well as the older ones of Huzinaga and Cantu. Projector-based embedding is simplified and accelerated by SSUB. Moreover, it directly relates to decoupling approaches for relativistic four-component many-electron theory. For a Fock operator based on the Dirac one-electron Hamiltonian, one would like to separate the so-called positronic (negative-energy) states from the electronic bound and continuum states. The exact two-component (X2C) approach developed for this purpose becomes a special case of the general SSUB framework and may therefore be viewed as a system-environment decoupling approach. Moreover, for SSUB, there exists no restriction with respect to the number of subsystems that are generated-in the limit, decoupling of all single-particle states is recovered, which represents exact diagonalization of the problem. The fact that a Fock operator depends on its eigenvectors poses challenges to all system-environment decoupling approaches and is discussed in terms of the SSUB framework. Apart from improved conceptual understanding, these relations bring about technical advances as developments in different fields can immediately cross-fertilize one another. As an important example, we discuss the atomic decomposition of the unitary block-diagonalization matrix in X2C-type approaches that can inspire approaches for the efficient partitioning of large total systems based on SSUB.
The concept of an effective Hamiltonian in Fock space is introduced. It is based on the division of the entire one-particle space into subspaces of ‘‘active’’ and ‘‘inactive’’ orbitals. The effective Fock space Hamiltonian has—for active model states—the same eigenvalues as the full Hamiltonian. The theory outlined in this context differs from that of paper I mainly in a different definition of the ‘‘diagonal part’’ of an operator, and in the fact that the ‘‘quasidegenerate case’’ applies throughout. The separation theorem, and as a consequence the connected diagram theorem, is shown to hold, in a more limited sense though, even for those normalizations where it did not in the context of universal wave and energy operators. Unlike in the theory of the ‘‘universal’’ operators of paper I the Fock space and n-particle Hilbert space approaches with analogous normalizations are no longer equivalent. In particular, the Primas normalization with a fully Lie-algebraic structure does not lead to a connected diagram expansion if it is formulated in n-particle Hilbert space, only so in a Fock space formulation. In n-particle Hilbert space with the present definition of the diagonal part of an operator the normalizations b (‘‘canonical’’) and c (‘‘Primas’’) happen to agree. As an alternative to the construction of the wave and energy operator W and L by perturbation theory the nonperturbative approach is presented as a generalization of the coupled-cluster method, in detail both in the intermediate and in the unitary normalization. In the unitary variant only a linear system for σ (the logarithm of the wave operator) has to be solved in order to get L correct through fifth order in perturbation theory with important contributions of higher orders included. A generalization of the Hartree–Fock method to Fock space theory is outlined, which guarantees stationarity of all (active) eigenstates with respect to one-particle transformations. A generalized electron pair theory is also defined. An analysis of the necessary computational steps shows that the nonperturbative approaches do not require significantly more computational effort than perturbation theory to the corresponding order. As a numerical example the H2 molecule in a small basis is discussed.
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.
Recently, some of the authors introduced the use of the Householder transformation as a simple and intuitive method for embedding local molecular fragments [see Sekaran et al., Phys. Rev. B 104, 035121 (2021) and Sekaran et al., Computation 10, 45 (2022)]. In this work, we present an extension of this approach to the more general case of multi-orbital fragments using the block version of the Householder transformation applied to the one-body reduced density matrix, unlocking the applicability to general quantum chemistry/condensed matter physics Hamiltonians. A step-by-step construction of the block Householder transformation is presented. Both physical and numerical areas of interest of the approach are highlighted. The specific mean-field (noninteracting) case is thoroughly detailed as it is shown that the embedding of a given N spin-orbital fragment leads to the generation of two separated sub-systems: (1) a 2N spin-orbitals "fragment+bath" cluster that exactly contains N electrons and (2) a remaining cluster's "environment" described by so-called core electrons. We illustrate the use of this transformation in different cases of embedding scheme for practical applications. We particularly focus on the extension of the previously introduced Local Potential Functional Embedding Theory and Householder-transformed Density Matrix Functional Embedding Theory to the case of multi-orbital fragments. These calculations are realized on different types of systems, such as model Hamiltonians (Hubbard rings) and ab initio molecular systems (hydrogen rings).
Decompositional theories describe the ways in which a global physical system can be split into subsystems, facilitating the study of how different possible partitions of a same system interplay, e.g. in terms of inclusions or signalling. In quantum theory, subsystems are usually framed as sub-C* algebras of the algebra of operators on the global system. However, most decompositional approaches have so far restricted their scope to the case of systems corresponding to factor algebras. We argue that this is a mistake: one should cater for the possibility for non-factor subsystems, arising for instance from symmetry considerations. Building on simple examples, we motivate and present a definition of partitions into an arbitrary number of parts, each of which is a possibly non-factor sub-C* algebra. We discuss its physical interpretation and study its properties, in particular with regards to the structure of algebras' centres. We prove that partitions, defined at the C*-algebraic level, can be represented in terms of a splitting of Hilbert spaces, using the framework of routed quantum circuits. For some partitions, however, such a representation necessarily retains a residual pseudo-nonlocality. We provide an example of this behaviour, given by the partition of a fermionic system into local modes.
In this paper, we discuss the extension of the recently introduced subsystem embedding subalgebra coupled cluster (SES-CC) formalism to unitary CC formalisms. In analogy to the standard single-reference SES-CC formalism, its unitary CC extension allows one to include the dynamical (outside the active space) correlation effects in an SES induced complete active space (CAS) effective Hamiltonian. In contrast to the standard single-reference SES-CC theory, the unitary CC approach results in a Hermitian form of the effective Hamiltonian. Additionally, for the double unitary CC (DUCC) formalism, the corresponding CAS eigenvalue problem provides a rigorous separation of external cluster amplitudes that describe dynamical correlation effects—used to define the effective Hamiltonian—from those corresponding to the internal (inside the active space) excitations that define the components of eigenvectors associated with the energy of the entire system. The proposed formalism can be viewed as an efficient way of downfolding many-electron Hamiltonian to the low-energy model represented by a particular choice of CAS. In principle, this technique can be extended to any type of CAS representing an arbitrary energy window of a quantum system. The Hermitian character of low-dimensional effective Hamiltonians makes them an ideal target for several types of full configuration interaction type eigensolvers. As an example, we also discuss the algebraic form of the perturbative expansions of the effective DUCC Hamiltonians corresponding to composite unitary CC theories and discuss possible algorithms for hybrid classical and quantum computing. Given growing interest in quantum computing, we provide energies for H2 and Be systems obtained with the quantum phase estimator algorithm available in the Quantum Development Kit for the approximate DUCC Hamiltonians.
A methodology is introduced for the systematic treatment of electron-electron correlations in solids and other interacting quantum N-particle systems. The method is developed within the framework of electronic structure theory (band theory) but, in contrast to conventional approaches, which are based on the single-particle picture, it is formulated within a many-particle picture in which n particles in d-dimensional phase space are treated as a single particle in a phase space of nd dimensions. In this phase space, interparticle interactions appear as external potentials allowing the treatment of the system of particles through the use of single-particle methods, while at the same time allowing a systematic, direct, and nonperturbative treatment of interparticle interactions. The method makes use of the invariance of the Hamiltonian describing an interacting-particle system under partitioning into subsystems of n particles. This treatment leads to exact results in the limit $n\ensuremath{\rightarrow}N$. Based on such partitioning, we propose a generalization of density functional theory and an appropriately defined local density approximation to treat the interactions between the n-particle units in a system of N>n particles. This approach yields n-particle correlated densities and n-particle states which can be used in an analysis of the electronic properties of materials, such as total energy, phase stability, electronic transport, and others. We use the formal construct of multiple-scattering theory to develop the method for the calculation of the two-particle electronic structure of a solid and the corresponding total energy of the ground state. We also illustrate some of the properties of the method in terms of a Hubbard model Hamiltonian on a linear ring. Various features of the method and further possible applications are presented in a discussion section.