The pseudopotential (PP) approximation is one of the most common techniques in computational chemistry. Despite its long history, the development of custom PPs has not tracked with the explosion of different density functional approximations (DFAs). As a result, the use of PPs with exchange/correlation models for which they were not developed is widespread, although this practice is known to be theoretically unsound. The extent of PP inconsistency errors (PPIEs) associated with this practice has not been systematically explored across the types of energy differences commonly evaluated in chemical applications. We evaluate PPIEs for a number of PPs and DFAs across 196 chemically relevant systems of both transition-metal and main-group elements, as represented by the W4-11, TMC34, and S22 data sets. Near the complete basis set limit, these PPs are found to cleanly approach all-electron (AE) results for noncovalent interactions but introduce root-mean-squared errors (RMSEs) upwards of 15 kcal mol^-1 into predictions of covalent bond energies for a number of popular DFAs. We achieve significant improvements through the use of empirical atom- and DFA-specific PP corrections, indicating considerable systematicity of the PPIEs. The results of this work have implications for chemical modeling in both molecular contexts and for DFA design, which we discuss.

COVID-19 is a devastating respiratory and inflammatory illness caused by a new coronavirus that is rapidly spreading throughout the human population. Over the past 12 months, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19, has already infected over 160 million (>20% located in United States) and killed more than 3.3 million people around the world (>20% deaths in USA). As we face one of the most challenging times in our recent history, there is an urgent need to identify drug candidates that can attack SARS-CoV-2 on multiple fronts. We have therefore initiated a computational dynamics drug pipeline using molecular modeling, structure simulation, docking and machine learning models to predict the inhibitory activity of several million compounds against two essential SARS-CoV-2 viral proteins and their host protein interactors-S/Ace2, Tmprss2, Cathepsins L and K, and Mpro-to prevent binding, membrane fusion and replication of the virus, respectively. All together, we generated an ensemble of structural conformations that increase high-quality docking outcomes to screen over >6 million compounds including all FDA-approved drugs, drugs under clinical trial (>3000) and an additional >30 million selected chemotypes from fragment libraries. Our results yielded an initial set of 350 high-value compounds from both new and FDA-approved compounds that can now be tested experimentally in appropriate biological model systems. We anticipate that our results will initiate screening campaigns and accelerate the discovery of COVID-19 treatments.

Unsupported non-bridged uranium-carbon double bonds have long been sought after in actinide chemistry as fundamental synthetic targets in the study of actinide-ligand multiple bonding. Here we report that, utilizing I_h(7)-C₈₀ fullerenes as nanocontainers, a diuranium carbide cluster, U=C=U, has been encapsulated and stabilized in the form of UCU@I_h(7)-C₈₀. This endohedral fullerene was prepared utilizing the Krätschmer-Huffman arc discharge method, and was then co-crystallized with nickel(II) octaethylporphyrin (Ni^II-OEP) to produce UCU@I_h(7)-C₈₀·[Ni^II-OEP] as single crystals. X-ray diffraction analysis reveals a cage-stabilized, carbide-bridged, bent UCU cluster with unexpectedly short uranium-carbon distances (2.03 Å) indicative of covalent U=C double-bond character. The quantum-chemical results suggest that both U atoms in the UCU unit have formal oxidation state of +5. The structural features of UCU@I_h(7)-C₈₀ and the covalent nature of the U(f¹)=C double bonds were further affirmed through various spectroscopic and theoretical analyses.