Understanding ion-transfer reactions in silver electrodissolution and electrodeposition from first-principles calculations and experiments.

The electrified aqueous/metal interface is critical in controlling the performance of energy conversion and storage devices, but an atomistic understanding of even basic interfacial electrochemical reactions challenges both experiment and computation. We report a combined simulation and experimental study of (reversible) ion-transfer reactions involved in anodic Ag corrosion/deposition, a model system for interfacial electrochemical processes generating or consuming ions. With the explicit modeling of the electrode potential and a hybrid implicit-explicit solvation model, the density functional theory calculations produce free energy curves predicting thermodynamics, kinetics, partial charge profiles, and reaction trajectories. The calculated (equilibrium) free energy barriers (0.2 eV), and their asymmetries, agree with experimental activation energies (0.4 eV) and transfer coefficients, which were extracted from temperature-dependent voltage-step experiments on Au-supported, Ag-nanocluster substrates. The use of Ag nanoclusters eliminates the convolution of the kinetics of Ag+(aq.) generation and transfer with those of nucleation or etch-pit formation. The results indicate that the barrier is controlled by the bias-dependent competition between partial solvation of the incipient ion, metal-metal bonding, and electrostatic stabilization by image charge, with the latter two factors weakened by stronger positive biases. We also report simulations of the bias-dependence of defect generation relevant to nucleating corrosion by removing an atom from a perfect Ag(100) surface, which is predicted to occur via a vacancy-adatom intermediate. Together, these experiments and calculations provide the first validated, accurate, molecular model of the central steps that govern the rates of important dissolution/deposition reactions broadly relevant across the energy sciences.

Chemical Bonding and the Role of Node-Induced Electron Confinement.

The chemical bond is the cornerstone of chemistry, providing a conceptual framework to understand and predict the behavior of molecules in complex systems. However, the fundamental origin of chemical bonding remains controversial and has been responsible for fierce debate over the past century. Here, we present a unified theory of bonding, using a separation of electron delocalization effects from orbital relaxation to identify three mechanisms [node-induced confinement (typically associated with Pauli repulsion, though more general), orbital contraction, and polarization] that each modulate kinetic energy during bond formation. Through analysis of a series of archetypal bonds, we show that an exquisite balance of energy-lowering delocalizing and localizing effects are dictated simply by atomic electron configurations, nodal structure, and electronegativities. The utility of this unified bonding theory is demonstrated by its application to explain observed trends in bond strengths throughout the periodic table, including main group and transition metal elements.

Mechanistic Insights into the Origins of Selectivity in a Cu-Catalyzed C-H Amidation Reaction.

The catalytic transformation of C-H to C-N bonds offers rapid access to fine chemicals and high-performance materials, but achieving high selectivity from undirected aminations of unactivated C(sp3)-H bonds remains an outstanding challenge. We report the origins of the reactivity and selectivity of a Cu-catalyzed C-H amidation of simple alkanes. Using a combination of experimental and computational mechanistic studies and energy decomposition techniques, we uncover a switch in mechanism from inner-sphere to outer-sphere coupling between alkyl radicals and the active Cu(II) catalyst with increasing substitution of the alkyl radical. The combination of computational predictions and detailed experimental validation shows that simultaneous minimization of both Cu-C covalency and alkyl radical size increases the rate of reductive elimination and that both strongly electron-donating and electron-withdrawing substituents on the catalyst accelerate the selectivity-determining C-N bond formation process as a result of a change in mechanism. These findings offer design principles for the development of improved catalyst scaffolds for radical C-H functionalization reactions.

Highly Accurate Prediction of NMR Chemical Shifts from Low-Level Quantum Mechanics Calculations Using Machine Learning.

