Dirac-Coulomb-Breit Molecular Mean-Field Exact-Two-Component Relativistic Equation-of-Motion Coupled-Cluster Theory.

We present a relativistic equation-of-motion coupled-cluster with single and double excitation formalism within the exact two-component framework (X2C-EOM-CCSD), where both scalar relativistic effects and spin-orbit coupling are variationally included at the reference level. Three different molecular mean-field treatments of relativistic corrections, including the one-electron, Dirac-Coulomb, and Dirac-Coulomb-Breit Hamiltonian, are considered in this work. Benchmark calculations include atomic excitations and fine-structure splittings arising from spin-orbit coupling. Comparison with experimental values and relativistic time-dependent density functional theory is also carried out. The computation of the oscillator strength using the relativistic X2C-EOM-CCSD approach allows for studies of spin-orbit-driven processes, such as the spontaneous phosphorescence lifetime.

Relativistic resolution-of-the-identity with Cholesky integral decomposition.

In this study, we present an efficient integral decomposition approach called the restricted-kinetic-balance resolution-of-the-identity (RKB-RI) algorithm, which utilizes a tunable RI method based on the Cholesky integral decomposition for in-core relativistic quantum chemistry calculations. The RKB-RI algorithm incorporates the restricted-kinetic-balance condition and offers a versatile framework for accurate computations. Notably, the Cholesky integral decomposition is employed not only to approximate symmetric large-component electron repulsion integrals but also those involving small-component basis functions. In addition to comprehensive error analysis, we investigate crucial conditions, such as the kinetic balance condition and variational stability, which underlie the applicability of Dirac relativistic electronic structure theory. We compare the computational cost of the RKB-RI approach with the full in-core method to assess its efficiency. To evaluate the accuracy and reliability of the RKB-RI method proposed in this work, we employ actinyl oxides as benchmark systems, leveraging their properties for validation purposes. This investigation provides valuable insights into the capabilities and performance of the RKB-RI algorithm and establishes its potential as a powerful tool in the field of relativistic quantum chemistry.

Electrochemical Strategy for Proton Relay Installation Enhances the Activity of a Hydrogen Evolution Electrocatalyst.

A new method to install a proton relay that enhances the reactivity near an active catalytic site for H2 production is reported, afforded by the electrochemical reduction and protonation of one of the ligands in the paddlewheel Rh2(II,II) hydrogen evolution complex, cis-[Rh2(DPhF)2(bncn)2]2+ (Rh-bncn; DPhF = N,N'-diphenylformamidinate, bncn = benzo[c]cinnoline). An electrochemical reversible prewave is observed for Rh-bncn at potentials more positive than the first bncn-centered reduction couple in the presence of strong acids, observed at -0.72 V vs Fc+/0 (Fc = ferrocene) in the cyclic voltammograms (CVs) in DMF (0.1 M TBAPF6). The origin of this prewave is shown to arise from a precatalytic transformation that originates from a concerted proton-electron transfer (CPET) event occurring at one of the bridging bncn ligands. Through electrochemical analysis, CV simulations, and electronic structure calculations, a reaction mechanism is elucidated. In this system, the electrochemically formed N-H bond on the reduced bncn ligand serves as a proton relay in the H2 formation reaction through a cooperative interligand pathway involving one of the bridging DPhF ligands after a second reduction step, accessible at approximately -1.15 V vs Fc+/0. Since calculations show that hydrogen evolution takes place at the bridging ligands and does not involve the dirhodium core, it is predicted that more abundant metal centers can be incorporated into this ligand scaffold, leading to new candidates for electrocatalytic hydrogen reduction. As such, this work delineates a new design strategy to incorporate proton relays in molecular bimetallic hydrogen evolution electrocatalysts to achieve higher efficiency.

Quantifying and reducing spin contamination in algebraic diagrammatic construction theory of charged excitations.

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 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.

Efficient implementation of the single-reference algebraic diagrammatic construction theory for charged excitations: Applications to the TEMPO radical and DNA base pairs.

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.

Recent developments in the PySCF program package.

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. This 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.

Third-order algebraic diagrammatic construction theory for electron attachment and ionization energies: Conventional and Green's function implementation.

We present implementation of second- and third-order algebraic diagrammatic construction (ADC) 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 single-particle 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 the 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 the efficiency of our GF-ADC implementation by computing core-level K-, L-, and M-shell ionization energies of a zinc atom without introducing the core-valence separation approximation. Finally, we use EA- and IP-ADC methods to compute the bandgaps of equally spaced hydrogen chains Hn 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.

Algebraic Diagrammatic Construction Theory for Simulating Charged Excited States and Photoelectron Spectra.

Charged excitations are electronic transitions that involve a change in the total charge of a molecule or material. Understanding the properties and reactivity of charged species requires insights from theoretical calculations that can accurately describe orbital relaxation and electron correlation effects in open-shell electronic states. In this Review, we describe the current state of algebraic diagrammatic construction (ADC) theory for simulating charged excitations and its recent developments. We start with a short overview of ADC formalism for the one-particle Green's function, including its single- and multireference formulations and extension to periodic systems. Next, we focus on the capabilities of ADC methods and discuss recent findings about their accuracy for calculating a wide range of excited-state properties. We conclude our Review by outlining possible directions for future developments of this theoretical approach.

Non-Dyson Algebraic Diagrammatic Construction Theory for Charged Excitations in Solids.

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).