Ivan Duchemin
email: ivan.duchemin@cea.fr

Xavier Blase

Robust analytic continuation approach to many-body GW calculations

The analytic continuation of the GW self-energy from the imaginary to the real energy axis is a central difficulty for approaches exploiting the favourable properties of response functions at imaginary frequencies. Within a scheme merging contourdeformation and analytic-continuation techniques, we show on the basis of extensive calculations for large molecular sets that it is preferable to perform an analytic continuation of the dynamically screened Coulomb potential W rather than the much more structured self-energy operator. The case of states lying far away from the gap, including core states, is addressed by generalizing the analytic continuation scheme, accounting further for quasiparticle lifetimes.

Introduction

The GW approximation [START_REF] Hedin | New Method for Calculating the One-Particle Green's Function with Application to the Electron-Gas Problem[END_REF][START_REF] Strinati | Dynamical Correlation Effects on the Quasiparticle Bloch States of a Covalent Crystal[END_REF][START_REF] Hybertsen | Electron correlation in semiconductors and insulators: Band gaps and quasiparticle energies[END_REF][START_REF] Godby | Self-energy operators and exchange-correlation potentials in semiconductors[END_REF][START_REF] Farid | Electron Correlation in the Solid State -Chapter 3[END_REF][START_REF] Aryasetiawan | [END_REF][START_REF] Onida | Electronic excitations: density-functional versus many-body Green's-function approaches[END_REF][START_REF] Martin | Interacting Electrons: Theory and Computational Approaches[END_REF] to the exchange-correlation self-energy stands nowadays as a very popular formalism to calculate charged excitations, as obtained e.g. using photoemission techniques, in metallic or semiconducting, finite size or periodic systems. Belonging to the family of Green's function many-body perturbation theories (MBPT), the related 1 arXiv:1912.06459v1 [cond-mat.mtrl-sci] 13 Dec 2019 formalism dates back to the mid-60s with a pioneering application to the electronic properties of the interacting homogeneous electron gas. [START_REF] Hedin | New Method for Calculating the One-Particle Green's Function with Application to the Electron-Gas Problem[END_REF] Nowadays, efficient implementations, including algorithms with cubic [START_REF] Foerster | An O(N3) implementation of Hedin's GW approximation for molecules[END_REF][START_REF] Liu | Cubic scaling GW : Towards fast quasiparticle calculations[END_REF][START_REF] Wilhelm | Toward GW Calculations on Thousands of Atoms[END_REF] or sub-cubic [START_REF] Neuhauser | Breaking the Theoretical Scaling Limit for Predicting Quasiparticle Energies: The Stochastic GW Approach[END_REF] scaling with system size, paved the way to applications to large systems with up to several hundred atoms. [START_REF] Wilhelm | Toward GW Calculations on Thousands of Atoms[END_REF][START_REF] Neuhauser | Breaking the Theoretical Scaling Limit for Predicting Quasiparticle Energies: The Stochastic GW Approach[END_REF][START_REF] Govoni | Large Scale GW Calculations[END_REF][START_REF] Li | Combining the Many-Body GW Formalism with Classical Polarizable Models: Insights on the Electronic Structure of Molecular Solids[END_REF][START_REF] Gao | Speeding up GW Calculations to Meet the Challenge of Large Scale Quasiparticle Predictions[END_REF] Besides providing accurate electronic energy levels, the GW formalism has been used further to evaluate the quasiparticle lifetime with respect to electron-electron scattering [START_REF] Campillo | Inelastic Lifetimes of Hot Electrons in Real Metals[END_REF][START_REF] Keyling | Comparison of the lifetime of excited electrons in noble metals[END_REF][START_REF] Marini | Quasiparticle band-structure effects on the d hole lifetimes of copper within the GW approximation[END_REF][START_REF] Yi | Quasiparticle band structures and lifetimes in noble metals using Gaussian orbital basis sets[END_REF] or provide more accurate electron-phonon coupling matrix elements as compared to DFT calculations performed with standard semilocal functionals. [START_REF] Lazzeri | Impact of the electron-electron correlation on phonon dispersion: Failure of LDA and GGA DFT functionals in graphene and graphite[END_REF][START_REF] Faber | Electron-phonon coupling in the C 60 fullerene within the many-body GW approach[END_REF][START_REF] Yin | Correlation-Enhanced Electron-Phonon Coupling: Applications of GW and Screened Hybrid Functional to Bismuthates, Chloronitrides, and Other High-T c Superconductors[END_REF][START_REF] Li | Electron-Phonon Coupling from Ab Initio Linear-Response Theory within the GW Method: Correlation-Enhanced Interactions and Superconductivity in Ba 1-x K x BiO 3[END_REF] At the core of Green's function MBPT stands the self-energy operator Σ(r, r ; E) that accounts for exchange and correlation interactions beyond the classical Hartree term. An important feature of such an operator, as compared e.g. to the DFT exchange correlation or the Fock exact exchange potentials, lies in its energy dependence. In short, the self-energy operator must be taken at the energy of the state it is acting on. Considering the specific case of the GW approximation, the self-energy reads: Σ(r, r ; E) = i 2π dω e iηω G(r, r ; E + ω)W (r, r ; ω)

where W is the time-ordered screened-Coulomb potential and G the time-ordered singleparticle Green's function (η is a small positive infinitesimal). The treatment of the energy integration involved in the construction of the GW self-energy operator stands as a central difficulty as it requires, in principle, the knowledge of the dynamically screened Coulomb potential along the real frequency axis where G and W present a large number of poles.

Such an integral can be performed analytically provided that one knows the exact spectral representation of the screened Coulomb operator. This representation can be obtained by calculating explicitly the poles and eigenvectors of the susceptibility operator, within e.g.

the random phase approximation (RPA), at the cost of diagonalizing a matrix of which the size grows quadratically with system size, namely a formally O(N 6 ) process. [START_REF] Tiago | Optical excitations in organic molecules, clusters, and defects studied by first-principles Green's function methods[END_REF][START_REF] Van Setten | The GW-Method for Quantum Chemistry Applications: Theory and Implementation[END_REF][START_REF] Bruneval | molgw 1: Many-body perturbation theory software for atoms, molecules, and clusters[END_REF][START_REF] Véril | Unphysical Discontinuities in GW Methods[END_REF] As a simplified but related technique, plasmon-pole approximations [START_REF] Hybertsen | Electron correlation in semiconductors and insulators: Band gaps and quasiparticle energies[END_REF][START_REF] Von Der Linden | Precise quasiparticle energies and Hartree-Fock bands of semiconductors and insulators[END_REF][START_REF] Godby | Metal-insulator transition in Kohn-Sham theory and quasiparticle theory[END_REF][START_REF] Engel | Generalized plasmon-pole model and plasmon band structures of crystals[END_REF][START_REF] Blase | Self-energy effects on the surface-state energies of H-Si(111)[END_REF][START_REF] Rohlfing | Efficient scheme for GW quasiparticle bandstructure calculations with applications to bulk Si and to the Si(001)-(2×1) surface[END_REF] provide a much simplified explicit spectral representation that has been used in the early days of the GW formalism. [START_REF] Hybertsen | Electron correlation in semiconductors and insulators: Band gaps and quasiparticle energies[END_REF][START_REF] Godby | Metal-insulator transition in Kohn-Sham theory and quasiparticle theory[END_REF] While plasmon-pole approximations are extremely efficient numerically, their accuracy can be challenged in systems where response functions are not clearly dominated by well defined plasmon peaks. This is the case e.g. in small molecular systems where discrete electron-hole excitations may not merge into broad collective modes.

Numerical integrations by explicitly calculating the susceptibility on discretized frequency grids provide alternative solutions. While integration along the real axis has been implemented in a few codes, [START_REF] Miyake | Efficient algorithm for calculating noninteracting frequency-dependent linear response functions[END_REF][START_REF] Shishkin | Implementation and performance of the frequency-dependent GW method within the PAW framework[END_REF][START_REF] Rostgaard | Fully self-consistent GW calculations for molecules[END_REF][START_REF] Liu | Numerical integration for ab initio many-electron self energy calculations within the GW approximation[END_REF] calculation of W along the imaginary axis -where all operators present a much smoother frequency dependence -has been also proposed in the very early days of GW calculations. [START_REF] Godby | Self-energy operators and exchange-correlation potentials in semiconductors[END_REF] Further, the so-called space-time technique, [START_REF] Wilhelm | Toward GW Calculations on Thousands of Atoms[END_REF][START_REF] Rojas | Space-Time Method for Ab Initio Calculations of Self-Energies and Dielectric Response Functions of Solids[END_REF][START_REF] Umari | Optimal representation of the polarization propagator for large-scale GW calculations[END_REF] with its recent extensions in quantum chemistry within the framework of interpolative density fitting or separable resolution-of-the-identity approaches, [START_REF] Lu | Compression of the electron repulsion integral tensor in tensor hypercontraction format with cubic scaling cost[END_REF][START_REF] Lu | Cubic scaling algorithms for RPA correlation using interpolative separable density fitting[END_REF][START_REF] Duchemin | Separable resolution-of-the-identity with all-electron Gaussian bases: Application to cubic-scaling RPA[END_REF] leads to the calculation of the susceptibility at imaginary times. Using Laplace transform, the susceptibility and related screened Coulomb potential can be then obtained along the imaginary-frequency axis.

From the knowledge of the screened Coulomb potential along the imaginary axis, one may calculate the self-energy operator at imaginary energies and proceed with its analytic continuation (AC) to the real axis where it is needed. [START_REF] Wilhelm | Toward GW Calculations on Thousands of Atoms[END_REF][START_REF] Rojas | Space-Time Method for Ab Initio Calculations of Self-Energies and Dielectric Response Functions of Solids[END_REF][42][START_REF] Umari | GW quasiparticle spectra from occupied states only[END_REF][START_REF] Nguyen | Improving accuracy and efficiency of calculations of photoemission spectra within the many-body perturbation theory[END_REF][START_REF] Pham | GW calculations using the spectral decomposition of the dielectric matrix: Verification, validation, and comparison of methods[END_REF][START_REF] Chu | All-electron self-consistent GW in the Matsubara-time domain: Implementation and benchmarks of semiconductors and insulators[END_REF][START_REF] Nabok | Accurate all-electron G 0 W 0 quasiparticle energies employing the full-potential augmented plane-wave method[END_REF] Such an approach is extremely elegant and numerically efficient, but difficulties exist concerning its accuracy in relation with the pole structure of the GW operator at real frequencies. [START_REF] Cancés | A mathematical analysis of the GW 0 method for computing electronic excited energies of molecules[END_REF] In a recent benchmark study of small molecular systems, [START_REF] Van Setten | Benchmarking G0W0 for Molecular Systems[END_REF] the need to increase the number of Padé approximants to a large number was documented for several systems. Even with more than a hundred approximants, pathological systems could not be treated with the desired accuracy in the most simple case of frontier orbitals.

