From CO2 to Methanol by Hybrid QM/MM Embedding

Electron Regulation of Single Indium Atoms at the Active Oxygen Vacancy of In2 O3 (110) for Production of Acetic Acid and Acetone through Direct Coupling of CH4 with CO2

The production of acetic acid and acetone from the direct coupling of CO₂ and CH₄ on the doped In₂ O₃ (110) surface has been studied by extensive first-principles calculations, and the Ga or Al substitution for the single In atom at the active oxygen vacancy of In₂ O₃ (110) can stabilize the reaction species and reduce the free energy barrier of the rate-limiting C-H activation for the conversion of CO₂ and CH₄ to acetic acid. Herein, the metal doping lowers the energy level of partially empty s and p orbitals of In₁ at the oxygen vacancy site and manipulates its electronic properties, resulting in the activity improvement. The stable intermediate with the newly-formed CH₃ COO* has the available In₁ site for subsequent CH₄ activation, which may initiate the direct C-C coupling of CH₃ COO* and CH₃ * to yield C3 species on the doped In₂ O₃ (110). These findings suggest that the metal doping of the active oxygen vacancy opens an avenue for the carbon-chain growth through heterogeneously catalytic coupling of CO₂ and CH₄ .

A highly selective and stable ZnO-ZrO2 solid solution catalyst for CO2 hydrogenation to methanol

Although methanol synthesis via CO hydrogenation has been industrialized, CO₂ hydrogenation to methanol still confronts great obstacles of low methanol selectivity and poor stability, particularly for supported metal catalysts under industrial conditions. We report a binary metal oxide, ZnO-ZrO₂ solid solution catalyst, which can achieve methanol selectivity of up to 86 to 91% with CO₂ single-pass conversion of more than 10% under reaction conditions of 5.0 MPa, 24,000 ml/(g hour), H₂/CO₂ = 3:1 to 4:1, 320° to 315°C. Experimental and theoretical results indicate that the synergetic effect between Zn and Zr sites results in the excellent performance. The ZnO-ZrO₂ solid solution catalyst shows high stability for at least 500 hours on stream and is also resistant to sintering at higher temperatures. Moreover, no deactivation is observed in the presence of 50 ppm SO₂ or H₂S in the reaction stream.

CO2 hydrogenation to methanol on supported Au catalysts under moderate reaction conditions: support and particle size effects

The potential of metal oxide supported Au catalysts for the formation of methanol from CO2 and H2 under conditions favorable for decentralized and local conversion, which could be concepts for chemical energy storage, was investigated. Significant differences in the catalytic activity and selectivity of Au/Al2 O3 , Au/TiO2 , AuZnO, and Au/ZrO2 catalysts for methanol formation under moderate reaction conditions at a pressure of 5 bar and temperatures between 220 and 240 °C demonstrate pronounced support effects. A high selectivity (>50 %) for methanol formation was obtained only for Au/ZnO. Furthermore, measurements on Au/ZnO samples with different Au particle sizes reveal distinct Au particle size effects: although the activity increases strongly with the decreasing particle size, the selectivity decreases. The consequences of these findings for the reaction mechanism and for the potential of Au/ZnO catalysts for chemical energy storage and a "green" methanol technology are discussed.

Point defects in ZnO

We have investigated intrinsic point defects in ZnO and extended this study to Li, Cu and Al impurity centres. Atomic and electronic structures as well as defect energies have been obtained for the main oxidation states of all defects using our embedded cluster hybrid quantum mechanical/molecular mechanical approach to the treatment of localised states in ionic solids. With these calculations we were able to explain the nature of a number of experimentally observed phenomena. We show that in zinc excess materials the energetics of zinc interstitial are very similar to those for oxygen vacancy formation. Our results also suggest assignments for a number of bands observed in photoluminescence and other spectroscopic studies of the material.

Adsorption of organosilanes at a Zn-terminated ZnO (0001) surface: molecular dynamics study

Four different organosilanes (octyltrihydroxysilane, butyltrihydroxysilane, aminopropyltrihydroxysilane, and thiolpropyltrihydroxysilane) adsorbed at a reconstructed Zn-terminated polar ZnO (0001) surface are studied via constant temperature (298 K) molecular dynamics simulations. Both single adsorbed silane molecules as well as adsorbed silane layers are modeled, and the energy, distance, orientation, and alignment of these adsorbates are analyzed. The adsorbed silane molecules exhibit behavior depending on the chemical nature of their tail (nonpolar or polar) as well as on the silane concentration at the solid surface (single adsorption or silane layer). In contrast to the O-terminated ZnO surface studied previously, now adsorption can only occur at the vacancies of this reconstructed crystal surface, thus leading to an arched structure of the liquid phase near the crystal surface. Nevertheless, both nonpolar and polar single adsorbed silanes show a similar orientation and alignment at the surface (orthogonal in the former, parallel in the latter case) as for the O-terminated ZnO surface, although the interaction energy with the surface is considerably increased for nonpolar silanes while it is nearly unaffected for the polar ones. For adsorbed silanes within silane layers, the difference to single adsorbed silanes depends on the polarity of the tail: nonpolar silanes again show an orthogonal alignment, while polar silanes exhibit two different orientations at the solid surface-a head and a tail down configuration. This leads to two completely different but nevertheless stable orientations of these silanes at the Zn-terminated ZnO surface.

