Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ground State Properties of the Wide Band Gap Semiconductor Beryllium Sulfide (BeS)

We report the results from self-consistent calculations of electronic, transport, and bulk properties of beryllium sulfide (BeS) in the zinc-blende phase, and employed an ab-initio local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). We obtained the ground state properties of zb-BeS with the Bagayoko, Zhao, and Williams (BZW) computational method, as enhanced by Ekuma and Franklin (BZW-EF). Our findings include the electronic energy bands, the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated band structure clearly shows that zb-BeS has an indirect energy band gap of 5.436 eV, from Γ to a point between Γ and X, for an experimental lattice constant of 4.863 Å. This is in excellent agreement with the experiment, unlike the findings of more than 15 previous density functional theory (DFT) calculations that did not perform the generalized minimization of the energy functional, required by the second DFT theorem, which is inherent to the implementation of our BZW-EF method.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Highly tunable band inversion in AB 2 X 4 (A=Ge, Sn, Pb; B=As, Sb, Bi; X=Se, Te) compounds

Topological materials have been discovered so far largely by searching for existing compounds in crystallographic databases, but there are potentially new topological materials with desirable features that have not been synthesized. One of the desirable features is high tunability resulting from the band inversion with a very small direct band gap, which can be tuned by changes in pressure or strain to induce a topological phase transition. Here, using density-functional theory (DFT) calculations, we have studied the septuple layered AB 2 X 4 series compounds, where A=(Ge, Sn and Pb), B=(As, Sb and Bi), and X=(Se and Te). With the DFT thermodynamic stability validated by the already-reported compounds in these series, we predict stable Se compounds, which are not found in crystallographic database. Among them, we find that GeBi 2 Se 4 and GeSb 2 Se 4 having a small direct band gap at the Z point are very close to a strong topological insulator, which can be tuned by a moderate pressure to induce the band inversion. Importantly, the topological features with the small direct band gap are well isolated in both momentum and energy windows, which offers high tunability for studying the topological phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution

Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding of Active Sites and Interconversion of Pd and PdO during CH 4 Oxidation

Pd-based catalysts are widely used in the oxidation of CH 4 and have a significant impact on global warming. However, understanding their active sites remains controversial, because interconversion between Pd and PdO occurs consecutively during the reaction. Understanding the intrinsic active sites under reaction conditions is critical for developing highly active and selective catalysts. In this study, we demonstrated that partially oxidized palladium (PdO x ) on the surface plays an important role for CH 4 oxidation. Regardless of whether the initial state of Pd corresponds to oxides or metallic clusters, the topmost surface is PdO x , which is formed during CH4 oxidation. A quantitative analysis using CO titration, diffuse reflectance infrared Fourier-transform spectroscopy, X-ray diffraction, and scanning transmission electron microscopy demonstrated that a surface PdO layer was formed on top of the metallic Pd clusters during the CH 4 oxidation reaction. Furthermore, the time-on-stream test of CH 4 oxidation revealed that the presence of the PdO layer on top of the metallic Pd clusters improves the catalytic activity. Our periodic density functional theory (DFT) calculations with a PdO x slab and nanoparticle models aided the elucidation of the structure of the experimental PdO particles, as well as the experimental C-O bands. The DFT results also revealed the formation of a PdO layer on the metallic Pd clusters. This study helps achieve a fundamental understanding of the active sites of Pd and PdO for CH 4 oxidation and provides insights into the development of active and durable Pd-based catalysts through molecular-level design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Temperature-dependent lattice dynamics in polycrystalline ZrN

Zirconium nitride holds promise for improving nuclear reactor components, where thermal conductivity is important, highlighting the need for a detailed understanding of its phonons. However, the phonon dispersion derived from previous inelastic neutron scattering data is in disagreement with modern theoretical calculations, pointing to the necessity for further experimental validation. Here, in this study, we use inelastic neutron scattering experiments on polycrystalline ZrN at 6, 100, 300, 400, and 600 K, and density functional theory (DFT) calculations to study the momentum-resolved scattering and the phonon density of states. Both experimental and theoretical results reveal a phonon cutoff energy near 65 meV, in stark disagreement with previous single-crystal neutron scattering measurements indicating a cutoff near 75 meV. The 65 meV cutoff energy is, however, in reasonable agreement with previous measurements on ZrN 0.9 . Consistency between our theoretical and experimental results confirms that ZrN is a wide phonon band gap system and that DFT provides a reliable description of its phonons.