Finally, in the contour deformation technique, [START_REF] Godby | Self-energy operators and exchange-correlation potentials in semiconductors[END_REF][START_REF] Govoni | Large Scale GW Calculations[END_REF][START_REF] Farid | GW approach to the calculation of electron self-energies in semiconductors[END_REF][START_REF] Lebègue | Implementation of an all-electron GW approximation based on the projector augmented wave method without plasmon pole approximation: Application to Si, SiC, AlAs, InAs, NaH, and KH[END_REF][START_REF] Blase | First-principles GW calculations for fullerenes, porphyrins, phtalocyanine, and other molecules of interest for organic photovoltaic applications[END_REF][START_REF] Golze | Core-Level Binding Energies from GW: An Efficient Full-Frequency Approach within a Localized Basis[END_REF] the integral along the real-axis is transformed into an integral along the imaginary axis. This integral must be complemented however by the contribution of the residues of W at the poles of the Green's function G that are shifted into the first/third quadrant by the energy at which the self-energy is required (see Fig. 1). The contour-deformation approach requires thus the knowledge of the screened Coulomb potential along both the imaginary and the real axis.

In the present study, we explore an approach merging the contour deformation technique with an analytic continuation of the screened Coulomb potential W to the complex plane for the evaluation of the needed residues. This leads to a very robust scheme as compared to the analytic continuation of the much more structured self-energy Σ. Our approach is validated by exploring its performances over the GW 100 test set [START_REF] Van Setten | Benchmarking G0W0 for Molecular Systems[END_REF][START_REF] Krause | Coupled-cluster reference values for the GW27 and GW100 test sets for the assessment of GW methods[END_REF][START_REF] Caruso | Benchmark of GW Approaches for the GW100 Test Set[END_REF][START_REF] Maggio | A Plane Wave Perspective for Small Molecules[END_REF][START_REF] Govoni | Comparison of Methods and Accuracy of Results Obtained with the WEST Code[END_REF] and a recent set of medium size molecules, [START_REF] Knight | Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods[END_REF] including systems that were shown to be pathological for the direct analytic continuation of the self-energy. We show that very few calculations of the susceptibility along the imaginary axis are required to reach the meV accuracy in the calculation of the quasiparticle energy of states lying close to the gap, in particular the highest occupied (HOMO) and lowest unoccupied (LUMO) molecular orbitals. Further, the GW self-energy and associated spectral functions can be efficiently and accurately obtained over a large energy range, covering e.g. the full valence manifold, provided that (a) a few additional reference W (ω) matrices are calculated over a coarse energy grid in the complex plane with (ω) spanning the targeted energy window and (b) the lifetime of quasiparticle states is accounted for. This approach is finally adapted to the case of core levels. The present study extends to the GW formalism applied to large sets of molecular systems the recent scheme proposed in the GT study of solid iron, [START_REF] Friedrich | Tetrahedron integration method for strongly varying functions: Application to the GT self-energy[END_REF] generalizing further the hybrid contour-deformation/analytic-continuation scheme to deep lying states.

Theory

We briefly recall here below the main characteristics of the general GW formalism as well as the contour deformation approach that we combine with an analytic continuation (AC) of the screened Coulomb potential. Details of the integration grid along the imaginary axis are provided, together with the technical details associated with the calculations presented in the Section "Results".

The GW formalism

In a standard GW calculation, one starts from an input time-ordered one-body Green's function typically built from {ε n , φ n } Kohn-Sham eigenstates, namely:

G(r, r ; ω) = n φ n (r)φ * n (r ) ω -ε n + iη × sgn(ε n -E F ) ( 2
)
where η is a small positive infinitesimal and E F the Fermi energy. The corresponding GW self-energy Σ(r, r ; E) is then generally divided into exact exchange Σ X (r, r ; E) and correlation-only Σ C (r, r ; E) contributions. The exact exchange operator built from the onebody eigenstates can be formulated as:

Σ X (r, r ; E) = i 2π +∞ -∞ dω e iηω G(r, r ; E + ω)V (r, r ) (3)
with V the bare Coulomb potential, so that by writing W = W -V , the correlation-only self-energy operator becomes:

Σ C (r, r ; E) = i 2π +∞ -∞ dω e iηω G(r, r ; E + ω) W (r, r ; ω) (4)
In order to highlight the desirable features of the present approach, one can take on general grounds the (generally unknown) spectral representation of W :

W (r, r ; ω) = λ w λ (r)w * λ (r ) ω -Ω λ + iη - w λ (r)w * λ (r ) ω + Ω λ -iη (5)
with {Ω λ , w λ (r)} the related (positive) pole energies and amplitudes, indexed by λ. Together with the spectral representation of G, the knowledge of the spectral representation of W allows the analytical integration of Eq. 4, yielding for the expression of Σ C (r, r ; E):

Σ C (r, r ; E) = i,λ φ i (r)w λ (r)w * λ (r )φ * i (r ) E -ε i + Ω λ -iη + a,λ φ a (r)w λ (r)w * λ (r )φ * a (r ) E -ε a -Ω λ + iη (6)
where i/a index occupied/unoccupied input (e.g. Kohn-Sham) eigenstates. As such, the resulting Σ C (r, r ; E) exhibits numerous poles along the real axis, distributed at energies {ε i -Ω λ } and {ε a + Ω λ }, summing up to the number of poles of G times the number of poles of W . Such a large number of poles may lead to potential difficulties with the standard AC approach that attempts to obtain Σ C (r, r ; E) from the Σ C (r, r ; iω) calculated along the imaginary frequency axis. Contrary to the input Green's function G, the spectral representation of W is generally unknown and one has to perform integration of Eq. 4 numerically. As a way to proceed, the contour-deformation approach transforms this latter integral along the real-frequency axis as an integral along the imaginary-frequency axis plus residues associated with the poles of G(E + ω) that have entered the first and third quadrant (see Fig. 1), namely:

The contour-deformation approach

+ + + + + + + + + + + + +
Σ C (r, r ; E) = - 1 2π +∞ -∞ dω G(r, r ; E + iω) W (r, r ; iω) - i φ i (r)φ * i (r ) W (r, r ; ε i -E)θ(ε i -E) + a φ a (r)φ * a (r ) W (r, r ; E -ε a )θ(E -ε a ) (7)
with θ the Heaviside function and ε i /ε a the occupied/unoccupied eigenstates used to build G. This reformulation is exact, leading to accurate results provided that the numerical treatment of the frequency integral is converged. In fact, excellent accuracy can be achieved with few quadrature points due to the absence of poles of W along the imaginary axis.

The main drawback of the contour-deformation approach is that the explicit calculation of the required W (r, r ; ε n -E) residues can become a bottleneck if the energy E at which the self-energy is required lands outside of the HOMO-LUMO gap, yielding potentially a large number of poles of G(E + ω) in the first/third quadrant. Indeed, the number of Kohn-Sham eigenstates within a given energy window grows linearly with system size, and so does the number of associated residues. Assuming a standard calculation of the independent-electron operator χ 0 (r, r ; ω) that grows as O(N 4 ) with system size, the total cost of calculating these residues grows thus as O(N 5 ). In addition, for a given system, as the energy window associated with the two Heaviside functions increases for energies located far away from the gap, the number of residues to be considered grows. Namely, the prefactor associated with the O(N 5 ) scaling increases as the states for which Σ C is calculated go deeper into the occupied or unoccupied manifold.

To bypass these difficulties, we explore an approach where the residues W (r, r ; ε n -E) are obtained by using an analytic continuation (AC) of the screened Coulomb potential W to the real-energy axis. This approach has several advantages: i) as argued and demonstrated in the previous section, the screened Coulomb potential is much less structured than the self-energy, considerably stabilising the AC procedure; ii) it can be done with no additional cost when continuation is performed from the imaginary axis only since the integration in Eq. 7 already requires the explicit calculation of W (iω) for a set of imaginary frequencies (Fig. 2a); iii) contrary to the case of Σ C (ω), the computational cost of W (ω) is constant all over the complex plane, so that extra reference points can be added away from the imaginary axis to increase the range of accuracy of the AC procedure without losing the global O(N 4 ) scaling of the whole calculation.

The merging of the contour deformation and analytic continuation of W following Eq. 7 is explained in the two following sections. We first introduce the specific quadrature along the imaginary axis for integration of the GW self energy. The advantage of performing the AC on the screened Coulomb potential to obtain efficiently and accurately the needed residues is then discussed.

Quadrature for the imaginary-axis integral

We start the presentation of our numerical treatment of Eq. 7 by describing the quadrature we adopt for the calculation of the imaginary-axis integral contribution, namely the integral part I(E) of φ m |Σ C (E)|φ m :

I(E) = dω φ m |G(r, r ; E + iω) W (r, r ; iω)|φ m (8)
where the integral ranges from -∞ to +∞. One can rewrite Eq. 8 using once more the formal functional forms of G and W , leading to:

I(E) = λ n I n λ (E) |w m,n λ | 2 (9)
with w m,n λ = φ m φ n |w λ and

I n λ (E) = dω 1 E + iω -ε n 1 iω -Ω λ - 1 iω + Ω λ (10)
The value of I n λ (E) depends on the sign of (Eε n ):

I n λ (E) =            -2iπ E-Ω λ -εn E < ε n -2iπ E+Ω λ -εn E > ε n (11)
Since the input energies ε n are known, contrary to the Ω λ poles of W (r, r ; ω), we can define an energy-specific quadrature Q(∆E) = ω k , z k (∆E) defined so as to minimise the error:

Ωmax Ω min dΩ k z k (∆E) iω k -Ω - z k (∆E) iω k + Ω - 2iπ ∆E + Ω (12)
Using such quadrature, the imaginary axis integral contribution I(E) can be estimated as

I(E) = n sgn(ε n -E)I n (E) (13)
with

I n (E) = k z k (|E -ε n |) φ m φ n | W (iω k )|φ m φ n (14)
Such a contribution only requires the knowledge of W (iω) along the imaginary axis where it is smooth. The quadrature bounds Ω min and Ω max are adjusted depending on the system, typically setting Ω min as a fraction (a half in this work) of the input (e.g. Kohn-Sham) gap

and Ω max as the maximum input (ε aε i ) transition energy. Let's emphasise here that the z k (∆E) weight are computed "on the fly" for each ε n contribution, while the w k sampling points are common to all these contributions and defined only once from the Ω min and Ω max parameters.

