Abstract:The Green's function coupled cluster (GFCC) method, originally proposed in the early 1990s, is a powerful many-body tool for computing and analyzing the electronic structure of molecular and periodic systems, especially when electrons of the system are strongly correlated. However, in order for the GFCC to become a method that may be routinely used in the electronic structure calculations, robust numerical techniques and approximations must be employed to reduce its extremely high computational overhead. In our recent studies, it has been demonstrated that the GFCC equations can be solved directly in the frequency domain using iterative linear solvers, which can be easily distributed in a massively parallel environment. In the present work, we demonstrate a successful application of model-order-reduction (MOR) techniques in the GFCC framework. Briefly speaking, for a frequency regime of interest that requires high-resolution descriptions of spectral function, instead of solving the GFCC linear equation of full dimension for every single frequency point of interest, an efficiently solvable linear system model of a reduced dimension may be built upon projecting the original GFCC linear system onto a subspace. From this reduced order model is obtained a reasonable approximation to the full dimensional GFCC linear equations in both interpolative and extrapolative spectral regions. Here, we show that the subspace can be properly constructed in an iterative manner from the auxiliary vectors of the GFCC linear equations at some selected frequencies within the spectral region of interest. During the iterations, the quality of the subspace, as well as the linear system model, can be systematically improved. The method is tested in this work in terms of the efficiency and accuracy of computing spectral functions for some typical molecular systems such as carbon monoxide, 1,3-butadiene, benzene, and adenine. To reach the same level of accuracy as that of the original GFCC method, the application of MOR in the GFCC method is able to significantly lower the original computational cost for the aforementioned molecules in designated frequency regimes. As a byproduct, the reduced order model obtained by this method is found to provide a high-quality initial guess, which improves the convergence rate for the existing iterative linear solver.

An O(N^3) implementation of Hedin's GW approximation for molecules

Scalable Molecular GW Calculations: Valence and Core Spectra

Low-order Scaling $G_0W_0$ by Pair Atomic Density Fitting

Low-Order Scaling G0W0 by Pair Atomic Density Fitting

Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions

Two-component $GW$ calculations: Cubic scaling implementation and comparison of vertex corrected and partially self-consistent $GW$ variants

Accelerating Time-Dependent Density Functional Theory and GW Calculations for Molecules and Nanoclusters with Symmetry Adapted Interpolative Separable Density Fitting

Realizing Effective Cubic-Scaling Coulomb Hole Plus Screened Exchange Approximation in Periodic Systems Via Interpolative Separable Density Fitting with a Plane-Wave Basis Set.

Scalable GW software for quasiparticle properties using OpenAtom

A Density-Based Basis-Set Incompleteness Correction for GW Methods

GW with hybrid functionals for large molecular systems

Exploiting a derivative discontinuity estimate for accurate $G_0W_0$ ionization potentials and electron affinities

Solving multi-pole challenges in the GW100 benchmark enables precise low-scaling GW calculations

Approximate Green's Function Coupled Cluster Method Employing Effective Dimension Reduction.

The $GW$ Approximation: A Quantum Chemistry Perspective

Fully selfconsistent GW calculations for molecules

Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules

Molgw 1: Many-body Perturbation Theory Software for Atoms, Molecules, and Clusters

All-electron periodic G 0 W 0 implementation with numerical atomic orbital basis functions: Algorithm and benchmarks

Quasi-Particle Self-Consistent $GW$ for Molecules

Accelerating Analytic-Continuation GW Calculations with a Laplace Transformation and Natural Auxiliary Functions