DFT↗
DOE OSTI2021

Inelastic Neutron Spectra of Uranium Tetrafluoride Hydrate, UF 4 (H 2 O) 2.5

Uranium tetrafluoride hydrate (UFH) is formed by immersing anhydrous UF4 under water for 12 h. UFH is therefore clearly a chemical species of environmental concern, as anhydrous UF4 is an intermediate uranium form in the nuclear fuel cycle. We use inelastic neutron scattering (INS) to probe the full vibrational spectra of UFH and its deuterated analogue in an effort to improve the fundamental understanding of its vibrational spectra. Coupled with density functional theory (DFT) calculations, the first for this compound, and full spectral modeling, we generate the complete vibrational spectra of UFH and compare them to prior optical spectroscopic results. In particular, the combination of DFT with INS allows us to identify multiple distinct chemical environments in the water bending and OH stretching regions. Whereas the water molecules directly bound to the U atoms execute OH stretching around 3600 cm –1 , a second class of H-bonded waters vibrate below 3000 cm –1 , an indicator of strong H bonding. In addition, a class of librational water modes are observed between 400 and 900 cm–1, which themselves can be separated in energy according to their chemical environments. Furthermore, measurements presented herein directly assist in the assignment of certain spectral features in the infrared spectrum and will inform future investigations of UFH for environmental or forensics purposes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Machine learning-based interatomic potential development and phase transition analysis of ferroelectric hafnium dioxide

The ferroelectric phase (𝑃⁢𝑐⁢𝑎⁢2 1 , which is in orthorhombic symmetry) of hafnium dioxide (HfO 2 ) has gained much attention due to its potential applications in nanoelectronics and advanced memory devices. However, its complex phase behavior under external stimuli, such as pressure and temperature, remains a subject of intense investigation. This study focuses on developing a machine learning-based interatomic potential (MLIP) that is trained with data from density-functional theory (DFT) calculations to simulate phase transitions and mechanical properties of HfO 2 . The developed MLIP predicts lattice parameters, equations of state, bulk and shear moduli, and elastic constants that closely align with DFT predictions for several phases and at various pressures. Once validated, the MLIP is used to investigate the phase transitions of ferroelectric HfO 2 (𝑃⁢𝑐⁢𝑎⁢2 1 ) under both isobaric and constant stress conditions at elevated temperatures ranging from 200 to 2500 K. We used several complementary methods, including local symmetry identification, radial distribution function, and x-ray diffraction characterization, to identify interesting phase transitions among several competitive hafnia phases predicted from our simulations. The suggested methods uniformly reveal that under pure deviatoric condition, the system favors a transition from the orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase to a tetragonal (𝑃⁢4 2 /𝑛⁢𝑚⁢𝑐) phase, whereas a zero stress condition drives the system from the 𝑃⁢𝑐⁢𝑎⁢2 1 phase to another orthorhombic (𝑃⁢𝑏⁢𝑐⁢𝑛) phase. These findings provide crucial insights into stress and temperature-induced phase behavior of hafnia, guiding future experimental and theoretical studies for optimizing hafnia-based ferroelectric devices.

Ferroelectric HfO2↗
DOE OSTI2022

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M –1 cm –1 ) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Density Functional Theory Guided Investigation of Ligand‐Induced Neptunyl‐Neptunyl Interactions