Anticipating on the results presented below, we collect in Table 1 the HOMO and LUMO energies obtained at the def2-QZVP G 0 W 0 @PBE level for the problematic MgO molecular system as a function of the number of {ω k } quadrature frequencies along the imaginary axis, demonstrating a high level of accuracy even with a limited set of quadrature points.

Analytical Continuation of W

Although there are numerous techniques to do so (Pade approximants, AAA, [START_REF] Nakatsukasa | The AAA Algorithm for Rational Approximation[END_REF] RKFit, [START_REF] Berljafa | The RKFIT Algorithm for Nonlinear Rational Approximation[END_REF] etc.) we choose here to perform the analytical continuation of the needed W (ω) residues:

W nm (ω) = φ m φ n | W (ε n -ω)|φ n φ m ( 15
)
using the continued fraction approach proposed in Ref. 62. Namely, given a set of reference points and associated values {(ω k , f k )}, the domain of the corresponding function is extended to the complex plane by means of a continued fraction:

f (z) := a1 1+ a2(z -ω 1 ) 1+ a3(z -ω 2 ) 1+ . . .
This functional form is constructed recursively so that at stage p, the coefficient a p is set for f (ω p ) to take exactly the value f p . Because of the (zω k ) factors, adding a reference point p doesn't affect the values of f (ω k ) for k < p. Such a form benefits in particular from recursion relations that allow for fast calculations of {a i } coefficients and fast evaluations at any complex z (see Ref. 62 Appendix for details). Anticipating on the two AC schemes presented below, we emphasise that there is no restriction on the location of the ω k reference points that can be located anywhere in the complex plane.

We slightly adapted to the present problem the algorithm of Ref. 62 in the following way:

i) since W is symmetric in frequency, namely W (ω) = W (-ω), the corresponding functional form should be:

s(z) := a1 1+ a2(z 2 -ω 2 1 ) 1+ a3(z 2 -ω 2 2 ) 1+ . . . , (16)
the recurrence relations of Ref. 62 being modified accordingly; ii) the inclusion order of the reference points in the continued fraction is determined so as to minimise the mean square error between the remaining references and continued expression at each recursion step; iii) a reference point is included in the continued fraction only if the AC result and the reference value differ by a significant amount (> 10 -8 × f ω=0 here). In practice these last two considerations greatly stabilise the AC procedure and allowed us to work with large set of reference points without experimenting trouble with numerical stability.

In this work, we thus propose to rely on the AC for the estimation of the W nm (Eε n ) residues arising during the calculation of the φ m |Σ C (E)|φ m correlation energy. For this, we experimented two different sampling strategies as exemplified in Figs. 2a-b, allowing to control the performance/accuracy ratio of the calculation.

The strategy of Fig. 2a considers reference points distributed along the imaginary axis only, taken identical to the points that were used to perform the numerical integration in Eq. 7. Once the screened Coulomb potential W is constructed for a selected {iω k } grid along the imaginary axis, the construction of the W nm values and their analytic continuations over the complex plane are obtained in negligible time. This matches the "traditional" analytic continuation approach relying on data points calculated only at imaginary frequencies. As shown below, very few sampling frequencies are necessary for an accurate analytic continuation of W in the vicinity of the imaginary axis, allowing to get accurate quasiparticle energies for states relatively close to the energy gap. Following the convergence data reported in Table 1 for the integral contribution to the self-energy, we will adopt in the following a limited nω=14 number of points along the imaginary axis. self-energies Σ(E) for E within ω max of the frontier orbitals. For sake of illustration, if one desires the quasiparticle correction to all valence states, with a valence bandwidth of about 30 eVs, (ω max /∆ω) amounts to 30 additional W matrices to be calculated. Due to the parity properties of W , the same calculated data points can be used to obtain the quasiparticle correction to all empty states within 30 eV from the gap. The same scheme will be adapted to core states in Section 6.

The ∆ω spacing can be also reduced to much smaller values to check the convergence and acquire reference quasiparticle energies. We verified that a spacing of ∆ω=125 meV leads to a converged analytic continuation scheme. This will be our strategy in the following to validate our calculations. In the following of the paper, such calculations will be referred to as our "reference" calculations with in parenthesis (Fig. 2b scheme, ∆ω=125 meV) as a reminder.

A parameter to be considered is the height η of the additional grid points with respect to the real axis. Experimentation leads us to consider that η = 1.5 × ∆ω, with ∆ω the grid spacing, was a reasonable value. We observed in particular that grid points lying too close to the real axis are detrimental to the AC precision, being probably too sensitive to the local pole structure of W (ω). On the other hand, having the reference point too distant from the AC region would decrease the AC accuracy.

To conclude this section on the analytic continuation of W , we note that the scheme of Fig. 2b does not allow to get an accurate AC expression for W (ω) everywhere in the complex plane, but only in the "vicinity" of the reference sampling points over which W (ω) is explicitly calculated. One therefore cannot assume that the resulting AC of W (ω) will be accurate at much higher energy along the real-axis. As such, we did not attempt to use this AC form to perform explicitly the integral of the G(E + ω)W (ω) integrand up to ω = +∞ along the real-axis following Eq. 1. In the contour deformation approach, W is only required over a limited energy range, namely close to the (ε n -E) poles of the G(E + ω) Green's function that entered the first/third quadrants.

Technical details

Benchmark calculations are first performed on the GW 100 molecular test set. [START_REF] Wilhelm | Toward GW Calculations on Thousands of Atoms[END_REF][START_REF] Van Setten | Benchmarking G0W0 for Molecular Systems[END_REF][START_REF] Krause | Coupled-cluster reference values for the GW27 and GW100 test sets for the assessment of GW methods[END_REF][START_REF] Caruso | Benchmark of GW Approaches for the GW100 Test Set[END_REF][START_REF] Maggio | A Plane Wave Perspective for Small Molecules[END_REF][START_REF] Govoni | Comparison of Methods and Accuracy of Results Obtained with the WEST Code[END_REF] Molecular structures are taken directly from the GW 100 original paper. [START_REF] Van Setten | Benchmarking G0W0 for Molecular Systems[END_REF] Out of the full GW 100 test set, we consider only the 93 systems not requiring pseudopotentials. For vynil-bromide and phenol, the corrected structures, [START_REF] Maggio | A Plane Wave Perspective for Small Molecules[END_REF] and corresponding modified data, are taken from the GW 100 website. 63 Such a test set encompasses rare gaz atoms, small diatomic molecules, etc. up to larger systems such as the DNA nucleobases. Our calculations are performed at the def2-QZVP level [START_REF] Weigend | Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy[END_REF] starting with input PBE [START_REF] Perdew | Generalized Gradient Approximation Made Simple[END_REF] Kohn-Sham eigenstates to allow direct comparison with the available TURBOMOLE [START_REF] Van Setten | The GW-Method for Quantum Chemistry Applications: Theory and Implementation[END_REF] (TM) and FHI-aims 66 localized basis set calculations from Ref. 49. The input Kohn-Sham calculations are obtained with the NWChem package. [START_REF] Valiev | NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations[END_REF] Coulomb-fitting resolution-of-identity (RI) [START_REF] Whitten | Coulombic potential energy integrals and approximations[END_REF] techniques are employed.

We adopt the auxiliary def2-QZVP-RIFIT basis sets [START_REF] Pritchard | New Basis Set Exchange: An Open, Up-to-Date Resource for the Molecular Sciences Community[END_REF] that we consider to be better converged than the compact and universal Weigend Coulomb fitting ones [START_REF] Weigend | Accurate Coulomb-fitting basis sets for H to Rn[END_REF] The second molecular test set we explore is composed of 24 organic acceptors as proposed by Marom and coworkers. [START_REF] Knight | Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods[END_REF] Such molecules have significant importance for e.g. hole-doping purposes and their large electronic affinity (EA) allows reliable measurements. For comparison with available CCSD(T) data, [START_REF] Richard | Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules I. Reference Data at the CCSD(T) Complete Basis Set Limit[END_REF] we adopt here the augmented aug-cc-pVTZ basis sets together with their corresponding auxiliary basis sets aug-cc-pVTZ-RI. [START_REF] Pritchard | New Basis Set Exchange: An Open, Up-to-Date Resource for the Molecular Sciences Community[END_REF] We compare in particular the ionization potential (IP) and electronic affinities (AE) obtained at the nonself-consistent G 0 W 0 @PBE0 level and "gap-self-consistent" eg-GW @PBE0 approach with the available CCSD(T) data. Calculations are performed with the newly developed BeDeft program (BeyondDFT) that originates from rewriting and extending the GW and Bethe-Salpeter Fiesta code. [START_REF] Li | Combining the Many-Body GW Formalism with Classical Polarizable Models: Insights on the Electronic Structure of Molecular Solids[END_REF][START_REF] Jacquemin | Benchmarking the BetheSalpeter Formalism on a Standard Organic Molecular Set[END_REF][START_REF] Duchemin | The BetheSalpeter formalism with polarisable continuum embedding: reconciling linear-response and state-specific features[END_REF] 4 Results

The GW 100 test set

We now explore the merits of performing the analytic continuation on the screened Coulomb potentiel W for the residues in the contour-deformation scheme by studying the def2-QZVP G 0 W 0 @PBE quasiparticule HOMO and LUMO energies of the GW 100 test set. We thus calculate the quasiparticle energies obtained by the AC of W starting from a very limited (nω=14) set of W (ω) matrices calculated at imaginary frequencies (scheme Fig. 2a) that we compare with the accurate reference calculations (Fig. 2b scheme, ∆ω=125 meV). The data are compiled in the Supporting Information (Tables S1 andS2).

The largest deviation with respect to the reference calculations amounts to 6 meV for the He atom HOMO for which the quasiparticle correction to the input Kohn-Sham energy is the largest ( 8 eVs). For such states with very large GW correction, the (Eε KS n ) energies at which residues must be calculated, with E E GW m the targeted quasiparticle energy, are lying far away from the imaginary axis. Besides He, 6 molecules show a discrepancy of the order of 1 to 2 meVs, all other systems displaying a deviation below the meV. In the case of the AEs, one system (He again) presents a 1 meV deviation. For sake of completeness, we also provide in the SI the def2-QZVP G 0 W 0 @PBE0 data that display even smaller differences between the two schemes due to reduced GW corrections. With a PBE0 starting point, the He HOMO displays the only discrepancy larger than a meV in relation again with a GW correction that is still substantial ( 6 eV).

To further analyse the remarkable success of the present analytic continuation approach based on very few calculated data points along the imaginary axis, we now focus on the BN, O 3 , BeO and MgO systems that were leading to difficulties when using the analytic continuation directly on the final self-energy matrix elements. Table 2 reports the corresponding ionization potential (IP) taken as the negative of the HOMO def2-QZVP G 0 W 0 @PBE value.

