Source author record

Dominika Zgid

Dominika Zgid appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Relativistic Self-Consistent $GW$: Exact Two-Component Formalism with One-Electron Approximation for Solids

We present a formulation of relativistic self-consistent $GW$ for solids based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and the interplay of relativistic effects with electron correlation without adjustable parameters. Our all-electron implementation is fully $ab$ $initio$ and does not require a pseudopotential constructed from atomic calculations. We examine the effect of the X2C1e approximation by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables the construction of higher order theories, such as embedding theories, on top of perturbative calculations.

preprint2021arXiv

Dynamical Self-energy Mapping (DSEM) for quantum computing

For noisy intermediate-scale quantum (NISQ) devices only a moderate number of qubits with a limited coherence is available thus enabling only shallow circuits and a few time evolution steps in the currently performed quantum computations. Here, we present how to bypass this challenge in practical molecular chemistry simulations on NISQ devices by employing a classical-quantum hybrid algorithm allowing us to produce a sparse Hamiltonian which contains only $\mathcal{O}(n^2)$ terms in a Gaussian orbital basis when compared to the $\mathcal{O}(n^4)$ terms of a standard Hamiltonian, where $n$ is the number of orbitals in the system. Classical part of this hybrid entails parameterization of the sparse, fictitious Hamiltonian in such a way that it recovers the self-energy of the original molecular system. Quantum machine then uses this fictitious Hamiltonian to calculate the self-energy of the system. We show that the developed hybrid algorithm yields very good total energies for small molecular test cases while reducing the depth of the quantum circuit by at least an order of magnitude when compared with simulations involving a full Hamiltonian.

preprint2021arXiv

Exploring Coupled Cluster Green's function as a method for treating system and environment in Green's function embedding methods

Within the self-energy embedding theory (SEET) framework, we study coupled cluster Green's function (GFCC) method in two different contexts: as a method to treat either the system or environment present in the embedding construction. Our study reveals that when GFCC is used to treat the environment we do not see improvement in total energies in comparison to the coupled cluster method itself. To rationalize this puzzling result, we analyze the performance of GFCC as an impurity solver with a series of transition metal oxides. These studies shed light on strength and weaknesses of such a solver and demonstrate that such a solver gives very accurate results when the size of the impurity is small. We investigate if it is possible to achieve a systematic accuracy of the embedding solution when we increase the size of the impurity problem. We found that in such a case, the performance of the solver worsens, both in terms of finding the ground state solution of the impurity problem as well as the self-energies produced. We concluded that increasing the rank of GFCC solver is necessary to be able to enlarge impurity problems and achieve a reliable accuracy. We also have shown that natural orbitals from weakly correlated perturbative methods are better suited than symmetrized atomic orbitals (SAO) when the total energy of the system is the target quantity.

preprint2021arXiv

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green's functions obtained from the self-consistent solution of the Dyson equation can be employed in evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green's function methods, because of non-convex grand potential functional and non-idempotent density matrix. Moreover, the inclusion of correlation effects in the form of the self-energy matrix and changing chemical potential or fluctuations in the number of particles can make the optimization problem more difficult. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS (CDIIS) and LCIIS, and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. The generalized CDIIS and LCIIS methods successfully converged restricted GF2 calculations for a number of strongly correlated systems, which could not be converged before. We also provide practical recommendations to guide convergence in such pathological cases.

preprint2020arXiv

Ab-Initio self-energy embedding for the photoemission spectra of NiO and MnO

The accurate ab-initio simulation of periodic solids with strong correlations is one of the grand challenges of condensed matter. While mature methods exist for weakly correlated solids, the ab-initio description of strongly correlated systems is an active field of research. In this work, we show results for the single particle spectral function of the two correlated $d$-electron solids NiO and MnO from self-energy embedding theory. Unlike earlier work, the theory does not use any adjustable parameters and is fully ab-initio, while being able to treat both the strong correlation and the non-local screening physics of these materials. We derive the method, discuss aspects of the embedding and choices of physically important orbitals, and compare our results to x-ray and angle-resolved photoemission spectroscopy as well as bremsstrahlung-isochromat spectroscopy.

preprint2020arXiv

Error propagation in the fully self-consistent stochastic second-order Green's function method

We present an implementation of a fully self-consistent finite temperature second order Green's function perturbation theory (GF2) within the diagrammatic Monte Carlo framework. In contrast to the previous implementations of stochastic GF2 ({\it J. Chem. Phys.},{\bf 151}, 044144 (2019)), the current self-consistent stochastic GF2 does not introduce a systematic bias of the resulting electronic energies. Instead, the introduced implementation accounts for the stochastic errors appearing during the solution of the Dyson equation. We present an extensive discussion of the error handling necessary in a self-consistent procedure resulting in dressed Green's function lines. We test our method on a series of simple molecular examples.

