Dominika Zgid

The density matrix renormalization group self-consistent field method: Orbital optimization with the density matrix renormalization group method in the active space

We present the density matrix renormalization group self-consistent field (DMRG-SCF) approach that is analogous to the complete active space self-consisted field (CASSCF) method but instead of using for the description of the active space the full configuration interaction (FCI) method, the DMRG-SCF uses the density matrix renormalization group (DMRG) method. The DMRG-SCF approach, similarly to CASSCF, properly describes the multiconfigurational character of the wave function but avoids the exponential scaling of the FCI method and replaces it with a polynomial scaling. Hence, calculations for a larger number of orbitals and electrons in the active space are possible since the DMRG method provides an efficient tool to automatically select from the full Hilbert space the many-body contracted basis states that are the most important for the description of the wave function.

On the spin and symmetry adaptation of the density matrix renormalization group method

We present a spin-adapted density matrix renormalization group (DMRG) algorithm designed to target spin and spatial symmetry states that can be difficult to obtain while using a non-spin-adapted algorithm. The algorithmic modifications that have to be introduced into the usual density matrix renormalization group scheme in order to spin adapt it are discussed, and it is demonstrated that the introduced modifications do not change the overall scaling of the method. The new approach is tested on HNCO, a model system, that has a singlet-triplet curve crossing between states of the same symmetry. The advantages of the spin-adapted DMRG scheme are discussed, and it is concluded that the spin-adapted DMRG method converges better in almost all cases and gives more parallel curves to the full configuration interaction result than the non-spin-adapted method. It is shown that the spin-adapted DMRG energies can be lower than the ones obtained from the non-spin-adapted scheme. Such a counterintuitive result is explained by noting that the spin-adapted method is not a special case of the non-spin-adapted one; consequently, the spin-adapted result is not an upper bound for the non-spin-adapted energy.

Obtaining the two-body density matrix in the density matrix renormalization group method

We present an approach that allows to produce the two-body density matrix during the density matrix renormalization group (DMRG) run without an additional increase in the current disk and memory requirements. The computational cost of producing the two-body density matrix is proportional to O(M3k2+M2k4). The method is based on the assumption that different elements of the two-body density matrix can be calculated during different steps of a sweep. Hence, it is desirable that the wave function at the convergence does not change during a sweep. We discuss the theoretical structure of the wave function ansatz used in DMRG, concluding that during the one-site DMRG procedure, the energy and the wave function are converging monotonically at every step of the sweep. Thus, the one-site algorithm provides an opportunity to obtain the two-body density matrix free from the N-representability problem. We explain the problem of local minima that may be encountered in the DMRG calculations. We discuss theoretically why and when the one- and two-site DMRG procedures may get stuck in a metastable solution, and we list practical solutions helping the minimization to avoid the local minima.

Interaction of NH(XΣ−3) with He: Potential energy surface, bound states, and collisional Zeeman relaxation

A detailed analysis of the He-NH((3)Sigma(-)) van der Waals complex is presented. We discuss ab initio calculations of the potential energy surface and fitting procedures with relevance to cold collisions, and we present accurate calculations of bound energy levels of the triatomic complex as well as collisional properties of NH molecules in a buffer gas of (3)He. The influence of the external magnetic field used to trap the NH molecules and the effect of the atom-molecule interaction potential on the collisionally induced Zeeman relaxation are explored. It is shown that minute variations of the interaction potential due to different fitting procedures may alter the Zeeman relaxation rate at ultralow temperatures by as much as 50%.

Chapter 7 The Density Matrix Renormalization Group in Quantum Chemistry

The density matrix renormalization group (DMRG) is an electronic structure method that has recently been applied to ab initio quantum chemistry. Even at this early stage, it has enabled the solution of many problems that would previously have been intractable with any other method, in particular, multireference problems with very large active spaces. Here we provide an expository introduction to the theory behind the DMRG and give a brief overview of some recent applications and developments in the context of quantum chemistry.

Dynamical mean-field theory from a quantum chemical perspective

We investigate the dynamical mean-field theory (DMFT) from a quantum chemical perspective. Dynamical mean-field theory offers a formalism to extend quantum chemical methods for finite systems to infinite periodic problems within a local correlation approximation. In addition, quantum chemical techniques can be used to construct new ab initio Hamiltonians and impurity solvers for DMFT. Here, we explore some ways in which these things may be achieved. First, we present an informal overview of dynamical mean-field theory to connect to quantum chemical language. Next, we describe an implementation of dynamical mean-field theory where we start from an ab initio Hartree-Fock Hamiltonian that avoids double counting issues present in many applications of DMFT. We then explore the use of the configuration interaction hierarchy in DMFT as an approximate solver for the impurity problem. We also investigate some numerical issues of convergence within DMFT. Our studies are carried out in the context of the cubic hydrogen model, a simple but challenging test for correlation methods. Finally, we finish with some conclusions for future directions.

A study of cumulant approximations to n -electron valence multireference perturbation theory

We investigate the possibility of reducing the complexity of multireference perturbation theory through cumulant based approximations to the high-order density matrices that appear in such theories. Our test cases show that while the cumulant approximated forms are degraded in accuracy relative to the parent theory and exhibit intruder state problems that must be carefully handled, they may provide a route to a simple estimation of dynamic correlation when the parent perturbation theory is infeasible. Nonetheless, further work is clearly needed on better approximations to the denominators in the perturbation theory.

Communication: The description of strong correlation within self-consistent Green's function second-order perturbation theory

We report an implementation of self-consistent Greens function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Greens function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H32 finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Greens function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.

Communication: Towards ab initio self-energy embedding theory in quantum chemistry

The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Greens function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.

Systematically improvable multiscale solver for correlated electron systems

The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe how a hybrid between self-consistent second order perturbation theory and exact diagonalization can be used as a multi-scale solver for such systems. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the 2D Hubbard model, as a benchmark, we show that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the behavior of our method in multiple regimes of interaction strengths and doping of the model. The algorithm avoids difficulties such as double counting corrections, frequency dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Greens function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.