Our data are compared to the published values obtained with the Turbomole and FHI-Aims [START_REF] Van Setten | Benchmarking G0W0 for Molecular Systems[END_REF] We report the TM calculations performed with full diagonalization of the (RPA) dielectric matrices and no resolution-of-the-identity approximation, namely exact calculations that we label TM-noRI. The AIMS-16 and AIMS-P128 data were based on the analytic continuation of the self-energy using N=16 and N=128 parameters in the Padé expansion of the self-energy matrix elements. Our results are given in the three last lines, using the notation beDeft-ref for the reference calculations (Fig. 2b scheme, ∆ω=125 meV) and beDeft-P14 for the efficient (Fig. 2a scheme, nω=14) approach. The renormalization Z factor is defined as 1/Z = 1 -∂Σ(ω)/∂ω where the derivative is taken at the quasiparticle energy.

For each system and approach, we provide when available two values associated with the largest Z factors. The existence of several possible solutions of the quasiparticle equation:

ε QP n = ε KS n + φ n |Σ GW (ε QP n ) -V XC |φ n ( 17
)
is illustrated in Figures 3 and4 (lower panels) for the specific MgO and BN cases. The straight red lines represent the To better understand the improved stability of the present scheme, we now analyse in

f (ω) = ω -ε KS H -φ H |Σ X -V XC
(ω) = ω -ε H -φ H |Σ X -V XC |φ H function.
The position of the G 0 W 0 @P BE HOMO quasiparticle energy is indicated by the vertical black dashed line.

eV originates from the (n=H-2) W nH (ω) contribution characterized by a single pole. As a result, the analytic continuation using only data points along the imaginary axis (Fig. 2a) does an excellent job with our limited (nω=14) imaginary frequencies grid, reproducing very faithfully the main pic close to ε QP H quasiparticle energy and leading to an accurate determination of the quasiparticle energy. The same quality of fit is obtained for the (n=H-1) contribution capturing accurately the pole at -8.6 eV. As a matter of fact, the number of reference points along the imaginary axis over which χ 0 (iω) is explicitly calculated can be reduced down to about nω=8 to retain the meV accuracy on the HOMO energy as shown previously in Table 1.

The ideal case of one-pole only per W nH (ω) matrix element cannot be a general situation as shown in the (n=H) case where 3 poles are observed in the reference calculation. In this specific case, the analytic continuation of W nH (iω) only captures the two broader pics, missing the sharpest one that is associated with a very small w mn λ weight. Clearly, the very same pole is captured in the case of (n=H-1) where it appears with a larger weight. As such, the resulting pole in the self-energy is well reproduced with only a small deficit of weight very close to the corresponding energy resonance. We note that the pole at about -3.4 eV reported in Fig. 13 We now turn in Fig. 4 to the case of the BN dimer for which a 128-Padé fit was required when attempting to perform a direct analytic continuation of the self-energy. The situation is similar to that of the MgO system with several poles dominating the φ H |Σ C (E)|φ H correlation energy dependence. However, an additional complication arises since the quasiparticle energy is related to the second pole of φ H |Σ C (E)|φ H (going away from the gap). Again, the poles of the self-energy are distributed amongst various W nH matrix elements. Remarkably, the analytic continuation with a limited nω=14 imaginary axis sampling reproduces very accurately the sharp structure at -9.9 eV. Such a very good fit is made possible by the fact that it is the only pole of the (n=H-2) contribution. One pole at -15.6 eV is missed in the (n=H) contribution but it is captured by the (n=H-1) contribution, appearing with good accuracy in the overall self-energy dispersion.

We summarize this Section by concluding that distributing the poles of the self-energy amongst several W nm (ω) functions that are individually continued significantly stabilizes the analytic continuation approach, resolving in particular all the difficulties encountered for the GW 100 test set with a remarkably compact nω=14 energy-grid along the imaginary axis. We note that this grid was optimized for the integral contribution along the imaginary axis, and not for the analytic continuation. Additional points along the imaginary axis may potentially lead to a better capture of poles in the contributing W nm (ω) residues but without improving our results for the GW 100 test set.

Acceptor molecules with large EA

We now turn to the set of Ref. 58 that contains 24 moderate size organic molecules relevant to organic electronics and light emitting devices. We compile in Table 3 our G 0 W 0 @PBE0 values for the IP and AE energies. Our calculations are performed at the aug-cc-pVTZ level for which selected reference CCSD(T) data 72 are available. As now well documented, non-self-consistent G 0 W 0 calculations starting with DFT Kohn-Sham eigenstates generated with limited amount of exact exchange (25% in the PBE0 case) lead to a significant underestimation of IP and AE values. The mean signed error (MSE) for the IP and AE amount to -0.32 eV and 0.32 eV, respectively, with maximum errors of -0.41 eV (dichlone IP) and 0.44 eV (TCNE AE).

As compared to our reference G 0 W 0 @PBE0 calculations (Fig. 2b scheme, ∆ω=125 meV), the results obtained with the AC from nω=14 data points along the imaginary-axis yield a maximum error of the order of 10 -2 meV over the entire set, demonstrating again the robustness of the present AC scheme. This error is even smaller than that obtained with the GW 100 test set, presumably because the GW corrections are much smaller in the present case. This results from the PBE0 starting point (instead of PBE) and the larger size of the systems considered. Such considerations indicate that the energies (Eε KS n ) at which the residues must be calculated are smaller, namely closer to the imaginary axis.

While various schemes were explored to improve G 0 W 0 @PBE0 data for this set of molecules, [START_REF] Knight | Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods[END_REF] including optimal DFT starting points generated with tuned range-separated hybrids, 76 selfconsistency on both eigenvalues and eigenstates (scGW ), or by including more diagrams such as the second-order screened exchange (SOSEX) term, we explore here a very simple approach, namely the partially self-consistent egGW scheme where only the value of the HOMO-LUMO gap is updated (-eg-stands for energy gap) in the spirit of a self-consistent scissor approximation. This is the simplest self-consistent scheme for approaches providing at reduced cost the energy of frontier orbitals. The egGW calculations reduce the mean absolute error (MAE) on the IP from 0.329 eV to 0.135 eV, as compared to the single-shot G 0 W 0 @PBE0 technique, and from 0.307 eV to 0.098 eV in the case of the AE.

Beyond frontier orbitals

The egGW scheme is a simplification of the popular evGW approach 52,77-79 where all eigenvalues are self-consistently updated. While the egGW approach only requires calculating the HOMO and LUMO quasiparticle energies, the full evGW approach requires calculating the energy of states located far away from the gap. Beyond the evGW approach, accessing reliably specific levels located far away from frontier orbitals is an important issue. One could invoke e.g. core states [START_REF] Golze | Core-Level Binding Energies from GW: An Efficient Full-Frequency Approach within a Localized Basis[END_REF][START_REF] Van Setten | Assessing GW Approaches for Predicting Core Level Binding Energies[END_REF] or the frontier orbitals of molecular systems deposited on metallic electrodes that are separated from the Fermi level by the "continuum" of metallic [START_REF] Neaton | Renormalization of Molecular Electronic Levels at Metal-Molecule Interfaces[END_REF][START_REF] Thygesen | Renormalization of Molecular Quasiparticle Levels at Metal-3 Acceptor molecular set (Ref. S2 ) ionization potential and electronic affinity Table S3: Ionization potentials (IP) and electronic affinities (EA) for the 24 molecules of Refs. S2,S3[END_REF] This is a severe test for analytic continuation techniques relying on data acquired only along the imaginary axis as we now discuss.

To illustrate the difficulty associated with deeper states, we extract from the previous set of moderate size molecules the paradigmatic TCNQ electron acceptor that belongs to an important family (F4TCNQ, F6TCNNQ, etc.) of p-type dopants in organic semiconductors. We select for illustration the n=HOMO-35th state located 18.5 eV below the gap (G 0 W 0 @PB0 value). We start again with reference calculations by using a very large set of energy grid points close to the real axis (Fig. 2b scheme, ∆ω=125 meV) on which W is explicitly calculated for the construction of an accurate AC to the complex plane.

Figure 5: Spectral function A GW n (ω) associated with the TCNQ n=HOMO-35 state located at about 18.5 eV below the gap (G 0 W 0 @P BE0 value). The spectral function is calculated along paths parallel to the real axis with a constant (ω) imaginary part of 10, 50, 200 and 400 meV. The blue dashed lines are the one-pole A f it n (ω) Lorentzian fit to the spectral function used to extract the quasiparticle energy E QP , lifetime Γ QP and renormalization factor Z QP .

We plot in Fig. 5 the corresponding spectral function :

A GW n (ω) = 1 π φ n |G(ω)|φ n with G = G 0 + G 0 (Σ GW -V XC )
G and G 0 the input PBE0 Green's function. More specifically, we represent A GW n (ω) along ω-lines parallel to the real-axis with a constant imaginary part (ω) ranging from 10 meV to 400 meV. For very small (ω) value ( 10meV), namely close to the real axis, the spectral function exhibits a very dense "forest" of peaks corresponding to a large density of solutions of the quasiparticle equation, all solutions associated with a very small Z factor. The presence of this very large density of poles for the self-energy in this energy range originates from the large number of poles of the susceptibility associated with transitions from occupied states to the quasi-continuum of empty states above the vacuum level.

Any analytic continuation from the imaginary axis will experience great difficulties to describe such a complex structure along (or very close) to the real axis. We emphasize that the poles of W in this energy range becomes denser and denser as the size of the Kohn-Sham basis set increases, following the densification of empty states above the vacuum level.

Eventually, all sharp contributions will merge into a broad structure forming the quasiparticle peak with a finite spectral weight.

This formation of a smooth quasiparticle peak can be achieved by increasing the (ω) value, as shown in Fig. 5. When (ω) becomes sufficiently large, the quasiparticle peak clearly appears together with the incoherent background at lower energy. To extract the quasiparticle energy, we fit the spectral function by a one-pole spectral (or Lorentzian) function:

A f it n (ω) = 1 π Z QP ω -(E QP + iΓ QP ) = Z QP π (ω) -Γ QP ( (ω) -E QP ) 2 + ( (ω) -Γ QP ) 2 (18)
resulting in the dashed-blue lines of Fig. 5. Together with the standard quasiparticle energies E QP and renormalization factor Z QP , one obtains the inverse lifetime Γ QP . The lorentzian fit leads to a stable quasiparticle energy nearly independent of the (ω) value. The renormal-ization factor Z QP and inverse lifetime Γ QP are also found to be reasonably stable provided that (ω) is chosen not too small.