Theoretical predictions of NMR chemical shifts from first-principles can greatly facilitate experimental interpretation and structure identification of molecules in gas, solution, and solid-state phases. However, accurate prediction of chemical shifts using the gold-standard coupled cluster with singles, doubles, and perturbative triple excitations [CCSD(T)] method with a complete basis set (CBS) can be prohibitively expensive. By contrast, machine learning (ML) methods offer inexpensive alternatives for chemical shift predictions but are hampered by generalization to molecules outside the original training set. Here, we propose several new ideas in ML of the chemical shift prediction for H, C, N, and O that first introduce a novel feature representation, based on the atomic chemical shielding tensors within a molecular environment using an inexpensive quantum mechanics (QM) method, and train it to predict NMR chemical shieldings of a high-level composite theory that approaches the accuracy of CCSD(T)/CBS. In addition, we train the ML model through a new progressive active learning workflow that reduces the total number of expensive high-level composite calculations required while allowing the model to continuously improve on unseen data. Furthermore, the algorithm provides an error estimation, signaling potential unreliability in predictions if the error is large. Finally, we introduce a novel approach to keep the rotational invariance of the features using tensor environment vectors (TEVs) that yields a ML model with the highest accuracy compared to a similar model using data augmentation. We illustrate the predictive capacity of the resulting inexpensive shift machine learning (iShiftML) models across several benchmarks, including unseen molecules in the NS372 data set, gas-phase experimental chemical shifts for small organic molecules, and much larger and more complex natural products in which we can accurately differentiate between subtle diastereomers based on chemical shift assignments.

Quantum chemical modeling of hydrogen binding in metal-organic frameworks: validation, insight, predictions and challenges.

A detailed chemical understanding of H2 interactions with binding sites in the nanoporous crystalline structure of metal-organic frameworks (MOFs) can lay a sound basis for the design of new sorbent materials. Computational quantum chemical calculations can aid in this quest. To set the stage, we review general thermodynamic considerations that control the usable storage capacity of a sorbent. We then discuss cluster modeling of H2 ligation at MOF binding sites using state-of-the-art density functional theory (DFT) calculations, and how the binding can be understood using energy decomposition analysis (EDA). Employing these tools, we illustrate the connections between the character of the MOF binding site and the associated adsorption thermodynamics using four experimentally characterized MOFs, highlighting the role of open metal sites (OMSs) in accessing binding strengths relevant to room temperature storage. The sorbents are MOF-5, with no open metal sites, Ni2(m-dobdc), containing Lewis acidic Ni(II) sites, Cu(I)-MFU-4l, containing π basic Cu(I) sites and V2Cl2.8(btdd), also containing π-basic V(II) sites. We next explore the potential for binding multiple H2 molecules at a single metal site, with thermodynamics useful for storage at ambient temperature; a materials design goal which has not yet been experimentally demonstrated. Computations on Ca2+ or Mg2+ bound to catecholate or Ca2+ bound to porphyrin show the potential for binding up to 4 H2; there is precedent for the inclusion of both catecholate and porphyrin motifs in MOFs. Turning to transition metals, we discuss the prediction that two H2 molecules can bind at V(II)-MFU-4l, a material that has been synthesized with solvent coordinated to the V(II) site. Additional calculations demonstrate binding three equivalents of hydrogen per OMS in Sc(I) or Ti(I)-exchanged MFU-4l. Overall, the results suggest promising prospects for experimentally realizing higher capacity hydrogen storage MOFs, if nontrivial synthetic and desolvation challenges can be overcome. Coupled with the unbounded chemical diversity of MOFs, there is ample scope for additional exploration and discovery.

Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone.

We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.

Selective Adsorption of Oxygen from Humid Air in a Metal-Organic Framework with Trigonal Pyramidal Copper(I) Sites.

High or enriched-purity O2 is used in numerous industries and is predominantly produced from the cryogenic distillation of air, an extremely capital- and energy-intensive process. There is significant interest in the development of new approaches for O2-selective air separations, including the use of metal-organic frameworks featuring coordinatively unsaturated metal sites that can selectively bind O2 over N2 via electron transfer. However, most of these materials exhibit appreciable and/or reversible O2 uptake only at low temperatures, and their open metal sites are also potential strong binding sites for the water present in air. Here, we study the framework CuI-MFU-4l (CuxZn5-xCl4-x(btdd)3; H2btdd = bis(1H-1,2,3-triazolo[4,5-b],[4',5'-i])dibenzo[1,4]dioxin), which binds O2 reversibly at ambient temperature. We develop an optimized synthesis for the material to access a high density of trigonal pyramidal CuI sites, and we show that this material reversibly captures O2 from air at 25 °C, even in the presence of water. When exposed to air up to 100% relative humidity, CuI-MFU-4l retains a constant O2 capacity over the course of repeated cycling under dynamic breakthrough conditions. While this material simultaneously adsorbs N2, differences in O2 and N2 desorption kinetics allow for the isolation of high-purity O2 (>99%) under relatively mild regeneration conditions. Spectroscopic, magnetic, and computational analyses reveal that O2 binds to the copper(I) sites to form copper(II)-superoxide moieties that exhibit temperature-dependent side-on and end-on binding modes. Overall, these results suggest that CuI-MFU-4l is a promising material for the separation of O2 from ambient air, even without dehumidification.

