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.

Development of a Week-Long Mathematics Intervention for Incoming Chemistry Graduate Students

A student-led mathematics bootcamp has been designed and implemented to help foster community building, improve confidence in mathematical skills, and provide mathematical resources for incoming physical chemistry doctoral students. The bootcamp is held immediately before the start of the first semester of graduate school and uses an active learning approach to review and practice undergraduate-level mathematics problems over 5 days in small student groups. This work includes the development and presentation of a new, publicly available mathematics curriculum for the bootcamp on select mathematics topics, including calculus, linear algebra, functions, differential equations, statistics, and coding in Python, aiming at improving students' confidence and learning experiences in graduate quantum mechanics and statistical physics courses. Surveys before and after the bootcamp showed an increase in students' confidence in problem-solving in key mathematical areas and social aspects of peer-led group learning. Qualitative and quantitative analyses demonstrate that the bootcamp reduced prior inequities in students' confidence metrics based on gender and mathematical background.

The Good, the Bad, and the Ugly: Pseudopotential Inconsistency Errors in Molecular Applications of Density Functional Theory

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.

Computational Studies of Rubber Ozonation Explain the Effectiveness of 6PPD as an Antidegradant and the Mechanism of Its Quinone Formation

The discovery that the commercial rubber antidegradant 6PPD reacts with ozone (O₃) to produce a highly toxic quinone (6PPDQ) spurred a significant research effort into nontoxic alternatives. This work has been hampered by lack of a detailed understanding of the mechanism of protection that 6PPD affords rubber compounds against ozone. Herein, we report high-level density functional theory studies into early steps of rubber and PPD (p-phenylenediamine) ozonation, identifying key steps that contribute to the antiozonant activity of PPDs. In this, we establish that our density functional theory approach can achieve chemical accuracy for many ozonation reactions, which are notoriously difficult to model. Using adiabatic energy decomposition analysis, we examine and dispel the notion that one-electron charge transfer initiates ozonation in these systems, as is sometimes argued. Instead, we find direct interaction between O₃ and the PPD aromatic ring is kinetically accessible and that this motif is more significant than interactions with PPD nitrogens. The former pathway results in a hydroxylated PPD intermediate, which reacts further with O₃ to afford 6PPD hydroquinone and, ultimately, 6PPDQ. This mechanism directly links the toxicity of 6PPDQ to the antiozonant function of 6PPD. These results have significant implications for development of alternative antiozonants, which are discussed.

Copper(III) Metallacyclopentadienes via Zirconocene Transfer and Reductive Elimination to an Isolable Phenanthrocyclobutadiene

Despite the widespread use of copper catalysis for the formation of C-C bonds, debate about the mechanism persists. Reductive elimination from Cu(III) is often invoked as a key step, yet examples of its direct observation from isolable complexes remain limited to only a few examples. Here, we demonstrate that incorporation of bulky mesityl (Mes) groups into the α-positions of a phenanthrene-appended zirconacyclopentadiene, Cp₂Zr(2,5-Mes₂-phenanthro[9,10]C₄), enables efficient oxidative transmetalation to the corresponding, formal Cu(III) metallacyclopentadiene dimer. The dimer was quantitatively converted to a structurally analogous anionic monomer [ⁿBu₄N]{Cl₂Cu(2,5-Mes₂-phenanthro[9,10]C₄)} upon treatment with [ⁿBu₄N][Cl]. Both metallacycles undergo quantitative reductive elimination upon heating to generate phenanthrocyclobutadiene and a Cu(I) species. Due to the steric protection provided by the mesityl groups, this cyclobutadiene was isolated and thoroughly characterized to reveal antiaromaticity comparable to that of free cyclobutadiene, which imbues it with a small highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap of 1.85 eV and accessible reduced and oxidized electronic states.

Optimized Pseudopotentials and Basis Sets for Semiempirical Density Functional Theory for Electrocatalysis Applications

A complete set of Goedecker, Teter, and Hutter (GTH) norm-conserving pseudopotentials (PPs) have been optimized, in conjunction with molecular optimized (MOLOPT) basis sets, for both the B97M-rV and ωB97X-V density functionals for members of the main-group elements and 3d and noble metals. The resulting small-core PPs and corresponding DZVP, TZVP, and TZV2P basis sets yield notable improvements compared to the original PBE defaults when validated against all electron calculations for redox reaction energies and geometries, binding energies, and vibrational Stark effects for metal monocarbonyls in vacuum. Further validation of the optimized PP/MOLOPT basis set combinations was performed using ab initio molecular dynamics simulations and shows greatly improved agreement with experimental trends for metal surface relaxations and the adsorption behavior of CO on solid metal surfaces.

Critical Role of Thermal Fluctuations for CO Binding on Electrocatalytic Metal Surfaces

This work considers the evaluation of density functional theory (DFT) when comparing against experimental observations of CO binding trends on the strong binding Pt(111) and intermediate binding Cu(111) and for weak binding Ag(111) and Au(111) surfaces important in electrocatalysis. By introducing thermal fluctuations using appropriate statistical mechanical NVT and NPT ensembles, we find that the RPBE and B97M-rV DFT functionals yield qualitatively better metal surface strain trends and CO enthalpies of binding for Cu(111) and Pt(111) than found at 0 K, thereby correcting the overbinding by 0.2 to 0.3 eV to yield better agreement with the enthalpies determined from experiment. The importance of dispersion effects are manifest for the weak CO binding Ag(111) and Au(111) surfaces at finite temperatures in which the RPBE functional does not bind CO at all, while the B97M-rV functional shows that the CO-metal interactions are a mixture of chemisorbed and physisorbed species with binding enthalpies that are within ∼0.05 eV of experiment. Across all M(111) surfaces, we show that the B97M-rV functional consistently predicts the correct atop site preference for all metals due to thermally induced surface distortions that preferentially favor the undercoordinated site. This study demonstrates the need to fully account for finite temperature fluctuations to make contact with the binding enthalpies from surface science experiments and electrocatalysis applications.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals

We assess four DFT functionals, RTPSS, RPBE, SCAN and B97M-rV, for surface interactions. We find that B97M-rV predicts the correct site preference for CO binding on Ag and Au while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.

Electronic structure calculations permit identification of the driving forces behind frequency shifts in transition metal monocarbonyls

Our direct DFT decomposition of CO frequency shifts updates the paradigm for metal carbonyl binding.