In the case of a well-defined quasiparticle peak (i.e. a single, well separated pole), such as the frontier orbitals discussed here above, this method yields the very same result as the one obtained by directly solving the quasiparticle equation (e.g. Fig. 3). This is illustrated in Fig 6 in the case of the TCNQ HOMO level where the E QP , Z QP and Γ QP are completely insensitive to the choice of (ω). In such cases, the fit yields a vanishingly small Γ QP inverse lifetime, as expected for states close to the gap. The present analysis suggests that studying the spectral function A GW n (ω) along energy lines slightly shifted away from the real-energy axis, rather than attempting to capture all details of the self-energy along the real-axis, provides a much simpler and stable way to extract quasiparticle energies, together with meaningful renormalization factor Z QP and inverse lifetime Γ QP . As exemplified here above, the obtained quasiparticle characteristics extracted through the one-pole fit strategy (Eq. 18) are very much insensitive to the choice of the (ω) value. In the following, we adopt this strategy with (ω) = 100 meV to evaluate the merits of the present analytic continuation approach.

We thus perform the same analysis for all molecules of the set, considering all states within 20 eV from the gap (PBE0 value), encompassing namely over 600 occupied and 2600 unoccupied levels. As above, we first establish reference (E QP , Z QP , Γ QP ) values using explicitly calculated W matrices over a dense grid of energy points close to the real-axis (Fig. 2b scheme, ∆ω=125 meV). From this large data set, W is analytically continuated to the (ω)= 100meV axis parallel to the real-axis along which the self-energy and spectral function A QP n (ω) are constructed. The Lorentzian fit of this accurate spectral function leads to reference quasiparticle data.

We now perform analytic continuations of W to the (ω)=100meV axis starting from much reduced sets of explicitly calculated W matrices. We start by the set of W matrices calculated over the nω=14 integration points along the imaginary-axis as in Sections 4.1 and 4.2. After Lorentzian fit of the resulting spectral function A QP n (ω), the extracted quasiparticle energy E QP n is compared to the reference values. We plot in Fig. 7(a) the resulting error in meV (log scale).

For frontier orbitals, namely the HOMO and LUMO levels, the errors are vanishingly small (≤ 0.01 meV). This is consistent with the results of Sections 4.1 and 4.2 revealing the accuracy of the analytic continuation of W using few data points obtained along the imaginary axis only. Further, this confirms that the extraction of the quasiparticle energies through the one-pole fit of the spectral function is stable and accurate for frontier orbitals.

Very consistently, recalculating the HOMO and LUMO energies of the GW 100 test set with the present scheme produces trifling errors as compared to the results of Section 4.1 where quasiparticle energies were obtained by solving the quasiparticle equation along the real axis. 58). In (a), the analytic continuation of W relies on the nω=14 data points calculated along the imaginary axis (see Fig. 2a). In (b) additional reference W matrices are calculated over a very coarse energy-grid parallel to the real-axis with a ∆ω =1 eV spacing (see Fig. 2b), resulting in a more accurate analytic continuation of W to the (ω) = 100 meV energy axis.

The red cross with a circle identifies the main outlier over the entire set (see text).

Similarly, we verify that the quasiparticle spectral weight Z QP stemming from the one-pole fit is identical to the renormalization factor obtained previously from the energy derivative of the self-energy at the quasiparticle energy.

Apart from frontier orbitals, the errors remain acceptable (below 10 meV) for states within one energy gap of frontier orbitals. In this energy range, the quasiparticle renormalization factor Z QP remains significant (≥ 0.65) indicating well defined quasiparticles. As such, the present analytic continuation scheme, relying on calculating W over a few points along the imaginary axis, is accurate not only for the HOMO and LUMO levels, but also for states located typically within one energy gap from frontier orbitals. As expected however, for states deeper in energy, the error increases, becoming occasionally close to an eV.

We now increase the accuracy for deeper state by explicitly calculating W over a very coarse grid of frequencies following the scheme of Fig. 2b but with a large ∆ω =1 eV energy spacing between grid points, leading to calculating W matrices over nω 28 additional points on top of the nω=14 points along the imaginary axis. The present scheme is thus intermediate between the "reference" calculations (∆ω =125 meV) and the analytic continuation of W from the imaginary-axis-only data points. Contrary to the self-energy operator Σ(ω), calculating the screened Coulomb potential W (ω) at frequencies with a large real-part does not present additional difficulties as compared to calculating W (ω) on the imaginary axis, offering in particular the very same cost and O(N 4 ) scaling with system size. 83 The error associated with this intermediate scheme is reported in Fig. 7(b). Clearly, adding 28 frequencies in the first-quadrant over which W is explicitly calculated allows to significantly reduce the error within the ±20 eV s range explored. In particular, the error for states located within an energy gap of frontier orbitals decreases dramatically to negligible values. Further, the errors fall systematically below 0.1 eV within the ±3 energy gap window.

Out of the 600 occupied states considered, one clear outlier presents a large discrepancy between the quasiparticle energy associated with the ∆ω = 1 eV coarse grid and the reference calculation (see red circled cross in Fig. 7). We plot in Fig. 8 the corresponding spectral function associated with the HOMO-15 state of nitrobenzonitrile located ∼ 6eV below the gap (G 0 W 0 value). We take this outlier as a good example of the limit of a one-pole fit procedure. As can be observed in Fig. 8 where the reference spectral function are provided, A GW n (ω) is dominated by two peaks of similar weight, each peak being characterized by a small Z-factor of about 0.3.

In both the reference (Fig. 8c) and imaginary axis only (Fig. 8a) calculations, the onepole fit selects the peak with the largest intensity. However, the AC from the imaginary axis reproduces both peaks with significant errors regarding their energy location and intensity. As a result, the lower energy structure is erroneously chosen as the quasiparticle peak. In the case of the analytic continuation based on calculating W over the coarse ∆ω=1 eV energy grid (Fig. 8b), although the structures of A GW n (ω) are well-captured (see dashed vertical lines), the fit procedure fails to select a single peak and averages both structures with a very broad lorentzian. The resulting Z factor exceeds unity, providing a good indication that the fit failed to identify the quasiparticle correctly.

We conclude from the present example that for deeper states, the AC technique can be made accurate with additional W reference points in the first-quadrant, even in situations where the spectral function is dominated by several structures. However, there may be difficulties in attempting to extract a well defined quasiparticle peak from the spectral function. Such an observation questions the concept of quasiparticle rather than the quality of the AC scheme. The error analysis for all occupied/unoccupied states is provided in the SI, showing no other outliers below the gap, and only 3 states out of 2600 with an error larger than 0.1 eV in the unoccupied manifold at larger energy ( 4 gaps value).

We finally plot in Fig. 9 the Γ QP (signed) inverse lifetime. As expected, we observe a significant increase as a function of the distance in energy from mid-gap, even though not following clearly the standard quadratic (E -E F ) 2 dependence associated with 3D

Fermi liquids. 84 Discussing whether the concept of quasiparticles remains valid for states located deep into the occupied/unoccupied manifolds is beyond the scope the present paper concerned with the analytic continuation approach. We can only recommend for very low lying states to analyze the full spectral function that can be accurately captured by the analytic continuation of W following the scheme of Fig. 2b, even in the limit of a very coarse ∆ω sampling. Too small Z factors may be chosen as an indication that difficulties potentially exist.

Core states

We conclude this exploration by considering core states, focusing on the H 2 O 1s energy level as a test case. For such core levels, the residuals W (ε i -E C ), where E C is some typical core energy, can be partitioned into two categories, one with ε i a core level leading With such reference points, we analytically continue W to the (ω)=100 meV line along which we plot the real-part of the correlation self-energy (Fig. 11a) and the spectral function (Fig. 11b). We chose the energy range considered in the Fig. 5 (full black lines) reference calculation obtained using Fig. 2b scheme with ∆ω = 125 meV ;

(dashed blue lines) calculation obtained using Fig. 10 scheme with ∆ω = 1 eV grid spacing;

(dashed pink lines) calculation obtained using Fig. 2a scheme with nω=14 points along the imaginary axis only. The correlation self-energies and spectral functions are plotted for (ω) = 100 meV . The Inset zooms the spectral functions around the quasiparticle energy.

As expected, the analytic continuation based on sampling data points along the imaginaryaxis only (pink dotted-dashed lines) fails in reproducing faithfully this complex self-energy structure. Still, we observe that with our limited nω=14 sampling of the imaginary axis, the resulting AC on W leads to a few structures in the vicinity of the QP peak. This contrasts with the results obtained with the AC of Σ C (Ref. 57) that does not capture any pole in this energy range, even with a very large (1024 pts) sampling set along the imaginary axis.

We did not attempt here to optimize the sampling strategy, by e.g. differentiating sampling grids and ∆ω spacing associated with valence and core residues. However, Fig. 11 demonstrates that Fig. 2b sampling strategy with a coarse ∆ω sampling of 1 eV already provides a well converged spectral function at moderate cost.

Conclusion

The contour-deformation approach to calculating the GW quasiparticle energies requires the knowledge of the screened Coulomb potential W along the imaginary frequency axis but also at specific energies on the real-frequency axis. For these latter contributions, we have shown that the analytic continuation from the imaginary to the real-frequency axis of the screened Coulomb potential W (ω) leads to a much more robust scheme as compared to current techniques involving the direct continuation of the self-energy Σ(E). As a result, the present approach allows to reach the meV accuracy for the calculation of the HOMO and LUMO levels energy of the GW 100 molecular set, [START_REF] Van Setten | Benchmarking G0W0 for Molecular Systems[END_REF] calculating only a very few (nω 14)

W (iω) matrices along the imaginary axis. In particular, all difficulties previously observed for selected systems (MgO, BN, O 3 , BeO) are now resolved. Even better results are obtained with a recent set of larger medium-size real-life acceptor molecules. [START_REF] Knight | Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods[END_REF] In this latter case, the error associated with the analytic continuation from the imaginary axis is reduced to about 10 -2 meV for frontier orbitals when starting with PBE0 input eigenstates.

The calculation of the GW correction to states located far away from the gap remains a challenging issue. The difficulties originate from the expected behaviour of such quasiparticle states with the decay of the main quasiparticle peak over a very dense forest of poles along the real axis, yielding a very large number of solutions of the quasiparticle equation associated with very small Z QP renormalization factors. The number and location of such peaks are unphysical, being related in particular to the finiteness of the Kohn-Sham basis set adopted to generate the input one-body eigenstates. However, as expected, these peaks merge rapidly into broad quasiparticle structures if one calculates the spectral function A GW n (ω) along ωlines parallel to the real-axis but with a finite (ω) imaginary part.

