Bridging the gap between insightful simplicity and successful complexity: From fundamental studies on model systems to technical catalysts

Letter to the Editor

A paper recently published entitled "Water crisis in Iran: A system dynamics approach on water, energy, food, land, and climate (WEFLC) nexus" (Barati et al., 2023). In the mentioned study, a WEFLC model is developed to analyze the water scarcity in Iran. Water crisis, as a complex and challenging issue, has different interdependencies in the context of socio-ecological systems (SES), making it an incorrigible issue. The original paper attempted to assess the water resource dynamics through a systemic lens and explore the impact of various driving forces of water resource planning and management on the water crisis. Iran is a well-studied country, especially around water-related problems. Many interesting facts and findings through the water scarcity analysis in the context of WEFLC are mentioned in the original paper. For instance, it is highlighted that "Mitigation and adaptation policies must be system-oriented and coherent at sectors." However, the original paper did not benefit enough from the previous studies and the full potential of available data. Moreover, some arguments contradict previous findings and, in some cases, are logically flawed. The original paper barely alludes to the nonlinear functional relationships among the components of WEFLC, the core expected component in complex system analysis. Incorrect problem statement formation, flawed methodology, insufficient information on the applied method, ambiguity in models' coupling or cohesion, lack of rational explanation, and inappropriate interpretations of abnormal findings may even mislead many readers. This paper aims to point out some concerns related to the problems mentioned above in the published study, with suggestions to improve the current study and methodological notes for future research.

Tuning the shape and crystal phase of TiO2 nanoparticles for catalysis

Synthesis of TiO₂ nanoparticles with tunable shape and crystal phase has attracted considerable attention for the design of highly efficient heterogeneous catalysts. Tailoring the shape of TiO₂, in the crystal phases of anatase, rutile, brookite and TiO₂(B), allows tuning of the atomic configurations on the dominantly exposed facets for maximizing the active sites and regulating the reaction route towards a specific channel for achieving high selectivity. Moreover, the shape and crystal phase of TiO₂ nanoparticles alter their interactions with metal species, which are commonly termed as strong metal-support interactions involving interfacial strain and charge transfer. On the other hand, metal particles, clusters and single atoms interact differently with TiO₂, because of the variation of the electronic structure, while the surface of TiO₂ determines the interfacial bonding via a geometric effect. The dynamic behavior of the metal-titania interfaces, driven by the chemisorption of the reactive molecules at elevated temperatures, also plays a decisive role in elaborating the structure-reactivity relationship.

In situandoperandoelectron microscopy in heterogeneous catalysis-insights into multi-scale chemical dynamics

This review features state-of-the-artin situandoperandoelectron microscopy (EM) studies of heterogeneous catalysts in gas and liquid environments during reaction. Heterogeneous catalysts are important materials for the efficient production of chemicals/fuels on an industrial scale and for energy conversion applications. They also play a central role in various emerging technologies that are needed to ensure a sustainable future for our society. Currently, the rational design of catalysts has largely been hampered by our lack of insight into the working structures that exist during reaction and their associated properties. However, elucidating the working state of catalysts is not trivial, because catalysts are metastable functional materials that adapt dynamically to a specific reaction condition. The structural or morphological alterations induced by chemical reactions can also vary locally. A complete description of their morphologies requires that the microscopic studies undertaken span several length scales. EMs, especially transmission electron microscopes, are powerful tools for studying the structure of catalysts at the nanoscale because of their high spatial resolution, relatively high temporal resolution, and complementary capabilities for chemical analysis. Furthermore, recent advances have enabled the direct observation of catalysts under realistic environmental conditions using specialized reaction cells. Here, we will critically discuss the importance of spatially-resolvedoperandomeasurements and the available experimental setups that enable (1) correlated studies where EM observations are complemented by separate measurements of reaction kinetics or spectroscopic analysis of chemical species during reaction or (2) real-time studies where the dynamics of catalysts are followed with EM and the catalytic performance is extracted directly from the reaction cell that is within the EM column or chamber. Examples of current research in this field will be presented. Challenges in the experimental application of these techniques and our perspectives on the field's future directions will also be discussed.

Surface-Dependent Activation of Model α-Al2 O3 -Supported P-Doped Hydrotreating Catalysts Prepared by Spin Coating

Requirements for improved catalytic formulations is continuously driving research in hydrotreating (HDT) catalysis for biomass upgrading and heteroatom removal for cleaner fuels. The present work proposes a surface-science approach for the understanding of the genesis of the active (sulfide) phase in model P-doped MoS₂ hydrotreating catalysts supported on α-Al₂ O₃ single crystals. This approach allows one to obtain a surface-dependent insight by varying the crystal orientations of the support. Model phosphorus-doped catalysts are prepared via spin-coating of Mo-P precursor solutions onto four α-Al₂ O₃ crystal orientations, C(0001), A(11 2 ‾ 0), M(10 1 ‾ 0) and R(1 1 ‾ 02) that exhibit different speciations of surface -OH. ³¹ P and ⁹⁵ Mo liquid-state NMR are used to give a comprehensive description of the Mo and P speciation of the phospho-molybdic precursor solution. The speciation of the deposition solution is then correlated with the genesis of the active MoS₂ phase. XPS quantification of the surface P/Mo ratio reveal a surface-dependent phosphate aggregation driven by the amount of free phosphates in solution. Phosphates aggregation decreases in the following order C(0001)≫M(10 1 ‾ 0)>A(11 2 ‾ 0), R(1 1 ‾ 02). This evolution can be rationalized by an increasing strength of phosphate/surface interactions on the different α-Al₂ O₃ surface orientations from the C(0001) to the R(1 1 ‾ 02) plane. Retardation of the sulfidation with temperature is observed for model catalysts with the highest phosphate dispersion on the surface (A(11 2 ‾ 0), R(1 1 ‾ 02)), suggesting that phosphorus strongly intervene in the genesis of the active phase through a close intimacy between phosphates and molybdates. The surface P/Mo ratio appears as a key descriptor to quantify this retarding effect. It is proposed that retardation of sulfidation is driven by two effects: i) a chemical inhibition through formation of hardly reducible mixed molybdo-phosphate structures and ii) a physical inhibition with phosphate clusters inhibiting the growth of MoS₂ . The surface-dependent phosphorus doping on model α-Al₂ O₃ supports can be used as a guide for the rational design of more efficient HDT catalysts on industrial γ-Al₂ O₃ carrier.

Recycling of CO2 by Hydrogenation of Carbonate Derivatives to Methanol: Tuning Copper-Oxide Promotion Effects in Supported Catalysts

The selective hydrogenation of organic carbonates to methanol is a relevant transformation to realize flexible processes for the recycling of waste CO₂ with renewable H₂ mediated by condensed carbon dioxide surrogates. Oxide-supported copper nanoparticles are promising solid catalysts for this selective hydrogenation. However, essential for their optimization is to rationalize the prominent impact of the oxide support on performance. Herein, the role of Lewis acid centers, exposed on the oxide support at the periphery of the Cu nanoparticles, was systematically assessed. For the hydrogenation of propylene carbonate, as a model cyclic carbonate, the conversion rate, the apparent activation energy, and the selectivity to methanol correlate with the Lewis acidity of the coordinatively unsaturated cationic sites exposed on the oxide support. Lewis sites of markedly low and high electron-withdrawing character promote unselective decarbonylation and decarboxylation reaction pathways, respectively. Supports exposing Lewis sites of intermediate acidity maximize the selectivity to methanol while inhibiting acid-catalyzed secondary reactions of the propanediol product, and thus enable its recovery in cyclic processes of CO₂ hydrogenation mediated by condensed carbonate derivatives. These findings help rationalize metal-support promotion effects that determine the performance of supported metal nanoparticles in this and other selective hydrogenation reactions of significance in the context of sustainable chemistry.

Catalytic CO oxidation on Pt under near ambient pressure: A NAP-LEEM study

A near ambient pressure low-energy electron microscope (NAP-LEEM) has recently been constructed, that allows in situ imaging of surfaces up to a pressure of 10^-1 mbar. Here we report on pattern formation in catalytic CO oxidation on a Pt(110) single crystal surface and on a polycrystalline Pt foil in the 10^-2 mbar range, operating the microscope in the mirror electron microscopy (MEM) and in the LEEM mode. Excitations localized at structural defects and spiral wave fragments have been observed.

Chemistry and Quantum Mechanics in 2019: Give Us Insight and Numbers

This Perspective revisits Charles Coulson’s famous statement from 1959 “give us insight not numbers” in which he pointed out that accurate computations and chemical understanding often do not go hand in hand. We argue that today, accurate wave function based first-principle calculations can be performed on large molecular systems, while tools are available to interpret the results of these calculations in chemical language. This leads us to modify Coulson’s statement to “give us insight and numbers”. Examples from organic, inorganic, organometallic and surface chemistry as well as molecular magnetism illustrate the points made.

Stand-alone polarization-modulation infrared reflection absorption spectroscopy instrument optimized for the study of catalytic processes at elevated pressures

This paper describes the design and construction of a compact, "user-friendly" polarization-modulation infrared reflection absorption spectroscopy (PM-IRRAS) instrument at the Center for Functional Nanomaterials (CFN) of Brookhaven National Laboratory, which allows studying surfaces at pressures ranging from ultra-high vacuum to 100 Torr. Surface infrared spectroscopy is ideally suited for studying these processes as the vibrational frequencies of the IR chromophores are sensitive to the nature of the bonding environment on the surface. Relying on the surface selection rules, by modulating the polarization of incident light, it is possible to separate the contributions from the isotropic gas or solution phase, from the surface bound species. A spectral frequency range between 1000 cm^-1 and 4000 cm^-1 can be acquired. While typical spectra with a good signal to noise ratio can be obtained at elevated pressures of gases in ∼2 min at 4 cm^-1 resolution, we have also acquired higher resolution spectra at 0.25 cm^-1 with longer acquisition times. By way of verification, CO uptake on a heavily oxidized Ru(0001) sample was studied. As part of this test study, the presence of CO adsorbed on Ru bridge sites was confirmed, in agreement with previous ambient pressure X ray photoelectron spectroscopy studies. In terms of instrument performance, it was also determined that the gas phase contribution from CO could be completely removed even up to pressures close to 100 Torr. A second test study demonstrated the use of the technique for studying morphological properties of a spin coated polymer on a conductive surface. Note that this is a novel application of this technique. In this experiment, the polarization of incident light was modulated manually (vs. through a photoelastic modulator). It was demonstrated, in good agreement with the literature, that the polymer chains preferentially lie parallel with the surface. This PM-IRRAS system is small, modular, and easily reconfigurable. It also features a "vacuum suitcase" that allows for the integration of the PM-IRRAS system with the rest of the suite of instrumentation at our laboratory available to external users through the CFN user proposal system.

Surface Adsorption Energetics Studied with "Gold Standard" Wave-Function-Based Ab Initio Methods: Small-Molecule Binding to TiO2(110)

Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T)) is widely considered to be the "gold standard" of ab initio quantum chemistry. Using the domain-based pair natural orbital local correlation concept (DLPNO-CCSD(T)), these calculations can be performed on systems with hundreds of atoms at an accuracy of ∼99.9% of the canonical CCSD(T) method. This allows for ab initio calculations providing reference adsorption energetics at solid surfaces with an accuracy approaching 1 kcal/mol. This is an invaluable asset, not least for the assessment of density functional theory (DFT) as the prevalent approach for large-scale production calculations in energy or catalysis applications. Here we use DLPNO-CCSD(T) with embedded cluster models to compute entire adsorbate potential energy surfaces for the binding of a set of prototypical closed-shell molecules (H₂O, NH₃, CH₄, CH₃OH, CO₂) to the rutile TiO₂(110) surface. The DLPNO-CCSD(T) calculations show excellent agreement with available experimental data, even for the "infamous" challenge of correctly predicting the CO₂ adsorption geometry. The numerical efficiency of the approach is within 1 order of magnitude of hybrid-level DFT calculations, hence blurring the borders between reference and production technique.