Carbon dioxide is a desired feedstock for platform molecules, such as carbon monoxide or higher hydrocarbons, from which we will be able to make many different useful, value-added chemicals. Its catalytic hydrogenation over abundant metals requires the amalgamation of theoretical knowledge with materials design. Here we leverage a theoretical understanding of structure sensitivity, along with a library of different supports, to tune the selectivity of methanation in the Power-to-Gas concept over nickel. For example, we show that carbon dioxide hydrogenation over nickel can and does form propane, and that activity and selectivity can be tuned by supporting different nickel particle sizes on various oxides. This theoretical and experimental toolbox is not only useful for the highly selective production of methane, but also provides new insights for carbon dioxide activation and subsequent carbon–carbon coupling towards value-added products thereby reducing the deleterious effects of this environmentally harmful molecule.

Adaptive quantum mechanical (QM) / molecular mechanical (MM) methods enable efficient molecular simulations of chemistry in solution by describing reactive subregions with an accurate many-body potential energy expression (QM) while the rest of the system is described in a more approximate manner (MM). As solvent molecules diffuse in and out of the reactive region, they are gradually included into (and excluded from) the many-body QM potential. It would be desirable to model such system 1 using an adaptive Hamiltonian, but so far it has resulted in distorted structures at the boundary between the two regions. Here, we propose a Hamiltonian scheme to describe adaptively solvent diffusion across a multi-scale boundary separating configurational potentials that cannot be expressed by a multi-body expansion. The adaptive expressions are entirely general, and complimentary to all standard (non-adaptive) QM/MM embedding schemes available. We demonstrate the validity of our approach on a system described by two different MM potentials (MM/MM'), in which long-range interactions are treated by many-body Ewald summation. Our Hamiltonian approach provides both energy conservation and the correct solvent structure everywhere in the system, thus enabling microcanonical adaptive QM/MM simulations that can be used to obtain vibrational spectra and dynamical properties.

Nucleophilic addition onto a carbonyl moiety is strongly affected by solvent, and correctly simulating this solvent effect is often beyond the capability of single-scale quantum mechanical (QM) models. This work explores multiscale approaches for the description of the reversible and highly solvent-sensitive nucleophilic N|···C=O bond formation in an Me2N–(CH2)3–CH=O molecule. In the first stage of this work, we rigorously compare and test four recent quantum mechanical/molecular mechanical (QM/MM) explicit solvation models, employing a QM description of water molecules in spherical regions around both the oxygen and the nitrogen atom of the solute. The accuracy of the models is benchmarked against a reference QM simulation, focusing on properties of the solvated Me2N–(CH2)3–CH=O molecule in its ring-closed form. In the second stage, we select one of the models (continuous adaptive QM/MM) and use it to obtain a reliable free energy profile for the N|···C bond formation reaction. We find that the dual-sphere approach allows the model to accurately account for solvent reorganization along the entire reaction path. In contrast, a simple microsolvation model cannot adapt to the changing conditions and provides an incorrect description of the reaction process.

Recently, an implementation of the specific reaction parameter (SRP) approach to density functional theory (DFT) was used to study several reactive scattering experiments of H2 on Cu(111). It was possible to obtain chemical accuracy (1 kcal/mol ≈ 4.2 kJ/mol), and therefore, accurately model this paradigmatic example of activated H2 dissociation on a metal surface. In this work, the SRP-DFT methodology is applied to the dissociation of hydrogen on a Pd(111) surface, in order to test whether the SRP-DFT approach is also applicable to non-activated H2-metal systems. In the calculations, the Born-Oppenheimer static surface approximations are used. A comparison to molecular beam sticking experiments, performed at incidence energies ≥125 meV, on H2 + Pd(111) suggested the PBE-vdW [where the Perdew, Burke, and Ernzerhof (PBE) correlation is replaced by van der Waals correlation] functional as a candidate SRP density functional describing the reactive scattering of H2 on Pd(111). Unfortunately, quantum dynamics calculations are not able to reproduce the molecular beam sticking results for incidence energies <125 meV. From a comparison to initial state-resolved (degeneracy averaged) sticking probabilities it seems clear that for H2 + Pd(111) dynamic trapping and steering effects are important, and that these effects are not yet well modeled with the potential energy surfaces considered here. Applying the SRP-DFT method to systems where H2 dissociation is non-activated remains difficult. It is suggested that a density functional that yields a broader barrier distribution and has more non-activated pathways than PBE-vdW (i.e., non-activated dissociation at some sites but similarly high barriers at the high energy end of the spectrum) should allow a more accurate description of the available experiments. Finally, it is suggested that new and better characterized molecular beam sticking experiments be done on H2 + Pd(111), to facilitate the development of a more accurate theoretical description of this system.

No abstract