For such broad structures, the analytic continuation of W from the imaginary-axis to the (ω) 100 meV axis yields errors below 10 meV for states within typically an energy gap of frontier orbitals, but the error can become as large as an eV further away in energy. In such situations, the analytic continuation scheme can be consolidated by calculating a few screened Coulomb potential matrices in the complex plane, increasing the accuracy of the analytic continuation in the vicinity of each added reference point. In practice, adding a very coarse grid (∆ω 1 eV) of reference W (ω) calculated along an ω-line parallel to the real axis, allows to bring the error of the present analytic continuation approach below 10 meV over a much larger energy range. Contrary to the self-energy operator, we emphasize that calculating W (ω) at frequencies anywhere in the complex plane bears the same computational cost.

The same strategy can be used for core states as demonstrated in the paradigmatic case of the H 2 O 1s level. Sampling separately the energy ranges spanned by core and valence poles allows for an efficient yet accurate calculations of the core state spectral function over a large energy window.

Difficulties remain occasionally with states far away from the gap for which there are no clear quasiparticle peaks with sizeable Z QP weight dominating the spectral function. In such situations, the analytic continuation of W using a coarse grid of data points in the complex plane remains accurate, but instabilities can occur when attempting to extract automatically one peak out of several equivalent structures. As such, the analytic continuation scheme is not to be blamed, but rather the concept of quasiparticles. Such cases question the use of self-consistent approaches where the quasiparticle energy of low-lying states is re-injected to build an updated Green's function G and screened potential W . These difficulties are well beyond the scope of the present paper.

Molecule Interfaces: Trends across Binding Regimes. Phys. Rev. Lett. 2009, 102, 046802.

(83) The standard technique relying on performing the analytic continuation of the selfenergy may also be improved by calculating reference Σ(ω) self-energies away from the imaginary-axis. We re-emphasize that calculating W always scales as O(N 4 ) for any frequency in the complex plane, while calculating the self-energy operator scales as O(N 5 ) for states away from the gap using the standard contour-deformation approach, and O(N

Analytic continuation error

Figure S1: Same as Fig. 7b of the main manuscript, but over the full energy range and with a lower Z = 0.2 cutoff value. For occupied states, no new outliers are observed. In the unoccupied manifold, the main outlier corresponds to the LUMO+139 of acridine located 16eV above the LUMO G 0 W 0 @PBE0 value.

Figure 1 :

 1 Figure 1: Contour-deformation technique. The poles of the G(E + ω) time-ordered Green's function for occupied/unoccupied states are indicated in blue/red. When the energy E is smaller than the HOMO energy, which happens when we seek the quasiparticle energy for occupied states, poles of G(E + ω) start entering the first quadrant (case of the figure).

Figure 2 :

 2 Figure 2: Schematic representation of two possible frequency samplings (dots) over which the screened Coulomb potential W is calculated explicitly for the construction of an Analytic Continuation of W (ω) in the complex plane: a) along the imaginary axis only and b) with additional points parallel to the real axis.

  used in conjunction with the def2-QZVP Kohn-Sham basis in the TM-RI calculations of Ref. 49. Following a recent exploration of resolution-of-identity techniques, auxiliary basis sets are considered in their Cartesian representation, leading to more accurate values as compared to the same basis set used in a spherical representation. 71 Core levels are included in the calculation of the susceptibility.

  |φ H functions where (ε KS H , φ H ) are the PBE Kohn-Sham HOMO eigenstates. The possible quasiparticle energies are found at the crossing of such lines with the φ H |Σ C (ω)|φ H energy-dependent correlation self-energy expectation value (full black and dashed blue lines). The selected values with the largest Z factors correspond to the two solutions yielding the smallest IP. Such solutions are labeled the -left-and -right-solutions in Ref. 49. 75The analysis of the data in Table2 clearlyindicates that our results (beDeft-P14) are within 20 meV of the reference TM-noRI HOMO values. For the present systems, the error associated with the (P14) analytic continuation is smaller than 2 meV (BeO) for the HOMO quasiparticle energies. As a result most of the discrepancy with TM-noRI data originates from the resolution-of-the-identity approximation. In comparison, the AIMS-P16 approach could not identify all possible solutions, yielding further errors of the order of 120 meV, 560 meV and 320 meV for BN, BeO and MgO, respectively. Increasing the number of Padé parameters up to 128 (AIMS-P128), requiring thus the calculation of the self-energy matrix elements at 128 imaginary frequencies, could restore excellent results. However, the AIMS-P128 calculations does not identify the correct solutions for the O 3 and MgO systems, namely does not lead presumably to an accurate enough determination of the Z-factors. Such an analysis emphasizes the difficulties associated with techniques attempting to perform an analytic continuation of the self-energy matrix elements, as compared to the present scheme where we perform an analytic continuation on individual screened Coulomb potential matrix elements.

Fig 3 -

 3 Fig 3-(Upper panels) the individual W nH (ω) matrix elements (see Eq.15) with n = (H, H -1, H -2) indexing the MgO HOMO, HOMO-1, and HOMO-2 eigenstates. Such states are selected by the two Heaviside functions in the contour-deformation formula. The reference W nH (ω) matrix elements (Fig. 2b scheme, ∆ω=125 meV) are plotted with black solid lines. As discussed above, each individual W nH (ω) matrix element presents much less structures than the corresponding φ H |Σ C (ω)|φ H correlation self-energy term represented in the lower panel. In particular, the pole close to the targeted MgO HOMO quasiparticle energy at -6.67

Figure 3 :

 3 Figure 3: (Upper panels) Energy dependence of the W nH = φ H φ n |W (r, r ; ωε n )|φ n φ H residues for MgO. The wavefunction φ H is the PBE HOMO eigenstate while the {ε n , φ n } with n = (H -1, H -2) are the corresponding HOMO-1 and HOMO-2 eigenstates. These contributions are compared to the overall φ H |Σ C (ω)|φ H correlation matrix element (Lower panel). The AC data obtained with nω=14 imaginary frequency data points (dashed blue lines) are compared to reference calculations (full black lines) obtained with the fine (∆ω=125 meV) sampling of Fig. 2b. The red line represents the f(ω) = ωε Hφ H |Σ X -V XC |φ H function.The position of the G 0 W 0 @P BE HOMO quasiparticle energy is indicated by the vertical black dashed line.

  of Ref. 49 originates from the n=LUMO W nH contribution.

Figure 4 :

 4 Figure 4: As in Fig. 3 but for the BN dimer.

Figure 6 :

 6 Figure 6: Same as Fig. 5 but for the TCNQ HOMO spectral function.

Figure 7 :

 7 Figure 7: Plot of the error induced by the analytic continuation as a function of the distance in energy of the considered states with respect to the gap edges (PBE0 input values) in units of the corresponding PBE0 gap. The color code indicates the Z-factor. Errors are given in meV in log scale. The graph compiles the data for all the molecules of the acceptors set (Ref.58). In (a), the analytic continuation of W relies on the nω=14 data points calculated along the imaginary axis (see Fig.2a). In (b) additional reference W matrices are calculated over a very coarse energy-grid parallel to the real-axis with a ∆ω =1 eV spacing (see Fig.2b), resulting in a more accurate analytic continuation of W to the (ω) = 100 meV energy axis. The red cross with a circle identifies the main outlier over the entire set (see text).

Figure 8 :

 8 Figure 8: Spectral function and one-pole quasiparticle (QP) fit associated the HOMO-15 of Nitrobenzonitrile, for different AC strategies: a) reference W are calculated along the imaginary axis only; b) Fig. 2b scheme with ∆ω = 1 eV grid spacing along real axis; c) Fig. 2b scheme with ∆ω = 125 meV grid spacing along real axis (taken as reference calculation). The spectral function is plotted at (ω) = 100 meV . This illustrate a typical pathological case for the quasiparticle energy definition.

Figure 9 :

 9 Figure 9: Inverse quasiparticle lifetime Γ QP (in eV) as obtained from the one-pole quasiparticle fit of the spectral function.

  to small |ε i -E C | energies and the other with ε i a valence level leading to very large (several hundred eVs) |ε i -E C | energies. Following our sampling strategy of Fig. 2(b), we verified that calculating W (ω) matrices for Re(ω) located in the large energy range between |ε i -E C | for core and valence states does not influence the quality of the analytic continuation of W where it is needed. This leads to the sampling scheme of Fig. which is a variation on the sampling scheme of Fig. 2(b). In practice, with a ∆ω spacing of 1 eV, we calculate explicitly 44 W (ω) matrices for 0 ≤ (ω) ≤ 45 eVs to allow plotting Σ(E) on a large energy range around the 1s core quasiparticle energy, while 73 frequencies are kept for sampling the residues associated with valence states.

Figure 10 :

 10 Figure 10: 1s core state W sampling strategy for H 2 O. The ∆ω spacing above the real axis is 1 eV , leading to a total of 131 reference frequencies for the analytical continuation of W residues required by the calculation of Σ C (ω) within a large 50 eV energy window.

  of Ref.53 where the full frequency and standard contour-deformation approaches were compared in the study of H 2 O 1s core state. Our def2-QZVP G 0 W 0 @PBE quasiparticle energy at -523.80 eV (Z QP = 0.20) is identical to the -523.8 eV def2-QZVP G 0 W 0 @PBE value of Ref.53.

Figure 11 :

 11 Figure 11: GW calculation for the φ 1s H 2 O core state: a) real part of the correlation selfenergy φ 1s |Σ C (ω)|φ 1s expectation value; b) associated spectral function. We compare:(full black lines) reference calculation obtained using Fig.2bscheme with ∆ω = 125 meV ; (dashed blue lines) calculation obtained using Fig.10scheme with ∆ω = 1 eV grid spacing; (dashed pink lines) calculation obtained using Fig.2ascheme with nω=14 points along the imaginary axis only. The correlation self-energies and spectral functions are plotted for (ω) = 100 meV . The Inset zooms the spectral functions around the quasiparticle energy.

  ) are associated with the fine (∆ω = 125 meV ) Fig. 2b sampling scheme. The (P-14) calculations use only the nω=14 reference W (iω) along the imaginary axis from Fig. 2a sampling scheme. Calculations are performed at the def2-QZVP level with Coulomb fitting resolution-of-the-identity using the corresponding def2-QZVP-RI auxiliary basis. Both PBE and PBE0 Kohn-Sham DFT starting points are presented. Digits differing between the (ref.) and (P-14) data are highlighted in red.

