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 (O3) 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 O3 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 O3 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, Cp2Zr(2,5-Mes2-phenanthro[9,10]C4), enables efficient oxidative transmetalation to the corresponding, formal Cu(III) metallacyclopentadiene dimer. The dimer was quantitatively converted to a structurally analogous anionic monomer [nBu4N]{Cl2Cu(2,5-Mes2-phenanthro[9,10]C4)} upon treatment with [nBu4N][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.

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.

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.

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

Density Functional Theory (DFT) is currently the most tractable choice of theoretical model used to understand the mechanistic pathways for electrocatalytic processes such as CO2 or CO reduction. Here, we assess the performance of two DFT functionals designed specifically to describe surface interactions, RTPSS and RPBE, as well as two popular meta-GGA functionals, SCAN and B97M-rV, that have not been a priori optimized for better interfacial properties. We assess all four functionals against available experimental data for prediction of bulk and bare surface properties on four electrocatalytically relevant metals, Au, Ag, Cu, and Pt, and for binding CO to surfaces of these metals. To partially mitigate issues such as thermal and anharmonic corrections associated with comparing computations with experiments, molecular benchmarks against high level quantum chemistry are reported for CO complexes with Au, Ag, Cu and Pt atoms, as well as the CO-water complex and the water dimer. Overall, we find that the surface modified RPBE functional performs reliably for many of the benchmarks examined here, and the meta-GGA functionals also show promising results. Specifically B97M-rV predicts the correct site preference for CO binding on Ag and Au (the only functional tested here to do so), 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.

We report the adiabatic energy decomposition analysis (EDA) of density functional theory (DFT) results, shedding light on the physical content of binding energies and carbon monoxide (CO) frequency (υCO) shifts in select first-row transition metal monocarbonyls (MCOs; M = Ti-, V-, Cr-, Co-, Ni-, Cu-, V, Cr, Mn, Ni, Cu, Zn, Cr+, Mn+, Fe+, Cu+, and Zn+). This approach allows for the direct decomposition of υCO, in contrast to previous studies of these systems. Neutral, anionic, and cationic systems are compared, and our results indicate that the relative importance of electrostatic interactions, intramolecular orbital polarization, and charge transfer can vary significantly with the charge and electron configuration of the metal participating in binding. Various anomalous systems are also discussed and incorporated into a general model of MCO binding. Electrostatic interactions and orbital polarization are found to promote blue shifts in υCO, while charge transfer effects encourage υCO red-shifting; previously reported values of υCO are found to be a result of a complex but quantifiable interplay between these physical components. Our computations indicate that CuCO- and ZnCO possess triplet ground states, and also that CrCO- exhibits a non-linear geometry, all in contrast to previous computational results. Advantages and limitations of this model as an approximation to more complicated systems, like those implicated in heterogeneous catalysis, are discussed. We also report benchmark results for MCO geometries, binding energies, and harmonic CO frequencies, and discuss the validity of single-reference wave function and DFT approaches to the study of these transition metal systems.

Reinterpreting the infrared spectrum of H + HCN: Methylene amidogen radical and its coproducts.

The methylene amidogen radical (H2CN) plays a role in high-energy material combustion and extraterresterial atmospheres. Recent theoretical work has struggled to match experimental assignments for its CN and antisymmetric CH2 stretching frequencies (ν2 and ν5), which were reported to occur at 1725 and 3103 cm-1. Herein, we compute the vibrational energy levels of this molecule by extrapolating quadruples-level coupled-cluster theory to the complete basis limit and adding corrections for vibrational anharmonicity. This level of theory predicts that ν2 and ν5 should occur at 1646 and 2892 cm-1, at odds with the experimental assignments. To investigate the possibility of defects in our theoretical treatment, we analyze the sensitivity of our approach to each of its contributing approximations. Our analysis suggests that the observed deviation from experiment is too large to be explained as an accumulation of errors, leading us to conclude that these transitions were misassigned. To help resolve this discrepancy, we investigate possible byproducts of the H + HCN reaction, which was the source of H2CN in the original experiment. In particular, we predict vibrational spectra for cis-HCNH, trans-HCNH, and H2CNH using high-level coupled-cluster computations. Based on these results, we reassign the transition at 1725 cm-1 to ν3 of trans-HCNH, yielding excellent agreement. Supporting this identification, we assign a known contaminant peak at 886 cm-1 to ν5 of the same conformer. Our computations suggest that the peak observed at 3103 cm-1, however, does not belong to any of the aforementioned species. To facilitate further investigation, we use structure and bonding arguments to narrow the range of possible candidates. These arguments lead us to tentatively put forth formaldazine [(H2CN)2] as a suggestion for further study, which we support with additional computations.