Abstract Actinyl‐actinyl interactions are particularly prevalent for the pentavalent neptunyl cation (Np(V)O 2 ) + where these interactions appear either as a T ‐ or D ‐shape (diamond‐shape). T ‐shaped interactions have been previously identified in high concentration Np(V) solutions containing simple anions (NO 3 − , ClO 4 − , Cl − ) whereas D ‐shaped have only been isolated in the solid‐state in the presence of carboxylate ligands. In this study, Density Functional Theory (DFT) calculations were paired with Raman spectroscopy to evaluate the formation of D ‐shaped interactions in the presence of aliphatic (R=H (formate), CH 3 (acetate), CH 2 CH 3 (propionate)) and aromatic (R=C 6 H 5 (benzoate), C 6 H 4 OH (4‐hydroxybenzoate), C 5 H 4 N (isonicotinate)) carboxylate ligands. DFT studies indicate that the Δ G to form hydrated T ‐ and D ‐shaped forms are not spontaneous but become so with the addition of the carboxylate ligands. Raman spectra of the Np(V) carboxylate solutions contained vibrational modes associated with the D ‐shaped interactions, but spectral changes observed over time indicate a dynamic system. Crystallization experiments from the Np(V) carboxylate systems confirmed the presence of D‐ shaped dimers for the aromatic carboxylates, suggesting that the choice of the anion in solution favors actinyl‐actinyl interactions even at low concentrations (≤20 mM) of Np(V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Neural network emulation of spontaneous fission

Large-scale computations of fission properties are an important ingredient for nuclear reaction network calculations simulating rapid neutron-capture process (the 𝑟 process) nucleosynthesis. Due to the large number of fissioning nuclei potentially contributing to the 𝑟 process, a microscopic description of fission based on nuclear density functional theory (DFT) is computationally challenging. Here, we explore the use of neural networks (NNs) to construct DFT emulators capable of predicting potential energy surfaces and collective inertia tensors across the whole nuclear chart, starting from a minimal set of DFT calculations. We use constrained Hartree-Fock-Bogoliubov (HFB) calculations to predict the potential energy and collective inertia tensor in the axial quadrupole and octupole collective coordinates, for a set of nuclei in the 𝑟-process region. We then employ NNs to emulate the HFB energy and collective inertia tensor across the considered region of the nuclear chart. Least-action pathways characterizing spontaneous fission half-lives and fragment yields are then obtained by means of the nudged elastic band method. The potential energy predicted by NNs agrees with the DFT value to within a root-mean-square error of 500 keV, and the collective inertia components agree to within an order of magnitude. These results are largely independent of the NN architecture. The exit points on the outer turning line are found to be well emulated. For the spontaneous fission half-lives the NN emulation provides values that are found to agree with the DFT predictions within a factor of 10 3 across more than 70 orders of magnitude. Neural networks are able to emulate the potential energy and collective inertia well enough to reasonably predict physical observables. Future directions of study, such as the inclusion of additional collective degrees of freedom and active learning, will improve the predictive power of microscopic theory and further enable large-scale fission studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Intricate short-range order in GeSn alloys revealed by atomistic simulations with highly accurate and efficient machine-learning potentials

GeSn alloys hold promise for silicon-compatible integrated applications in electronics, photonics, and topological quantum devices. However, understanding their intricate structures using density functional theory (DFT) calculations is hindered by spatiotemporal constraints. To overcome this limitation, we develop highly accurate and efficient machine-learning interatomic potentials based on a neuroevolution potential approach with farthest point sampling on a comprehensive DFT data set. The application of the developed machine-learning potential in large-scale atomistic simulations bridges the spatiotemporal gap between modeling and advanced characterizations, and facilitates the discovery of structural intricacies in GeSn alloys. Through extensive statistical sampling, we identify a type of short-range order (SRO) that is distinguished by both its structural signature and electronic band gap from the SRO structure previously predicted. Modeling based on a large simulation cell reveals the coexistence of nano SRO domains with various degrees of ordering, demonstrating a complex spatial heterogeneity of SRO structure. Finally, our study not only reinforces the significance of fine-level structural information in alloys, but it also constitutes an effective framework for exploring SRO in a broad range of complex alloys based on highly accurate and effective machine-learning potentials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Characterizing Defects Inside Hexagonal Boron Nitride Using Random Telegraph Signals in van der Waals 2D Transistors

Single-crystal hexagonal boron nitride (hBN) is used extensively in many two-dimensional electronic and quantum devices, where defects significantly impact performance. Therefore, characterizing and engineering hBN defects are crucial for advancing these technologies. Here, we examine the capture and emission dynamics of defects in hBN by utilizing low-frequency noise (LFN) spectroscopy in hBN-encapsulated and graphene-contacted MoS 2 field-effect transistors (FETs). The low disorder of this heterostructure allows the detection of random telegraph signals (RTS) in large device dimensions of 100 μm 2 at cryogenic temperatures. Analysis of gate bias- and temperature-dependent LFN data indicate that RTS originates from a single trap species within hBN. By performing multi-space density functional theory (MS-DFT) calculations on a gated defective hBN/MoS 2 heterostructure model, we assign substitutional carbon atoms in boron sites as the atomistic origin of RTS. This study demonstrates the utility of LFN spectroscopy combined with MS-DFT analysis on a low-disorder all-vdW FET as a powerful means for characterizing the atomistic defects in single-crystal hBN.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

A critical examination of compound stability predictions from machine-learned formation energies

Machine learning has emerged as a novel tool for the efficient prediction of material properties, and claims have been made that machine-learned models for the formation energy of compounds can approach the accuracy of Density Functional Theory (DFT). The models tested in this work include five recently published compositional models, a baseline model using stoichiometry alone, and a structural model. By testing seven machine learning models for formation energy on stability predictions using the Materials Project database of DFT calculations for 85,014 unique chemical compositions, we show that while formation energies can indeed be predicted well, all compositional models perform poorly on predicting the stability of compounds, making them considerably less useful than DFT for the discovery and design of new solids. Most critically, in sparse chemical spaces where few stoichiometries have stable compounds, only the structural model is capable of efficiently detecting which materials are stable. The nonincremental improvement of structural models compared with compositional models is noteworthy and encourages the use of structural models for materials discovery, with the constraint that for any new composition, the ground-state structure is not known a priori. This work demonstrates that accurate predictions of formation energy do not imply accurate predictions of stability, emphasizing the importance of assessing model performance on stability predictions, for which we provide a set of publicly available tests.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-pressure characterization of Ag 3 AuTe 2 : Implications for strain-induced band tuning

Recent band structure calculations have suggested the potential for band tuning in the chiral semiconductor Ag 3 AuTe 2 to zero upon application of negative strain. In this study, we report on the synthesis of polycrystalline Ag 3 AuTe 2 and investigate its transport and optical properties and mechanical compressibility. Transport measurements reveal the semiconducting behavior of Ag 3 AuTe 2 with high resistivity and an activation energy E a of 0.2 eV. The optical bandgap determined by diffuse reflectance measurements is about three times wider than the experimental E a ⁠. Despite the difference, both experimental gaps fall within the range of predicted bandgaps by our first-principles density functional theory (DFT) calculations employing the Perdew–Burke–Ernzerhof and modified Becke–Johnson methods. Furthermore, our DFT simulations predict a progressive narrowing of the bandgap under compressive strain, with a full closure expected at a strain of –4% relative to the lattice parameter. To evaluate the feasibility of gap tunability at such substantial strain, the high-pressure behavior of Ag 3 AuTe 2 was investigated by in situ high-pressure x-ray diffraction up to 47 GPa. Mechanical compression beyond 4% resulted in a pressure-induced structural transformation, indicating the possibility of substantial gap modulation under extreme compression conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning for impurity charge-state transition levels in semiconductors from elemental properties using multi-fidelity datasets

Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). Further, the model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II–VI and III–V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III–V and II–VI systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Rhodium Single–Atom Catalysts on Titania for Reverse Water Gas Shift Reaction Explored by First Principles Mechanistic Analysis and Compared to Nanoclusters

