Tailored Platinum Group Metal/Spinel Oxide Catalysts for Dynamically Enhanced Methane Oxidation.

A combined experimental and molecular modeling study identifies a family of spinel oxides that in combination with PGM (platinum group metals) provide enhanced methane oxidation activity. With a reduction in greenhouse gas (GHG) emissions urgently needed, there is renewed interest in the use of natural gas vehicles (NGVs) and engines (NGEs) for transportation, commerce, and industrial applications. NGVs and NGEs emit less CO2 than their petroleum-derived counterparts but may emit uncombusted methane, an even more potent GHG. For stoichiometric engines, methane oxidation catalysts containing PGM and spinel oxide in layered architectures offer increased methane oxidation activity and lower light-off temperatures (T50). The reducible spinel oxide has direct and indirect roles that are effectively described by the bulk oxygen vacancy formation energy (Evac). We apply density functional theory (DFT) to identify several earth-abundant, cobalt-rich spinel oxides with favorable Evac, shown to correlate with dynamic oxygen storage capacity (DOSC) and CO and H2 oxidation activity. We experimentally rank-order the DFT-identified spinel oxides in combination with Pt+Pd for their methane oxidation activity measurements, under both time-invariant and modulated feed conditions. We show good agreement between the activity and the DFT-computed reducibility of the spinel oxide. The findings suggest spinel reducibility is a key factor in achieving enhanced low-temperature methane conversion, enabled through a balance of methane activation on the PGM sites and subsequent oxidation of the intermediates and byproducts on spinel oxides. In agreement with its computationally predicted Evac, NiCo2O4 was confirmed to have the highest DOSC and lowest T50 among the tested spinel samples.

Spooky correlations and unusual van der Waals forces between gapless and near-gapless molecules.

We consider the zero-temperature van der Waals (vdW) interaction between two molecules, each of which has a zero or near-zero electronic gap between a ground state and the first excited state, using a toy model molecule (equilateral H3) as an example. We show that the van der Waals energy between two ground state molecules falls off as D-3 instead of the usual D-6 dependence, when the molecules are separated by distance D. We show that this is caused by a perfect "spooky" correlation between the two fluctuating electric dipoles. The phenomenon is related to, but not the same as, the "resonant" interaction between an electronically excited and a ground state molecule introduced by Eisenschitz and London in 1930. It is also an example of "type C van der Waals non-additivity" recently introduced by one of us [J. F. Dobson, Int. J. Quantum Chem. 114, 1157 (2014)]. Our toy molecule H3 is not stable, but symmetry considerations suggest that a similar vdW phenomenon may be observable, despite Jahn-Teller effects, in molecules with a discrete rotational symmetry and broken inversion symmetry, such as certain metal atom clusters. The motion of the nuclei will need to be included for a definitive analysis of such cases, however.

Investigation of the active sites of rhodium sulfide for hydrogen evolution/oxidation using carbon monoxide as a probe.

Carbon monoxide (CO) was observed to decrease the activity for hydrogen evolution, hydrogen oxidation, and H2-D2 exchange on rhodium sulfide, platinum, and rhodium metal. The temperature at which the CO was desorbed from the catalyst surface (detected by recovery in the H2-D2 exchange activity of the catalyst) was used as a descriptor for the CO binding energy to the active site. The differences in the CO desorption temperature between the different catalysts showed that the rhodium sulfide active site is not metallic rhodium. Using density functional theory, the binding energy of CO to the Rh sites in rhodium sulfide is found comparable to the binding energy on Pt. Coupled with experiment this supports the proposition that rhodium rather than sulfur atoms in the rhodium sulfide are the active site for the hydrogen reaction. This would indicate the active sites for hydrogen evolution/oxidation as well as oxygen reduction (determined by other groups using X-ray absorption spectroscopy) may be the same.

Transition metal sulfide hydrogen evolution catalysts for hydrobromic acid electrolysis.

Mixed metal sulfides containing combinations of W, Fe, Mo, Ni, and Ru were synthesized and screened for activity and stability for the hydrogen evolution reaction (HER) in aqueous hydrobromic acid (HBr). Co- and Ni-substituted RuS(2) were identified as potentially active HER electrocatalysts by high-throughput screening (HTS), and the specific compositions Co(0.4)Ru(0.6)S(2) and Ni(0.6)Ru(0.4)S(2) were identified by optimization. Hydrogen evolution activity of Co(0.4)Ru(0.6)S(2) in HBr is greater than RuS(2) or CoS(2) and comparable to Pt and commercial Rh(x)S(y). Structural and morphological characterizations of the Co-substituted RuS(2) suggest that the nanoparticulate solids are a homogeneous solid solution with a pyrite crystal structure. No phase separation is detected for Co substitutions below 30% by X-ray diffraction. In 0.5 M HBr electrolyte, the Co-Ru electrode material synthesized with 30% Co rapidly lost approximately 34% of the initial loading of Co; thereafter, it was observed to exhibit stable activity for HER with no further loss of Co. Density functional theory calculations indicate that the S(2)(2-) sites are the most important for HER and the presence of Co influences the S(2)(2-) sites such that the hydrogen binding energy at sufficiently high hydrogen coverage is decreased compared to ruthenium sulfide. Although showing high HER activity in a flow cell, the reverse reaction of hydrogen oxidation is slow on the RuS(2) catalysts tested when compared to platinum and rhodium sulfide, leaving rhodium sulfide as the only suitable tested material for a regenerative HBr cell due its stability compared to platinum.