Table 1 :

 1 Convergence of the solutions of the quasiparticle equation for the HOMO and LUMO levels in MgO (in eV, with Z factors in parenthesis) as a function of the number nω of {ω k } quadrature frequencies along the imaginary axis. The three solutions with largest Z are provided in the case of the HOMO level. The uncertainties stemming from the residue calculations using AC techniques are cast away by using a large and constant number of reference frequencies within the complex plane (see text Section 2.4).

	nω	E 1 H	E 2 H	E 3 H	E L
	6 -6.6538(0.29) -7.0835(0.19) -10.9733(0.10) -1.8752(0.80)
	8 -6.6696(0.28) -7.0943(0.20) -10.9787(0.10) -1.8907(0.80)
	10 -6.6711(0.28) -7.0953(0.20) -10.9793(0.10) -1.8922(0.80)
	12 -6.6712(0.28) -7.0954(0.20) -10.9793(0.10) -1.8923(0.80)
	14 -6.6714(0.28) -7.0956(0.20) -10.9794(0.10) -1.8924(0.80)

Table 2 :

 2 Ionization potential for BN, O 3 , BeO and MgO taken as the negative of the HOMO quasiparticle energy at the def2-QZVP G 0 W 0 @PBE level. The AIMPS-P16, AIMS-P128 and TM (no-RI) data are taken from Ref.49. The TM data should be considered as the reference def2-QZVP G 0 W 0 @PBE values. For each system and method, we provide when available the "left/right" solutions as defined in the Table4of Ref. 49 and in the main text. The bold values select the true quasiparticle energies associated with the largest Z value. The three last lines correspond to the present study, with the label "P14" for Fig.2ascheme and "ref" for Fig. 2b scheme with ∆ω=125 meV. Z values at second decimal place do not change between ref. and P-14 calculations. Energies are in eV. * We take the value from the GW 100 website that differs by 10 meV from the Table 4 value of Ref. 49.

		BN	O 3	BeO	MgO
	AIMS-P16	--/11.15	11.96/11.39	9.07/--	6.79/--
	AIMS-P128	11.67/11.03	11.96/11.39	9.63/8.58	7.11/6.68
	TM-noRI	11.67/11.01 *	11.95/11.39	9.63/8.62	7.09/6.66
	beDeft-ref	11.667/11.010	11.967/11.391	9.633/8.608	7.096/6.671
	beDeft-P14	11.671/11.010	11.967/11.392	9.635/8.608	7.096/6.671
	Z	0.28/0.48	0.40/0.25	0.47/0.15	0.20/0.28
	codes.

Table 3 :

 3 Ionization potential (IP) and electronic affinities (AE) aug-cc-pVTZ G 0 W 0 @PBE0 data for Ref. 58 test set. The CCSD(T) data are from Ref.72. The IP and AE yielding the largest discrepancy (Max.Err.) with the CCSD(T) data are highlighted. The Maximum error (MaxErr), the mean absolute error (MAE) and mean signed error (MSE) are given with respect to CCSD(T). The same data are reproduced in the SI at the third decimal place for reference. We further provide for information the data associated with the simple "gap-self-consistent" (egGW ) calculations.

		CCSD(T)	G 0 W 0		egGW
		IP	EA	IP	EA	IP	EA
	Anthracene	7.47 0.26 7.16 0.56	7.30 0.37
	Acridine	N/A N/A 7.64 0.90	7.81 0.72
	Phenazine	8.42 1.03 8.05 1.32	8.25 1.15
	Azulene	7.49 0.48 7.20 0.72	7.36 0.54
	Benzoquinone	10.17 1.46 9.81 1.75 10.41 1.54
	Naphthalenedione	9.79 1.39 9.41 1.68 10.01 1.45
	Dichlone	9.89 1.82 9.48 2.12	9.72 1.90
	F4-benzoquinone 11.04 2.18 10.68 2.48 10.95 2.29
	Cl4-benzoquinone 10.12 2.36 9.79 2.67 10.03 2.46
	Nitrobenzene	10.14 0.44 9.85 0.75 10.06 0.48
	F4-benzenedicar. 10.66 1.51 10.32 1.84 10.57 1.67
	Dinitrobenzoni.	11.07 1.67 10.81 2.02 11.06 1.75
	Nitrobenzoni.	10.55 1.20 10.26 1.56 10.49 1.32
	Benzonitrile	9.88 -0.29 9.57 -0.02	9.78 -0.24
	Fumaronitrile	11.40 0.89 11.03 1.21 11.28 0.98
	mDCNB	10.37 0.54 10.08 0.86 10.30 0.63
	TCNE	11.90 2.94 11.52 3.34 11.77 3.17
	TCNQ	9.49 3.23 9.20 3.67	9.35 3.56
	Maleicanhydride	11.27 0.92 10.85 1.15 11.45 0.92
	Phthalimide	10.02 0.54 9.75 0.79 10.36 0.56
	phthalicanhydride 10.48 0.78 10.18 1.04 10.39 0.81
	Cl4-isobenz.	N/A N/A 9.65 1.84	9.87 1.61
	NDCA	N/A N/A 8.78 1.46	8.94 1.27
	bodipy	N/A N/A 7.88 1.76	8.00 1.61
	MaxErr			-0.41 0.44	0.34 0.33
	MAE			0.33 0.31	0.14 0.10
	MSE			-0.33 0.31 -0.04 0.10
	states.