Doubles Connected Moments Expansion: A Tractable Approximate Horn-Weinstein Approach for Quantum Chemistry.

Ab initio methods based on the second-order and higher connected moments, or cumulants, of a reference function have seen limited use in the determination of correlation energies of chemical systems over the years. Moment-based methods have remained unattractive relative to more ubiquitous methods, such as perturbation theory and coupled cluster theory, due in part to the intractable cost of assembling moments of high-order and poor performance of low-order expansions. Many of the traditional quantum chemical methodologies can be recast as a selective summation of perturbative contributions to their energy; using this familiar structure as a guide in selecting terms, we develop a scheme to approximate connected moments limited to double excitations. The tractable Doubles Connected Moments [DCM(N)] approximation is developed and tested against a multitude of common single-reference methods to determine its efficacy in the determination of the correlation energy of model systems and small molecules. The DCM(N) sequence of energies exhibits smooth convergence toward limiting values in the range of N = 11-14, with compute costs that scale as a noniterative O(M6) with molecule size, M. Numerical tests on correlation energy recovery for 55 small molecules comprising the G1 test set in the cc-pVDZ basis show that DCM(N) strongly outperforms MP2 and even CCD with a Hartree-Fock reference. When using an approximate Brueckner reference from orbital-optimized (oo) MP2, the resulting oo:DCM(N) energies converge to values more accurate than CCSD for 49 of 55 molecules. The qualitative success of the method in regions where strong correlation effects begin to dominate, even while maintaining spin purity, suggests this may be a good starting point in the development of methodologies for the description of strongly correlated or spin-contaminated systems while maintaining a tractable single-reference formalism.

Exactly Fragment Additive Breakdown of Polarization for Energy Decomposition Analysis Based on the Self-Consistent Field for Molecular Interactions.

Energy decomposition analysis (EDA) is a useful method to unravel intermolecular interaction energy into chemically meaningful components such as geometric distortion, frozen interactions, polarization, and charge transfer. A further decomposition of the polarization (POL) and charge transfer (CT) energy into fragment-wise contributions would be useful to understand the significance of each fragment during these two processes. To complement the existing exact pairwise decomposition of the CT term, this work describes the formulation and implementation of a nonperturbative polarization analysis that decomposes the POL energy into an exactly fragment-wise additive sum based on the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA). These fragment-wise contributions can be further decomposed into chemically intuitive molecular orbital pairs using complementary occupied-virtual pair (COVP) analysis. A very useful phase convention is established for each COVP such that constructive interference of occupied and virtual corresponds to electron flow into that region, while destructive interference corresponds to electron outflow. A range of model problems are used to demonstrate that the polarization process is typically a collective behavior of the electrons that is quite different from the charge transfer process. This provides another reason in addition to their different distance dependence on fragment separation for separating these two processes in EDA.

Fantastical excited state optimized structures and where to find them.

The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.

Earth Mover's Distance as a Metric to Evaluate the Extent of Charge Transfer in Excitations Using Discretized Real-Space Densities.