Dispersion interaction in hydrogen-chain models.

We have investigated the dispersion interaction in hydrogen chain models via density functional theory-based symmetry-adapted perturbation theory using the asymptotically corrected PBE0 energy functional. The quasimetallic and the insulating prototype systems were chosen to be hydrogen chains with equally and alternately spaced H(2) units, respectively. The dependence of the dispersion energy on the chain length for quasimetallic and insulating cases has been determined for two chains arranged either in pointing or in parallel geometries. The results are compared with those previously calculated from a continuum coupled-plasmon approach [Phys. Rev. B 77, 075436 (2008)]. The interaction energy has also been modeled by pairwise summations over short fragments of the chains, demonstrating the failure of the additivity principle for the quasimetallic case, while confirming that the additivity is a qualitatively reasonable hypothesis for the insulating case.

Correlation Energy Expressions from the Adiabatic-Connection Fluctuation-Dissipation Theorem Approach.

We explore several random phase approximation (RPA) correlation energy variants within the adiabatic-connection fluctuation-dissipation theorem approach. These variants differ in the way the exchange interactions are treated. One of these variants, named dRPA-II, is original to this work and closely resembles the second-order screened exchange (SOSEX) method. We discuss and clarify the connections among different RPA formulations. We derive the spin-adapted forms of all the variants for closed-shell systems and test them on a few atomic and molecular systems with and without range separation of the electron-electron interaction.

Aurophilic Interactions from Wave Function, Symmetry-Adapted Perturbation Theory, and Rangehybrid Approaches.

The aurophilic interaction is examined in three model systems Au2((3)Σg(+)), (AuH)2, and (HAuPH3)2 which contain interactions of pairs of the Au centers in the oxidation state (I). Several methods are employed ranging from wave function theory-based (WFT) approaches to symmetry-adapted perturbation theory (SAPT) and range-separated hybrid (RSH) density functional theory (DFT) methods. The most promising and accurate approach consists of a combination of the DFT and WFT approaches in the RSH framework. In this combination the short-range DFT handles the slow convergence of the correlation cusp, whereas the long-range WFT is best suited for the long-range correlation. Of the three tested RSH DFT methods, the one which uses a short-range exchange functional based on the Ernzerhof-Perdew exchange hole model with a range-separation parameter of 0.4 bohr(-1) seems to be the best candidate for treatment of gold. In combination with the long-range coupled cluster singles, doubles, and noniterative triples [CCSD(T)] treatment it places the strength of aurophilic bonding in (HAuPH3)2 at 5.7 kcal/mol at R = 3.09 Å. This value is somewhat larger than our best purely WFT result based on CCSD(T), 4.95 kcal/mol (R = 3.1 Å), and considerably smaller than the Hartree-Fock+dispersion value of 7.4 kcal/mol (R = 2.9 Å). The 5.7 kcal/mol estimate fits reasonably well within the prediction of the empirical relationship proposed by Schwerdtfeger et al. (J. Am. Chem. Soc.1998, 120, 6587). A direct computation of dispersion energy, including exchange corrections, results in values of ca. -9 kcal/mol for Au2((3)Σg(+)) and (AuH)2 and -13 kcal/mol for (HAuPH3)2 at the distance of a typical aurophilic bond, R = 3.0 Å.

On the equivalence of ring-coupled cluster and adiabatic connection fluctuation-dissipation theorem random phase approximation correlation energy expressions.

The correlation energy in the direct random phase approximation (dRPA) can be written, among other possibilities, either in terms of the interaction strength averaged correlation density matrix, or in terms of the coupled cluster doubles amplitudes obtained in the direct ring approximation (drCCD). Although the corresponding dRPA correlation density matrix on the one hand, and the drCCD amplitude matrix on the other hand, differ significantly, they yield identical energies. Similarly, the analogous RPA and rCCD correlation energies calculated from antisymmetrized two-electron integrals are identical to each other despite very different underlying working equations. In the present communication, a direct correspondence between amplitudes and densities is established and investigated with perturbation theory arguments. Our analysis also sheds some light on the properties of recently proposed RPA/rCCD variants which use antisymmetrized integrals in part of the equations and nonantisymmetrized integrals in others.