Adsorption dynamics of CO2 on Zn-ZnO(0001): a molecular beam study

Presented are initial S(0) and coverage Theta dependent, S(Theta), adsorption probability measurements, respectively, of CO(2) adsorption on the polar Zn-terminated surface of ZnO, parametric in the impact energy E(i), the surface temperature T(s), the impact angle alpha(i), varied along the [001] azimuth, the CO(2) flux, and the density of defects, chi(Ar(+)), as varied by rare gas ion sputtering. S(0) decreases linearly from 0.72 to 0.25 within E(i)=0.12-1.33 eV and is independent of T(s). Above E(i)=0.56 eV, S(0) decreases by approximately 0.2 with increasing alpha(i). The shape of S(Theta) curves is consistent with precursor-mediated adsorption (Kisliuk shape, i.e., S approximately const) for low E(i); above E(i)=0.56 eV, however, a turnover to adsorbate-assisted adsorption (S increases with Theta) has been observed. The initial slope of S(Theta) curves decreases thereby with increasing alpha(i), chi(Ar(+)), and T(s), i.e., the adsorbate-assisted adsorption is most distinct for normal impact on the pristine surface at low T(s) and is independent of the CO(2) flux. The S(Theta) curves have been parametrized by analytic precursor models and Monte Carlo simulations have been conducted as well. The temperature dependence of the saturation coverage shows two structures which could be assigned to adsorption on pristine and intrinsic defect sites, respectively, in agreement with a prior thermal desorption spectroscopy study. The heat of adsorption E(d) for the pristine sites amounts to 34.0-5.4Theta, whereas for adsorption on the intrinsic defect sites E(d) of approximately 43.6 kJ/mol could be estimated. Thus, a kinetic structure-activity relationship was present.

Adsorption dynamics of CO2 on copper-precovered ZnO(0001)–Zn: A molecular-beam scattering and thermal-desorption spectroscopy study

Initial, S(0), as well as coverage-dependent adsorption probability measurements, S(Theta), have been conducted at normal impact angle and as a function of the impact energy of CO(2), E(i), adsorption temperature, T(s), and copper precoverage, Theta(Cu) (at 300 K). S(0), which decreased from approximately 0.4 exponentially to approximately 0.05 with E(i) was independent of Theta(Cu). Astonishingly, S(0) for Cu on ZnO(0001)-Zn is smaller than for the clean support which indicates a chemical modification of the support by the Cu deposits. S(Theta) curves consist of two regimes, a Kisliuk-type and Langmuirian-type section. The first is consistent with capture zone models; the second may indicate direct adsorption of CO(2) on the Cu cluster. The thermal-desorption (TDS) curves for Cu on ZnO(0001)-Zn consist of two structures with binding energies of 26 and approximately 40 kJmol (nu=1 x 10(13) ls). The TDS results indicate that CO(2) populates predominantly the Cu deposits and the rim along the Cu nanoparticles. No indications for CO(2) dissociation could be obtained with Auger electron spectroscopy.

The nature of the oxidation states of gold on ZnO

The interaction between gold in the 0, i, ii and iii oxidation states and the zinc-terminated ZnO(0001) surface is studied via the QM/MM electronic embedding method using density functional theory. The surface sites considered are the vacant zinc interstitial surface site (VZISS) and the bulk-terminated island site (BTIS). We find that on the VZISS, only Au(0) and Au(i) are stable oxidation states. However, all clusters of i to iii oxidation states are stable as substitutionals for Zn2+ in the bulk terminated island site. Au(OH)(x) complexes (x= 1-3) can adsorb exothermically onto the VZISS, indicating that higher oxidation states of gold can be stabilised at this site in the presence of hydroxyl groups. CO is used as a probe molecule to study the reactivity of Au in different oxidation states in VZISS and BTIS. In all cases, we find that the strongest binding of CO is to surface Au(i). Furthermore, CO binding onto Au(0) is stronger when the gold atom is adsorbed onto the VZISS compared to CO binding onto a gas phase neutral gold atom. These results indicate that the nature of the oxidation states of Au on ZnO(0001) will depend on the type of adsorption site. The role of ZnO in Au/ZnO catalysts is not, therefore, merely to disperse gold atoms/particles, but to also modify their electronic properties.

A combined molecular dynamics+quantum mechanics method for investigation of dynamic effects on local surface structures

A combined molecular dynamics (MD)+quantum mechanics (QM) method for studying processes on ionic surfaces is presented. Through the combination of classical MD and ab initio embedded-cluster calculations, this method allows the modeling of surface processes involving both the structural and dynamic features of the substrate, even for large-scale systems. The embedding approach used to link the information from the MD simulation to the cluster calculation is presented, and rigorous tests have been carried out to ensure the feasibility of the method. The electrostatic potential and electron density resulting from our embedded-cluster model have been compared with periodic slab results, and confirm the satisfying quality of our embedding scheme as well as the importance of applying embedding in our combined MD+QM approach. We show that a highly accurate representation of the Madelung potential becomes a prerequisite when the embedded-cluster approach is applied to temperature-distorted surface snapshots from the MD simulation.