Anharmonic Correction to Free Energy Barriers from DFT-Based Molecular Dynamics Using Constrained Thermodynamic Integration.

For the calculation of anharmonic contributions to free energy barriers, constrained thermodynamic λ-path integration (λ-TI) from a harmonic reference force field to density functional theory is presented as an alternative to the established Blue Moon ensemble method (ξ-TI), in which free energy gradients along the reaction coordinate ξ are integrated. With good agreement in all cases, the λ-TI method is benchmarked against the ξ-TI method for several reactions, including the internal CH3 group rotation in ethane, a nucleophilic substitution of CH3Cl, a retro-Diels-Alder reaction, and a proton transfer in zeolite H-SSZ-13. An advantage of λ-TI is that one can use virtually any reference state to compute anharmonic contributions to reaction free energies or free energy barriers. This is particularly relevant for catalysis, where it is now possible to compute anharmonic corrections to the free energy of a transition state relative to any reference, for example, the most stable state of the active site and the reactants in the gas phase. This is in contrast to ξ-TI, where free energy barriers can only be computed relative to an initial state with all reactants coadsorbed. Finally, the Bennett acceptance ratio method combined with λ-TI is demonstrated to reduce the number of required integration grid points with tolerable accuracy, favoring thus λ-TI over ξ-TI in terms of computational efficiency.

Anharmonic Correction to Adsorption Free Energy from DFT-Based MD Using Thermodynamic Integration.

Adsorption processes are often governed by weak interactions for which the estimation of entropy contributions by means of the harmonic approximation is prone to be inaccurate. Thermodynamic integration (TI) from the harmonic to the fully interacting system (λ-path integration) can be used to compute anharmonic corrections. Here, we combine TI with (curvilinear) internal coordinates in periodic systems to make the formalism available in computational studies. Our implementation of ab initio molecular dynamics in VASP is independent of the reaction path and can be thus applied to study adsorption processes relative to the gas phase and does hence provide a useful tool for computational catalysis. We discuss the application of the approach on three model systems for which exact semianalytical solutions exist and illustrate and quantify the importance of anharmonic vibrations, hindered rotations, and hindered translations (dissociation). Eventually, we apply the method to study the adsorption of small adsorbates in a zeolite (H-SSZ-13).

Prospects of Heterogeneous Hydroformylation with Supported Single Atom Catalysts.

The potential of oxide-supported rhodium single atom catalysts (SACs) for heterogeneous hydroformylation was investigated both theoretically and experimentally. Using high-level domain-based local-pair natural orbital coupled cluster singles doubles with perturbative triples contribution (DLPNO-CCSD(T)) calculations, both stability and catalytic activity were investigated for Rh single atoms on different oxide surfaces. Atomically dispersed, supported Rh catalysts were synthesized on MgO and CeO2. While the CeO2-supported rhodium catalyst is found to be highly active, this is not the case for MgO, most likely due to increased confinement, as determined by extended X-ray absorption fine structure spectroscopy (EXAFS), that diminishes the reactivity of Rh complexes on MgO. This agrees well with our computational investigation, where we find that rhodium carbonyl hydride complexes on flat oxide surfaces such as CeO2(111) have catalytic activities comparable to those of molecular complexes. For a step edge on a MgO(301) surface, however, calculations show a significantly reduced catalytic activity. At the same time, calculations predict that stronger adsorption at the higher coordinated adsorption site leads to a more stable catalyst. Keeping the balance between stability and activity appears to be the main challenge for oxide supported Rh hydroformylation catalysts. In addition to the chemical bonding between rhodium complex and support, the confinement experienced by the active site plays an important role for the catalytic activity.

One-Pot Cooperation of Single-Atom Rh and Ru Solid Catalysts for a Selective Tandem Olefin Isomerization-Hydrosilylation Process.

Realizing the full potential of oxide-supported single-atom metal catalysts (SACs) is key to successfully bridge the gap between the fields of homogeneous and heterogeneous catalysis. Here we show that the one-pot combination of Ru1 /CeO2 and Rh1 /CeO2 SACs enables a highly selective olefin isomerization-hydrosilylation tandem process, hitherto restricted to molecular catalysts in solution. Individually, monoatomic Ru and Rh sites show a remarkable reaction specificity for olefin double-bond migration and anti-Markovnikov α-olefin hydrosilylation, respectively. First-principles DFT calculations ascribe such selectivity to differences in the binding strength of the olefin substrate to the monoatomic metal centers. The single-pot cooperation of the two SACs allows the production of terminal organosilane compounds with high regio-selectivity (>95 %) even from industrially-relevant complex mixtures of terminal and internal olefins, alongside a straightforward catalyst recycling and reuse. These results demonstrate the significance of oxide-supported single-atom metal catalysts in tandem catalytic reactions, which are central for the intensification of chemical processes.

Thermodynamics and reaction mechanism of urea decomposition.

The selective catalytic reduction technique for automotive applications depends on ammonia production from a urea-water solution via thermolysis and hydrolysis. In this process, undesired liquid and solid by-products are formed in the exhaust pipe. The formation and decomposition of these by-products have been studied by thermogravimetric analysis and differential scanning calorimetry. A new reaction scheme is proposed that emphasizes the role of thermodynamic equilibrium of the reactants in liquid and solid phases. Thermodynamic data for triuret have been refined. The observed phenomenon of liquefaction and re-solidification of biuret in the temperature range of 193-230 °C is explained by formation of a eutectic mixture with urea.