The thermocatalytic reduction of CO2 by H2 often proceeds via two competing reaction mechanisms—the reverse water gas shift reaction (rWGSR, CO2 + H2 ⇌ CO + H2O) and methanation (CO2 + 4H2 ⇌ CH4 + 2H2O). Atomically dispersed Rh1 catalysts on TiO2 show high selectivity toward the rWGSR compared with larger Rh nanoclusters, but the origin of this size-dependent selectivity has not been fully explained. Here we report density functional theory (DFT) calculations and microkinetic simulations that clarify the Rh1 active sites and rWGSR pathway on anatase TiO2(101), as well as the high rWGSR selectivity of Rh1 compared with supported Rhx (x = 2–8 atoms) nanoclusters. DFT-computed formation energies, vibrational frequency analysis, and microkinetic modeling suggest three plausible active sites: Rh1 on titania (Rh1/TiO2(101)), Rh1 with a nearby hydroxyl group (Rh1OH/TiO2(101)), and Rh1 near an oxygen vacancy at a three-fold coordinated site (Rh1 near O3cvac). Predicted turnover frequencies and apparent activation barriers for Rh1 indicate a faster reaction involving CO2 dissociation assisted by a support oxygen vacancy via Rh1 near O3cvac, as well as slower reactions involving Rh1OH/TiO2(101) or Rh1/TiO2(101) through a COOH intermediate. Furthermore, these Rh1 sites are selective toward CO rather than CH4 because of the weak adsorption of CO, large barrier for C-O bond dissociation, and the lack of nearby metal sites for H2 dissociation, in contrast to Rhx nanoclusters, including Rh2 dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