PYSCF is a Python-based general-purpose electronic structure platform that both supports first-principles simulations of molecules and solids, as well as accelerates the development of new methodology and complex computational workflows. The present paper explains the design and philosophy behind PYSCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PYSCF as a development environment. We then summarize the capabilities of PYSCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PYSCF across the domains of quantum chemistry, materials science, machine learning and quantum information science.

We present implementation of second-and third-order algebraic diagrammatic construction theory for efficient and accurate computations of molecular electron affinities (EA), ionization potentials (IP), and densities of states (EA-/IP-ADC(n), n = 2, 3). Our work utilizes the non-Dyson formulation of ADC for the singleparticle propagator and reports working equations and benchmark results for the EA-ADC(2) and EA-ADC(3) approximations. We describe two algorithms for solving EA-/IP-ADC equations: (i) conventional algorithm that uses iterative diagonalization techniques to compute low-energy EA, IP, and density of states, and (ii) Green's function algorithm (GF-ADC) that solves a system of linear equations to compute density of states directly for a specified spectral region. To assess accuracy of EA-ADC(2) and EA-ADC(3), we benchmark their performance for a set of atoms, small molecules, and five DNA/RNA nucleobases. As our next step, we demonstrate efficiency of our GF-ADC implementation by computing core-level K-, L-, and M -shell ionization energies of a zinc atom without introducing core-valence separation approximation. Finally, we use EA-and IP-ADC methods to compute band gaps of equally-spaced hydrogen chains H n with n up to 150, providing their estimates near thermodynamic limit. Our results demonstrate that EA-/IP-ADC(n) (n = 2, 3) methods are efficient and accurate alternatives to widely used electronic structure methods for simulations of electron attachment and ionization properties.

We present an efficient implementation of the second- and third-order single-reference algebraic diagrammatic construction (ADC) theory for electron attachment and ionization energies and spectra [EA/IP-ADC( n), n = 2, 3]. Our new EA/IP-ADC program features spin adaptation for closed-shell systems, density fitting for efficient handling of the two-electron integral tensors, and vectorized and parallel implementation of tensor contractions. We demonstrate capabilities of our efficient implementation by applying the EA/IP-ADC( n) (n = 2, 3) methods to compute the photoelectron spectrum of the (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) radical, as well as the vertical and adiabatic electron affinities of TEMPO and two DNA base pairs (guanine–cytosine and adenine–thymine). The spectra and electron affinities computed using large diffuse basis sets with up to 1028 molecular orbitals are found to be in good agreement with the best available results from the experiment and theoretical simulations.

Algebraic diagrammatic construction (ADC) theory is a computationally efficient and accurate approach for simulating electronic excitations in chemical systems. However, for the simulations of excited states in molecules with unpaired electrons the performance of ADC methods can be affected by the spin contamination in unrestricted Hartree-Fock (UHF) reference wavefunctions. In this work, we benchmark the accuracy of ADC methods for electron attachment and ionization of open-shell molecules with the UHF reference orbitals (EA/IP-ADC/UHF) and develop an approach to quantify the spin contamination in the charged excited states. Following this assessment, we demonstrate that the spin contamination can be reduced by combining EA/IP-ADC with the reference orbitals from restricted open-shell Hartree-Fock (ROHF) or orbital-optimized Møller-Plesset perturbation (OMP) theories. Our numerical results demonstrate that for open-shell systems with strong spin contamination in the UHF reference the third-order EA/IP-ADC methods with the ROHF or OMP reference orbitals are similar in accuracy to equation-of-motion coupled cluster theory with single and double excitations.

We present the first implementation and applications of non-Dyson algebraic diagrammatic construction theory for charged excitations in three-dimensional periodic solids (EA/IP-ADC). The EA/IP-ADC approach has a computational cost similar to the ground-state Møller–Plesset perturbation theory, enabling efficient calculations of a variety of crystalline excited-state properties (e.g., band structure, band gap, density of states) sampled in the Brillouin zone. We use EA/IP-ADC to compute the quasiparticle band structures and band gaps of several materials (from large-gap atomic and ionic solids to small-gap semiconductors) and analyze the errors of EA/IP-ADC approximations up to the third order in perturbation theory. Our work also reports the first-ever calculations of ground-state properties (equation-of-state and lattice constants) of three-dimensional crystalline systems using a periodic implementation of third-order Møller–Plesset perturbation theory (MP3).