preprint2020arXiv

Legendre-spectral Dyson equation solver with super-exponential convergence

Quantum many-body systems in thermal equilibrium can be described by the imaginary time Green's function formalism. However, the treatment of large molecular or solid ab inito problems with a fully realistic Hamiltonian in large basis sets is hampered by the storage of the Green's function and the precision of the solution of the Dyson equation. We present a Legendre-spectral algorithm for solving the Dyson equation that addresses both of these issues. By formulating the algorithm in Legendre coefficient space, our method inherits the known faster-than-exponential convergence of the Green's function's Legendre series expansion. In this basis, the fast recursive method for Legendre polynomial convolution, enables us to develop a Dyson equation solver with quadratic scaling. We present benchmarks of the algorithm by computing the dissociation energy of the helium dimer He$_2$ within dressed second-order perturbation theory. For this system, the application of the Legendre spectral algorithm allows us to achieve an energy accuracy of $10^{-9} E_h$ with only a few hundred expansion coefficients.

preprint2019arXiv

Direct comparison of many-body methods for realistic electronic Hamiltonians

A large collaboration carefully benchmarks 20 first principles many-body electronic structure methods on a test set of 7 transition metal atoms, and their ions and monoxides. Good agreement is attained between the 3 systematically converged methods, resulting in experiment-free reference values. These reference values are used to assess the accuracy of modern emerging and scalable approaches to the many-electron problem. The most accurate methods obtain energies indistinguishable from experimental results, with the agreement mainly limited by the experimental uncertainties. Comparison between methods enables a unique perspective on calculations of many-body systems of electrons.

preprint2016arXiv

Efficient temperature-dependent Green's function methods for realistic systems: using cubic spline interpolation to approximate Matsubara Green's functions

The popular, stable, robust and computationally inexpensive cubic spline interpolation algorithm is adopted and used for finite temperature Green's function calculations of realistic systems. We demonstrate that with appropriate modifications the temperature dependence can be preserved while the Green's function grid size can be reduced by about two orders of magnitude by replacing the standard Matsubara frequency grid with a sparser grid and a set of interpolation coefficients. We benchmarked the accuracy of our algorithm as a function of a single parameter sensitive to the shape of the Green's function. Through numerous examples, we confirmed that our algorithm can be utilized in a systematically improvable, controlled, and black-box manner and highly accurate one- and two-body energies and one-particle density matrices can be obtained using only around 5% of the original grid points. Additionally, we established that to improve accuracy by an order of magnitude, the number of grid points needs to be doubled, whereas for the Matsubara frequency grid an order of magnitude more grid points must be used. This suggests that realistic calculations with large basis sets that were previously out of reach because they required enormous grid sizes may now become feasible.

preprint2016arXiv

Exploring connections between statistical mechanics and Green's functions for realistic systems. Temperature dependent electronic entropy and internal energy from a self-consistent second-order Green's function

Including finite-temperature effects from the electronic degrees of freedom in electronic structure calculations of semiconductors and metals is desired; however, in practice it remains exceedingly difficult when using zero-temperature methods, since these methods require an explicit evaluation of multiple excited states in order to account for any finite-temperature effects. Using a Matsubara Green's function formalism remains a viable alternative, since in this formalism it is easier to include thermal effects and to connect the dynamic quantities such as the self-energy with static thermodynamic quantities such as the Helmholtz energy, entropy, and internal energy. However, despite the promising properties of this formalism, little is know about the multiple solutions of the non-linear equations present in the self-consistent Matsubara formalism and only a few cases involving a full Coulomb Hamiltonian were investigated in the past. Here, to shed some light onto the iterative nature of the Green's function solutions, we self-consistently evaluate the thermodynamic quantities for a one-dimensional (1D) hydrogen solid at various interatomic separations and temperatures using the self-energy approximated to second-order (GF2). At many points in the phase diagram of this system, multiple phases such as a metal and an insulator exist, and we are able to determine the most stable phase from the analysis of Helmholtz energies. Additionally, we show the evolution of the spectrum of 1D boron nitride (BN) to demonstrate that GF2 is capable of qualitatively describing the temperature effects influencing the size of the band gap.

preprint2016arXiv

Rigorous ab initio quantum embedding for quantum chemistry using Green's function theory: screened interaction, non-local self-energy relaxation, orbital basis, and chemical accuracy

We present a detailed discussion of self-energy embedding theory (SEET) which is a quantum embedding scheme allowing us to describe a chosen subsystem very accurately while keeping the description of the environment at a lower cost. We apply SEET to molecular examples where commonly our chosen subsystem is made out of a set of strongly correlated orbitals while the weakly correlated orbitals constitute an environment. Such a self-energy separation is very general and to make this procedure applicable to multiple systems a detailed and practical procedure for the evaluation of the system and environment self-energy is necessary. We list all the intricacies for one of the possible procedures while focusing our discussion on many practical implementation aspects such as the choice of best orbital basis, impurity solver, and many steps necessary to reach chemical accuracy. Finally, on a set of carefully chosen molecular examples, we demonstrate that SEET which is a controlled, systematically improvable Green's function method can be as accurate as established wavefunction quantum chemistry methods.

preprint2016arXiv

Stochastic self-consistent Green's function second-order perturbation theory (sGF2)

The second-order Green's function method (GF2) was shown recently to be an accurate self-consistent approach for electronic structure of correlated systems since the self-energy accounts for both the weak and some of the strong correlation. The numerical scaling of GF2 is quite steep however, $O({N^5})$ (where the pre-factor is often hundreds), effectively preventing its application to large systems. Here, we develop a stochastic approach to GF2 (sGF2) where the self-energy is evaluated by a random-vector decomposition of Green's functions so that the dominant part of the calculation scales quasi linearly with system size. A study of hydrogen chains shows that the resulting approach is numerically efficient and accurate, as the stochastic errors are very small, 0.05% of the correlation energy for large systems with only a moderate computational effort. The method also yields automatically efficient MP2 energies and is automatically temperature dependent.

preprint2015arXiv

Efficient temperature-dependent Green's functions methods for realistic systems: compact grids for orthogonal polynomial transforms

The temperature-dependent Matsubara Green's function that is used to describe temperature-dependent behavior is expressed on a numerical grid. While such a grid usually has a couple of hundred points for low-energy model systems, for realistic systems in large basis sets the size of an accurate grid can be tens of thousands of points, constituting a severe computational and memory bottleneck. In this paper, we determine efficient imaginary time grids for the temperature-dependent Matsubara Green's function formalism that can be used for calculations on realistic systems. We show that due to the use of orthogonal polynomial transform, we can restrict the imaginary time grid to few hundred points and reach micro-Hartree accuracy in the electronic energy evaluation. Moreover, we show that only a limited number of orthogonal polynomial expansion coefficients are necessary to preserve accuracy when working with a dual representation of Green's function or self-energy and transforming between the imaginary time and Matsubara frequency domain.

preprint2015arXiv

Fractional charge and spin errors in self-consistent Green's function theory

We examine fractional charge and spin errors in self-consistent Green's function theory within a second-order approximation (GF2). For GF2 it is known that the summation of diagrams resulting from the self-consistent solution of the Dyson equation removes the divergences pathological to second-order Moller-Plesset theory (MP2) for strong correlations. In the language often used in density functional theory contexts, this means GF2 has a greatly reduced fractional spin error relative to MP2. The natural question then is what effect, if any, does the Dyson summation have on the fractional charge error in GF2? To this end we generalize our previous implementation of GF2 to open-shell systems and analyze its fractional spin and charge errors. We find that like MP2, GF2 possesses only a very small fractional charge error, and consequently little many electron self-interaction error. This shows that GF2 improves on the critical failings of MP2, but without altering the positive features that make it desirable. Furthermore, we find that GF2 has both less fractional charge and fractional spin errors than typical hybrid density functionals as well as random phase approximation with exchange.

preprint2015arXiv

Ionization potentials and electron affinities from the extended Koopmans' theorem in self-consistent Green's function theory

One-body Green's function theories implemented on the real frequency axis offer a natural formalism for the unbiased theoretical determination of quasiparticle spectra in molecules and solids. Self-consistent Green's function methods employing the imaginary axis formalism on the other hand can benefit from the iterative implicit resummation of higher order diagrams that are not included when only the first iteration is performed. Unfortunately, the imaginary axis Green's function does not give direct access to the desired quasiparticle spectra, which undermines its utility. To this end we investigate how reliably one can calculate quasiparticle spectra from the Extended Koopmans' Theorem (EKT) applied to the imaginary time Green's function in a second order approximation (GF2). We find that EKT in conjunction with GF2 yields IPs and EAs that systematically underestimate experimental and accurate coupled-cluster reference values for a variety of molecules and atoms. This establishes that the EKT allows one to utilize the computational advantages of an imaginary axis implementation, while still being able to acquire real axis spectral properties. Because the EKT requires negligible computational effort, and can be used with a Green's function from any level of theory, we conclude that it is a potentially very useful tool for the systematic study of quasiparticle spectra in realistic systems.

preprint2015arXiv

Second-order Green's function perturbation theory for periodic systems

Despite recent advances, systematic quantitative treatment of the electron correlation problem in extended systems remains a formidable task. Systematically improvable Green's function methods capable of quantitatively describing weak and at least qualitatively strong correlations appear promising candidates for computational treatment of periodic systems. We present a periodic implementation of temperature-dependent self-consistent 2nd-order Green's function method (GF2), where the self-energy is evaluated in the basis of atomic orbitals. Evaluating the real-space self-energy in atomic orbitals and solving the Dyson equation in $\mathbf{k}$-space are the key components of a computationally feasible algorithm. We apply this technique to the 1D hydrogen lattice - a prototypical crystalline system with a realistic Hamiltonian. By analyzing the behavior of the spectral functions, natural occupations, and self-energies, we claim that GF2 is able to recover metallic, band insulating, and at least qualitatively Mott regimes. We observe that the iterative nature of GF2 is essential to the emergence of the metallic and Mott phases.

preprint2015arXiv

Systematically improvable multi-scale 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 Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.

preprint2015arXiv

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 Green's function (GF2) method is used to describe the non-local correlations, while the full configuration interaction (FCI) 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 (NEVPT2) 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.

preprint2014arXiv

Local Hamiltonians for quantitative Green's function embedding methods

Embedding calculations that find approximate solutions to the Schrödinger equation for large molecules and realistic solids are performed commonly in a three step procedure involving (i) construction of a model system with effective interactions approximating the low energy physics of the initial realistic system, (ii) mapping the model system onto an impurity Hamiltonian, and (iii) solving the impurity problem. We have developed a novel procedure for parametrizing the impurity Hamiltonian that avoids the mathematically uncontrolled step of constructing the low energy model system. Instead, the impurity Hamiltonian is immediately parametrized to recover the self-energy of the realistic system in the limit of high frequencies or short time. The effective interactions parametrizing the fictitious impurity Hamiltonian are local to the embedded regions, and include all the non-local interactions present in the original realistic Hamiltonian in an implicit way. We show that this impurity Hamiltonian can lead to excellent total energies and self-energies that approximate the quantities of the initial realistic system very well. Moreover, we show that as long as the effective impurity Hamiltonian parametrization is designed to recover the self-energy of the initial realistic system for high frequencies, we can expect a good total energy and self-energy. Finally, we propose two practical ways of evaluating effective integrals for parametrizing impurity models.

preprint2014arXiv

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

We report an implementation of self-consistent Green's 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 Green's 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 H$_{32}$ 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 Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.

preprint2012arXiv

Truncated Configuration Interaction expansions as solvers for correlated quantum impurity models and dynamical mean field theory

The development of polynomial cost solvers for correlated quantum impurity models, with controllable errors, is a central challenge in quantum many-body physics, where these models find applications ranging from nano-science to the dynamical mean-field theory (DMFT). Here we describe how configuration interaction (CI) approximations to exact diagonalization (ED) may be used as solvers in DMFT. CI approximations retain the main advantages of ED, such as the ability to treat general interactions and off-diagonal hybridizations and to obtain real spectral information, but are of polynomial cost. Furthermore, their errors can be controlled by monitoring the convergence of physical quantities as a function of the CI hierarchy. Using benchmark DMFT applications, such as single-site DMFT of the 1D Hubbard model and $2\times 2$ cluster DMFT of the 2D Hubbard model, we show that CI approximations allow us to obtain near-exact ED results for a tiny fraction of the cost. This is true over the entire range of interaction strengths including "difficult" regimes, such as in the pseudogap phase of the 2D Hubbard model. We use the ability of CI approximations to treat large numbers of orbitals to demonstrate convergence of the bath representation in the $2\times 2$ cluster DMFT using a 24 bath orbital representation. CI approximations thus form a promising route to extend ED to problems that are currently difficult to study using other solvers such as continuous-time quantum Monte Carlo, including impurity models with large numbers of orbitals and general interactions.

preprint2010arXiv

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.

preprint2010arXiv

Embedding theory for excited states with inclusion of self-consistent environment screening

We present a general embedding theory of electronic excitations of a relatively small, localized system in contact with an extended, chemically complex environment. We demonstrate how to include the screening response of the environment into highly accurate electronic structure calculation of the localized system by means of an effective interaction between the electrons, which contains only screening processes occurring in the environment. For the common case of a localized system which constitutes an inhomogeneity in an otherwise homogeneous system, such as a defect in a crystal, we show how matrix elements of the environment-screened interaction can be calculated from density-functional calculations of the homogeneous system only. We apply our embedding theory to the calculation of excitation energies in crystalline ethylene.