Table S1 :

 S1 [START_REF] Aryasetiawan | [END_REF] ) if the exact spectral representation of W is obtained by full diagonalization Ionization potential (defined as -E HOM O ) for the GW 100 test set as obtained with the present analytic continuation scheme. Reference calculations (ref.

	GW100 IPs Continuation of GW100 IPs Continuation of GW100 IPs
		G 0 W 0 @PBE G 0 W 0 @PBE G 0 W 0 @PBE	G 0 W 0 @PBE0 G 0 W 0 @PBE0 G 0 W 0 @PBE0
	i Formula i Formula i Formula	ref. ref. ref.	P-14 P-14 P-14	ref. ref. ref.	P-14 P-14 P-14
	1 He 47 NH 3 89 C 6 H 5 OH	23.477 23.483 23.995 23.999 10.315 10.315 10.698 10.698 8.263 8.263 8.502 8.502
	2 Ne 48 HN 3 90 C 6 H 5 NH 2	20.376 20.376 20.879 20.879 10.395 10.395 10.686 10.686 7.644 7.644 7.914 7.914
	3 Ar 49 PH 3 91 C 5 H 5 N	15.132 15.132 15.417 15.417 10.273 10.273 10.487 10.487 9.040 9.040 9.600 9.600
	4 Kr 50 AsH 3 92 C 5 H 5 N 5 O	13.571 13.571 13.781 13.781 10.123 10.123 10.300 10.300 7.691 7.691 7.976 7.976
	6 H2 51 SH 2 93 C 5 H 5 N 5 O	15.815 15.815 16.182 16.182 10.031 10.031 10.248 10.248 7.975 7.975 8.256 8.256
	7 Li 2 52 FH 94 C 4 H 5 N 3 O	4.992 15.302 15.302 15.748 15.748 4.992 5.247 5.247 8.287 8.287 8.683 8.683
	8 Na 2 53 ClH 95 C 5 H 6 N 2 O 2	4.837 12.246 12.246 12.501 12.501 4.837 4.966 4.966 8.707 8.707 9.057 9.057
	9 Na 4 54 LiF 96 C 4 H 4 N 2 O 2	4.120 9.946 9.223	4.120 9.947 10.829 10.829 4.233 4.233 9.223 9.457 9.457
	10 Na 6 55 F 2 Mg 97 CH 4 N 2 O	4.242 12.323 12.323 13.188 13.188 4.242 4.356 4.356 9.321 9.321 9.908 9.908
	11 K 2 56 TiF 4 99 Cu 2	3.980 13.896 13.896 14.868 14.868 3.981 4.072 4.072 7.529 7.529 7.551 7.551
	13 N 2 57 AlF 3 100 CuCN	14.889 14.889 15.429 15.429 14.252 14.252 14.950 14.950 9.423 9.425 10.253 10.253
	of the susceptibility. 14 P 2 58 BF	10.205 10.205 10.373 10.373 10.562 10.562 10.909 10.909
	15 As 2 59 SF 4	9.471 12.117 12.117 12.558 12.558 9.471 9.599 9.599
	16 F 2 60 BrK (84) Nozieres, P.; Pines, D. Theory Of Quantum Liquids; Advanced Books Classics; Avalon 14.962 14.962 15.529 15.529 7.304 7.305 7.946 7.946 17 Cl 2 11.101 11.101 11.384 11.384 61 GaCl 9.554 9.554 9.719 9.719
	Publishing, 1999. 18 Br 2 62 NaCl	10.217 10.217 10.451 10.451 8.098 8.099 8.835 8.835
	20 CH 4 63 MgCl 2	13.927 13.927 14.267 14.267 10.991 10.991 11.462 11.462
	21 C 2 H 6 65 BN	12.365 12.365 12.676 12.676 11.010 11.010 11.543 11.543
	22 C 3 H 8 66 NCH	11.795 11.795 12.106 12.106 13.211 13.211 13.525 13.525
	23 C 4 H 10 67 PN	11.490 11.490 11.804 11.804 11.136 11.136 11.688 11.688
	24 C 2 H 4 68 H 2 NNH 2	10.325 10.325 10.514 10.514 9.281 9.281 9.647 9.647
	25 C 2 H 2 69 H 2 CO	11.020 11.020 11.271 11.271 10.328 10.328 10.784 10.784
	26 C 4 70 CH 4 O	10.781 10.781 11.156 11.156 10.562 10.562 10.974 10.974
	27 C 3 H 6 71 C 2 H 6 O	10.555 10.555 10.821 10.821 10.157 10.157 10.604 10.604
	28 C 6 H 6 72 C 2 H 4 O	8.985 9.548	8.985 9.548 10.094 10.094 9.196 9.196
	29 C 8 H 8 73 C 4 H 10 O	8.055 9.319	8.055 9.319	8.284 9.777	8.284 9.777
	30 C 5 H 6 74 CH 2 O 2	8.349 10.730 10.730 11.292 11.292 8.349 8.551 8.551
	31 C 2 H 3 F 75 HOOH	10.197 10.197 10.441 10.441 10.987 10.987 11.444 11.444
	32 C 2 H 3 Cl 76 H 2 O	9.761 11.973 11.973 12.387 12.387 9.761 10.008 10.008
	33 C 2 H 3 Br 77 CO 2	8.990 13.250 13.250 13.662 13.662 8.990 9.184 9.184
	33 C 2 H 3 Br 78 CS 2	9.491 9.748	9.491 9.748	9.728 9.984	9.728 9.984
	35 CF 4 79 CSO	15.366 15.366 15.995 15.995 10.906 10.906 11.170 11.170
	36 CCl 4 80 CSeO	10.975 10.975 11.401 11.401 10.199 10.199 10.413 10.413
	37 CBr 4 81 CO	9.899 13.571 13.571 14.120 14.120 9.899 10.272 10.272
	39 SiH 4 82 O 3	12.310 12.310 12.699 12.699 11.967 11.967 12.563 12.563
	40 GeH 4 83 SO 2	12.022 12.022 12.358 12.358 11.823 11.823 12.272 12.272
	41 Si 2 H 6 84 BeO	10.309 10.309 10.585 10.585 9.633 9.635 9.644 9.644
	42 Si 5 H 12 85 MgO	8.934 6.671	8.934 6.671	9.220 7.457	9.220 7.457
	43 LiH 86 C 7 H 8	6.552 8.612	6.552 8.612	7.547 8.818	7.547 8.818
	44 KH 87 C 8 H 10	4.862 8.553	4.862 8.553	5.660 8.766	5.660 8.766
	45 BH 3 88 C 6 F 6	12.867 12.867 13.225 13.225 9.493 9.493 9.893 9.893
	46 B 2 H 6 89 C 6 H 5 OH	11.835 11.835 12.218 12.218 8.366 8.366 8.609 8.609

Table S2 :

 S2 Same as in TableS1but for the electronic affinities (defined as -E LU M O ).

	Continuation of GW100 EAs Continuation of GW100 EAs
		G 0 W 0 @PBE G 0 W 0 @PBE	G 0 W 0 @PBE0 G 0 W 0 @PBE0
	i Formula i Formula	ref. ref.	P-14 P-14	ref. ref.	P-14 P-14
	18 Br 2 62 NaCl	1.396 0.394	1.396 0.394	1.224 0.433	1.224 0.433
	20 CH 4 63 MgCl 2	-2.450 0.430	-2.450 0.430	-2.486 0.333	-2.486 0.333
	21 C 2 H 6 65 BN	-2.294 3.947	-2.294 3.947	-2.357 3.760	-2.357 3.760
	22 C 3 H 8 66 NCH	-2.195 -2.577	-2.195 -2.577	-2.278 -2.674	-2.278 -2.674
	23 C 4 H 10 67 PN	-2.145 0.202	-2.145 0.202	-2.244 0.084	-2.244 0.084
	24 C 2 H 4 68 H 2 NNH 2	-2.020 -1.993	-2.020 -1.993	-2.189 -2.062	-2.189 -2.062
	25 C 2 H 2 69 H 2 CO	-2.863 -0.959	-2.863 -0.959	-2.965 -1.192	-2.965 -1.192
	26 C 4 70 CH 4 O	2.936 -2.247	2.936 -2.247	2.710 -2.264	2.710 -2.264
	27 C 3 H 6 71 C 2 H 6 O	-2.453 -2.076	-2.453 -2.076	-2.525 -2.130	-2.525 -2.130
	28 C 6 H 6 72 C 2 H 4 O	-1.087 -1.050	-1.087 -1.050	-1.300 -1.375	-1.300 -1.375
	29 C 8 H 8 73 C 4 H 10 O	-0.059 -2.100	-0.059 -2.100	-0.316 -2.199	-0.316 -2.199
	30 C 5 H 6 74 CH 2 O 2	-1.040 -1.910	-1.040 -1.910	-1.274 -2.172	-1.274 -2.172
	31 C 2 H 3 F 75 HOOH	-2.146 -2.349	-2.146 -2.349	-2.317 -2.548	-2.317 -2.548
	32 C 2 H 3 Cl 76 H 2 O	-1.425 -2.370	-1.425 -2.370	-1.623 -2.377	-1.623 -2.377
	33 C 2 H 3 Br 77 CO 2	-1.380 -2.497	-1.380 -2.497	-1.571 -2.509	-1.571 -2.509
	33 C 2 H 3 Br 78 CS 2	-1.257 0.197	-1.257 0.197	-1.472 0.104	-1.472 0.104
	35 CF 4 79 CSO	-4.413 -1.215	-4.413 -1.215	-4.479 -1.360	-4.479 -1.360
	36 CCl 4 80 CSeO	0.008 -0.868	0.008 -0.868	-0.245 -1.024	-0.245 -1.024
	37 CBr 4 81 CO	1.079 -0.671	1.079 -0.671	0.856 -0.797	0.856 -0.797
	39 SiH 4 82 O 3	-2.508 2.296	-2.508 2.296	-2.560 2.212	-2.560 2.212
	40 GeH 4 83 SO 2	-2.300 1.002	-2.300 1.002	-2.391 0.883	-2.391 0.883
	41 Si 2 H 6 84 BeO 2 GW 100 test set electronic affinities -1.686 -1.686 -1.889 2.493 2.493 2.161 42 Si 5 H 12 -0.161 -0.161 -0.416 85 MgO 1.893 1.893 1.737	-1.889 2.161 -0.416 1.737
	43 LiH 86 C 7 H 8	0.072 -1.016	0.072 -1.016	0.092 -1.237	0.092 -1.237
	44 KH 87 C 8 H 10	0.180 -1.044	GW100 EAs 0.180 0.154 -1.044 -1.279	0.154 -1.279
	45 BH 3 88 C 6 F 6	G 0 W 0 @PBE -0.115 -0.115 -0.657 -0.657	G 0 W 0 @PBE0 -0.263 -0.263 -0.633 -0.633
	i Formula 46 B 2 H 6 89 C 6 H 5 OH	ref. -0.838 -0.959	P-14 -0.838 -0.959	ref. -1.014 -1.155	P-14 -1.014 -1.155
	1 He 47 NH 3 89 C 6 H 5 OH	-11.006 -11.007 -11.019 -11.020 -2.312 -2.312 -2.345 -2.345 -1.013 -1.013 -1.211 -1.211
	2 Ne 48 HN 3 90 C 6 H 5 NH 2	-11.643 -11.643 -11.674 -11.674 -1.399 -1.399 -1.549 -1.549 -1.148 -1.148 -1.358 -1.358
	3 Ar 49 PH 3 91 C 5 H 5 N	-8.113 -2.496 -0.514	-8.113 -2.496 -0.514	-8.166 -2.096 -0.735	-8.166 -2.096 -0.735
	4 Kr 50 AsH 3 92 C 5 H 5 N 5 O	-7.635 -2.320 -0.749	-7.635 -2.320 -0.749	-7.689 -2.015 -1.008	-7.689 -2.015 -1.008
	6 H2 51 SH 2 93 C 5 H 5 N 5 O	-3.504 -2.557 -0.474	-3.504 -2.557 -0.474	-3.436 -2.024 -0.735	-3.436 -2.024 -0.735
	7 Li 2 52 FH 94 C 4 H 5 N 3 O	0.626 -2.543 -0.263	0.626 -2.543 -0.263	0.409 -2.483 -0.476	0.409 -2.483 -0.476
	8 Na 2 53 ClH 95 C 5 H 6 N 2 O 2	0.550 -2.064 -0.057	0.550 -2.064 -0.057	0.421 -2.051 -0.282	0.421 -2.051 -0.282
	9 Na 4 54 LiF 96 C 4 H 4 N 2 O 2	1.009 -0.091 -0.012	1.009 -0.091 -0.012	0.779 0.049 -0.229	0.779 0.049 -0.229
	10 Na 6 55 F 2 Mg 97 CH 4 N 2 O	0.971 0.138 -1.621	0.971 0.138 -1.621	0.762 0.084 -1.702	0.762 0.084 -1.702
	11 K 2 56 TiF 4 99 Cu 2	0.647 0.598 0.756	0.647 0.598 0.756	0.508 0.675 0.595	0.508 0.675 0.595
	13 N 2 57 AlF 3 100 CuCN	-2.449 -0.157 1.647	-2.449 -0.157 1.647	-2.586 -0.271 1.264	-2.586 -0.271 1.264
	14 P 2 58 BF	0.724 -1.216	0.724 -1.216	0.608 -1.298	0.608 -1.298
	15 As 2 59 SF 4	0.847 -0.377	0.847 -0.377	0.753 -0.543	0.753 -0.543
	16 F 2 60 BrK	0.704 0.307	0.704 0.307	0.308 0.362	0.308 0.362
	17 Cl 2 61 GaCl	0.887 0.020	0.887 0.020	0.687 -0.062	0.687 -0.062

Supporting Information Available

The following files are available free of charge.

We provide in the Supporting Information the calculated ionization potential and electronic affinity for the GW 100 test set and the 24 medium size acceptor molecular set introduced in Ref. 58. We further provide the AC errors as in Fig. 7 but on a larger energy range that contains all occupied and unoccupied states.

Graphical TOC Entry

The analytic continuation of the screened Coulomb potential from few selected data points yields accurate GW self-energies over large energy ranges.

Supplementary Information :

Robust analytic continuation approach to many-body GW calculations Ivan Duchemin * , † and Xavier Blase ‡ †Univ. Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France ‡Univ. Grenoble Alpes, CNRS, Inst NEEL, F-38042 Grenoble, France E-mail: ivan.duchemin@cea.fr

We provide here below the ionization potentials (Section 1) and electronic affinities (Section 2) for the GW100 test set, S1 comparing reference calculations (Fig. 2b scheme, ∆ω = 125meV ) with data stemming from the analytic continuation of W relying on a reduced nω=14 data points sampling along the imaginary axis (Fig. 2a scheme). We consider only the 93 molecules not requiring pseudopotentials. Both def2-QZVP G 0 W 0 @PBE and def2-QZVP G 0 W 0 @PBE0 data are presented. Similarly, we list the aug-cc-pVTZ G 0 W 0 @PBE0 ionization potential and electronic affinities (Section 3) for the set of 24 medium size acceptor molecules from Refs. S2,S3. Finally, we reproduce Fig. 7b of the main manuscript, but over the full energy range and a lower Z = 0.2 cutoff value.

1 GW 100 test set ionisation potentials 1


est destinée au dépôt

  Comparaison of the biological potencies of the two tilapia prolactins on tiPRL receptor transfected in human fibroblastes (293) cell line Frédéric Sohm, A. Pezet, Olivier Sandra, M. Edery, Patrick Prunet, P. Kelly, B. Demeneix, A. de Luze To cite this version: Frédéric Sohm, A. Pezet, Olivier Sandra, M. Edery, Patrick Prunet, et al.. Comparaison of the biological potencies of the two tilapia prolactins on tiPRL receptor transfected in human fibroblastes (293) cell line. 10. Conference of European Comparative Endocrinologists from Molecular to Ingrative Biology, Sep 1996, Rouen, France. Masson, Annales d'Endocrinologie, 57 (4) supplément, 1996. hal-


est destinée au dépôt

Molecular characterization of prolactin receptor (PRLR) in fish

Olivier Sandra, Philippe Le Rouzic, Benoît Aupérin, M. Edery, P.A. Kelly,

Patrick Prunet

To cite this version: