An efficient approach to calculate approximate pure-state and transition reduced density matrices in the framework of the multireference relativistic Fock-space coupled cluster (FS CC) theory is proposed. The method is based on the effective operator formalism and consists of the direct substitution of the FS CC Ansatz for a wave operator into the effective operator expression with the subsequent truncation of expansion at the terms quadratic in cluster amplitudes. The final density matrix is defined by active-space density matrices of different ranks ‘dressed’ with contributions from cluster operators. The method gives a connected expression for pure-state density matrices, provided that the intermediate normalisation condition is fulfilled. Moreover, under some additional assumptions, the connectivity can also be ensured for calculated transition property matrix elements and natural transition spinors. The developed technique allows for fast and accurate calculations of one-particle reduced density matrices for a wide range of electronic states. A pilot application of the new technique to construct averaged atomic natural orbital (ANO) basis sets for fully relativistic electronic structure calculations is presented.
An efficient approach to calculate approximate pure-state and transition reduced density matrices in the framework of the multireference relativistic Fock-space coupled cluster (FS CC) theory is proposed. The method is based on the effective operator formalism and consists of the direct substitution of the FS CC Ansatz for a wave operator into the effective operator expression with the subsequent truncation of expansion at the terms quadratic in cluster amplitudes. The final density matrix is defined by active-space density matrices of different ranks ‘dressed’ with contributions from cluster operators. The method gives a connected expression for pure-state density matrices, provided that the intermediate normalisation condition is fulfilled. Moreover, under some additional assumptions, the connectivity can also be ensured for calculated transition property matrix elements and natural transition spinors. The developed technique allows for fast and accurate calculations of one-particle reduced density matrices for a wide range of electronic states. A pilot application of the new technique to construct averaged atomic natural orbital (ANO) basis sets for fully relativistic electronic structure calculations is presented.
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.
The problem of finding the spectrum of the Fock-space Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{} for a system of many fermions is analyzed. The quasiparticle formalism is employed, with ``holes'' and ``particles'' defined with respect to some N-particle determinantal wave function (the model vacuum). The basic idea behind the proposed approach is to perform a similarity transformation of Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{}, such that the resulting effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{} is, unlike \ifmmode \hat{H}\else \^{H}\fi{}, a quasiparticle-number--conserving operator. It is shown that eigenvalues of \ifmmode \hat{G}\else \^{G}\fi{}, corresponding to small numbers of quasiparticles (0,1,2) can be easily calculated. This is equivalent to finding eigenvalues of \ifmmode \hat{H}\else \^{H}\fi{} for certain states of N, N\ifmmode\pm\else\textpm\fi{}1, and N\ifmmode\pm\else\textpm\fi{}2 particles. The construction of the operator transforming \ifmmode \hat{H}\else \^{H}\fi{} into \ifmmode \hat{G}\else \^{G}\fi{} (the wave operator) stems from an analysis of the structure of the algebra of linear operators acting in a (finite-dimensional) Fock space. The exponential Ansatz for the wave operator is used, resulting in a generalization of the coupled-cluster (CC) method of Coester [Nucl. Phys. 7, 421 (1958)]. The generalized CC equations determining the wave operator, and equations determining the effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{}, are presented in a diagrammatic form. An effort has been made to obtain a concise notation for expressing these equations in an algebraic form. Approximation schemes, necessary for practical applications of the proposed method, are also studied.
In this paper, we report the relativistic Fock space multireference coupled-cluster method for atomic structure calculations. We use the no-pair Dirac-Coulomb-Breit Hamiltonian, together with a finite B-spline basis set to expand the large and small components of the Dirac wave function. Our method is applied to calculate ionization energies, reduced matrix elements, lifetimes, and polarizabilities for many states of atomic francium. To evaluate uncertainties of our results and investigate the role of electron correlation effects, we carry out calculations using approximated models at different levels. The quality of our calculations is assessed by comparing with available experimental results, showing a good agreement. In addition, the tune-out wavelengths of the ground state in the range of 340--800 nm, the magic wavelengths for the transition $7s\ensuremath{-}8s$ in the range of 800--1500 nm and the transition $7s\ensuremath{-}7p$ in the range of 600--1500 nm are determined by evaluating the dynamic polarizabilities of the 7s,8s, and 7p states for a linearly polarized light. These tune-out and magic wavelengths may be useful for laser cooling and trapping of the Fr atom, and for related high-precision trapping measurements.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
A fully spin-adapted approach to many-electron density matrices is developed in the context of the unitary group approach to many-electron systems. An explicit expression for the single-electron spin-density operator, as a polynomial of degree two in the orbital U(n) generators, is derived for the case of spin-independent systems. Extensions to spin-dependent systems are also considered, leading to the appearance of total-spin transition densities, whose general properties are investigated. A corresponding formalism for the two-electron density matrix, which is capable of further generalization, is also developed. The results of this paper, together with recent developments on the matrix elements of the U(2n) generators in the electronic Gel’fand basis, afford a versatile method for the direct calculation of one- and two-body density matrices in the unitary group approach framework.
We report an extension of our previous development that incorporated quantum-chemical density matrix renormalization group (DMRG) into the complete active space second-order perturbation theory (CASPT2) [Y. Kurashige and T. Yanai, J. Chem. Phys. 135, 094104 (2011)]. In the previous study, the combined theory, referred to as DMRG-CASPT2, was built upon the use of pseudo-canonical molecular orbitals (PCMOs) for one-electron basis. Within the PCMO basis, the construction of the four-particle reduced density matrix (4-RDM) using DMRG can be greatly facilitated because of simplicity in the multiplication of 4-RDM and diagonal Fock matrix in the CASPT2 equation. In this work, we develop an approach to use more suited orbital basis in DMRG-CASPT2 calculations, e.g., localized molecular orbitals, in order to extend the domain of applicability. Because the multiplication of 4-RDM and generalized Fock matrix is no longer simple in general orbitals, an approximation is made to it using the cumulant reconstruction neglecting higher-particle cumulants. Also, we present the details of the algorithm to compute 3-RDM of the DMRG wavefunction as an extension of the 2-RDM algorithm of Zgid et al. [J. Chem. Phys. 128, 144115 (2008)] and Chan et al. [J. Chem. Phys. 128, 144117 (2008)]. The performance of the extended DMRG-CASPT2 approach was examined for large-scale multireference systems, such as low-lying excited states of long-chain polyenes and isomerization potential of {[Cu(NH3)3]2O2}(2+).