This paper presents a novel theoretical measure, µEMD, based on the earth mover's distance (EMD), for quantifying the density shift caused by electronic excitations in molecules. As input, the EMD metric uses only the discretized ground- and excited-state electron densities in real space, rendering it compatible with almost all electronic structure methods used to calculate excited states. The EMD metric is compared against other popular theoretical metrics for describing the extent of electron-hole separation in a wide range of excited states (valence, Rydberg, charge transfer, etc.). The results showcase the EMD metric's effectiveness across all excitation types and suggest that it is useful as an additional tool to characterize electronic excitations. The study also reveals that µEMD can function as a promising diagnostic tool for predicting the failure of pure exchange-correlation functionals. Specifically, we show statistical relationships among the functional-driven errors, the exact exchange content within the functional, and the magnitude of µEMD values.

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory.

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter (α) was set to recover the exact dissociation limit of H2 in a minimal basis set. Alternatively α can be viewed as a regularization parameter, where each value of α represents a valid variant of BW-s2, which we denote as BW-s2(α). In this work, we semi-empirically optimize α for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal α parameter (α = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2(α) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2(α) is iterative O(N5), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

The Energetic Origins of Pi-Pi Contacts in Proteins.

Accurate potential energy models of proteins must describe the many different types of noncovalent interactions that contribute to a protein's stability and structure. Pi-pi contacts are ubiquitous structural motifs in all proteins, occurring between aromatic and nonaromatic residues and play a nontrivial role in protein folding and in the formation of biomolecular condensates. Guided by a geometric criterion for isolating pi-pi contacts from classical molecular dynamics simulations of proteins, we use quantum mechanical energy decomposition analysis to determine the molecular interactions that stabilize different pi-pi contact motifs. We find that neutral pi-pi interactions in proteins are dominated by Pauli repulsion and London dispersion rather than repulsive quadrupole electrostatics, which is central to the textbook Hunter-Sanders model. This results in a notable lack of variability in the interaction profiles of neutral pi-pi contacts even with extreme changes in the dielectric medium, explaining the prevalence of pi-stacked arrangements in and between proteins. We also find interactions involving pi-containing anions and cations to be extremely malleable, interacting like neutral pi-pi contacts in polar media and like typical ion-pi interactions in nonpolar environments. Like-charged pairs such as arginine-arginine contacts are particularly sensitive to the polarity of their immediate surroundings and exhibit canonical pi-pi stacking behavior only if the interaction is mediated by environmental effects, such as aqueous solvation.

Local Second-Order Møller-Plesset Theory with a Single Threshold Using Orthogonal Virtual Orbitals: Theory, Implementation, and Assessment.

It has long been clear that electron correlation methods exhibit unphysical compute scalings with molecular size, which has motivated the development of local correlation methods to discard effectively zero contributions in a controlled way to yield an approximate correlation energy. The ideal local correlation method should have a single numerical threshold that controls the dropping of terms with the ability to have that threshold set small enough so that the correlation energy is reproduced to enough significant figures such that the result is chemically identical. This work reports such a method for the second-order Møller-Plesset (MP2) theory. The theory, implementation, and testing of this local MP2 theory are reported. Thresholds ranging from 10-5 to 10-8 and basis sets ranging from split valence plus polarization through to quadruple-ζ are assessed for local MP2 calculations on a range of molecules, including linear chains and molecules with two- and three-dimensional character. The implementation is shared memory parallel via OpenMP and yields roughly 50% parallel efficiency with 16 cores for a large job. Considerable efforts were made to minimize memory demands, which increased as thresholds were tightened. A variety of relative energy calculations are presented as a function of threshold to provide some guidance to users on how to obtain adequate precision at a low compute cost. It is particularly clear that derivative properties require tighter thresholds in order to achieve an adequate precision.

Bridging physical intuition and hardware efficiency for correlated electronic states: the local unitary cluster Jastrow ansatz for electronic structure.

A prominent goal in quantum chemistry is to solve the molecular electronic structure problem for ground state energy with high accuracy. While classical quantum chemistry is a relatively mature field, the accurate and scalable prediction of strongly correlated states found, e.g., in bond breaking and polynuclear transition metal compounds remains an open problem. Within the context of a variational quantum eigensolver, we propose a new family of ansatzes which provides a more physically appropriate description of strongly correlated electrons than a unitary coupled cluster with single and double excitations (qUCCSD), with vastly reduced quantum resource requirements. Specifically, we present a set of local approximations to the unitary cluster Jastrow wavefunction motivated by Hubbard physics. As in the case of qUCCSD, exactly computing the energy scales factorially with system size on classical computers but polynomially on quantum devices. The local unitary cluster Jastrow ansatz removes the need for SWAP gates, can be tailored to arbitrary qubit topologies (e.g., square, hex, and heavy-hex), and is well-suited to take advantage of continuous sets of quantum gates recently realized on superconducting devices with tunable couplers. The proposed family of ansatzes demonstrates that hardware efficiency and physical transparency are not mutually exclusive; indeed, chemical and physical intuition regarding electron correlation can illuminate a useful path towards hardware-friendly quantum circuits.

Greater transferability and accuracy of norm-conserving pseudopotentials using nonlinear core corrections.

We present an investigation into the transferability of pseudopotentials (PPs) with a nonlinear core correction (NLCC) using the Goedecker, Teter, and Hutter (GTH) protocol across a range of pure GGA, meta-GGA and hybrid functionals, and their impact on thermochemical and non-thermochemical properties. The GTH-NLCC PP for the PBE density functional demonstrates remarkable transferability to the PBE0 and ωB97X-V exchange-correlation functionals, and relative to no NLCC, improves agreement significantly for thermochemical benchmarks compared to all-electron calculations. On the other hand, the B97M-rV meta-GGA functional performs poorly with the PBE-derived GTH-NLCC PP, which is mitigated by reoptimizing the NLCC parameters for this specific functional. The findings reveal that atomization energies exhibit the greatest improvements from use of the NLCC, which thus provides an important correction needed for covalent interactions relevant to applications involving chemical reactivity. Finally we test the NLCC-GTH PPs when combined with medium-size TZV2P molecularly optimized (MOLOPT) basis sets which are typically utilized in condensed phase simulations, and show that they lead to consistently good results when compared to all-electron calculations for atomization energies, ionization potentials, barrier heights, and non-covalent interactions, but lead to somewhat larger errors for electron affinities.

A Perspective on Sustainable Computational Chemistry Software Development and Integration.

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

When Is a Bond Broken? The Polarizability Perspective.

The question of when a chemical bond can be said to be broken is of fundamental chemical interest but has not been widely studied. Herein we propose that the maxima of static polarizability along bond dissociation coordinates naturally define cutoff points for bond rupture, as they represent the onset of localization of shared electron density into constituent fragments. Examples of computed polarizability maxima over the course of bond cleavage in main-group and transition metal compounds are provided, across covalent, dative and charge-shift bonds. The behavior along reaction paths is also considered. Overall, the static polarizability is found to be a sensitive reporter of electronic structure reorganization associated with bond stretching, and thus can serve as a metric for describing bond cleavage (or diagnose the absence of a chemical bond).

Indolo[2,3-b]quinoxaline as a Low Reduction Potential and High Stability Anolyte Scaffold for Nonaqueous Redox Flow Batteries.

Redox flow batteries (RFBs) are a promising stationary energy storage technology for leveling power supply from intermittent renewable energy sources with demand. A central objective for the development of practical, scalable RFBs is to identify affordable and high-performance redox-active molecules as storage materials. Herein, we report the design, synthesis, and evaluation of a new organic scaffold, indolo[2,3-b]quinoxaline, for highly stable, low-reduction potential, and high-solubility anolytes for nonaqueous redox flow batteries (NARFBs). The mixture of 2- and 3-(tert-butyl)-6-(2-methoxyethyl)-6H-indolo[2,3-b]quinoxaline exhibits a low reduction potential (-2.01 V vs Fc/Fc+), high solubility (>2.7 M in acetonitrile), and remarkable stability (99.86% capacity retention over 49.5 h (202 cycles) of H-cell cycling). This anolyte was paired with N-(2-(2-methoxyethoxy)-ethyl)phenothiazine (MEEPT) to achieve a 2.3 V all-organic NARFB exhibiting 95.8% capacity retention over 75.1 h (120 cycles) of cycling.

Even Faster Exact Exchange for Solids via Tensor Hypercontraction.

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (Nk) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, NISDF) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of NISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. We perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.