On the Capabilities of Transition Metal Carbides for Carbon Capture and Utilization Technologies.

The search for cheap and active materials for the capture and activation of CO2 has led to many efforts aimed at developing new catalysts. In this context, earth-abundant transition metal carbides (TMCs) have emerged as promising candidates, garnering increased attention in recent decades due to their exceptional refractory properties and resistance to sintering, coking, and sulfur poisoning. In this work, we assess the use of Group 5 TMCs (VC, NbC, and TaC) as potential materials for carbon capture and sequestration/utilization technologies by combining experimental characterization techniques, first-principles-based multiscale modeling, vibrational analysis, and catalytic experiments. Our findings reveal that the stoichiometric phase of VC exhibits weak interactions with CO2, displaying an inability to adsorb or dissociate it. However, VC often exhibits the presence of surface carbon vacancies, leading to significant activation of CO2 at room temperature and facilitating its catalytic hydrogenation. In contrast, stoichiometric NbC and TaC phases exhibit stronger interactions with CO2, capable of adsorbing and even breaking of CO2 at low temperatures, particularly notable in the case of TaC. Nevertheless, NbC and TaC demonstrate poor catalytic performance for CO2 hydrogenation. This work suggests Group 5 TMCs as potential materials for CO2 abatement, emphasizes the importance of surface vacancies in enhancing catalytic activity and adsorption capability, and provides a reference for using the infrared spectra as a unique identifier to detect oxy-carbide phases or surface C vacancies within Group 5 TMCs.

Chemical bonding and electronic properties along Group 13 metal oxides.

CONTEXT: The present work provides a systematic theoretical analysis of the nature of the chemical bond in Al2O3, Ga2O3, and In2O3 group 13 cubic crystal structure metal oxides. The influence of the functional in the resulting band gap is assessed. The topological analysis of the electron density provides unambiguous information about the degree of ionicity along the group which is linearly correlated with the band gap values and with the cost of forming a single oxygen vacancy. Overall, this study offers a comprehensive insight into the electronic structure of metal oxides and their interrelations. This will help researchers to harness information effectively, boosting the development of novel metal oxide catalysts or innovative methodologies for their preparation. METHODS: Periodic density functional theory was used to predict the atomic structure of the materials of interest. Structure optimization was carried out using the PBE functional, using a plane wave basis set and the PAW representation of the atomic cores, using the VASP code. Next, the electronic properties were computed by carrying out single point calculations employing PBE, PBE + U functionals using VASP and also with PBE and the hybrid HSE06 functionals using the FHI-AIMS software. For the hybrid HSE06, the impact of the screening parameter, ω, and mixing parameter, α, on the calculated band gap has also been assessed.

Atomic and Electronic Structures of Co-Doped In2O3 from Experiment and Theory.

The synthesis and properties of stoichiometric, reduced, and Co-doped In2O3 are described in the light of several experimental techniques, including X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), ultraviolet (UV)-visible spectroscopy, porosimetry, and density functional theory (DFT) methods on appropriate models. DFT-based calculations provide an accurate prediction of the atomic and electronic structure of these systems. The computed lattice parameter is linearly correlated with the experimental result in the Co concentration ranging from 1.0 to 5.0%. For higher Co concentrations, the theoretical-experimental analysis of the results indicates that the dopant is likely to be preferentially present at surface sites. The analysis of the electronic structure supports the experimental assignment of Co2+ for the doped material. Experiments and theory find that the presence of Co has a limited effect on the material band gap.

On the CO2 Harvesting from N2 Using Grazyne Membranes.

The separation of carbon dioxide (CO2) from nitrogen (N2) is at the core of any global warming remediation technology aimed at reducing the CO2 content in the atmosphere. Chemical membranes designed to differentially permeate both molecules have become quite appealing due to their simple use, although many membrane-based separations stand out as a promising solution for CO2 separation. These are environmentally friendly, with high active surface areas, compact design, easy to maintain and cost-effective, although the field is still growing due to the difficulties in the CO2/N2 separation. The present study poses grazynes, two-dimensional C-based materials with sp and sp2 C atoms, aligned along stripes, as suited membranes for the CO2/N2 separation. The combination of density functional theory (DFT) and molecular dynamics (MD) simulations allow tackling the energetics, kinetics, and dynamics of the membrane effectiveness of grazynes with engineered pores for such a separation in a holistic fashion. The explored grazynes are capable of physisorbing CO2 and N2, thus avoiding material poisoning by molecular decoration, while the diffusion of CO2 through the pores is found to be rapid, yet easier than that of N2, in the rate order of the s-1 in the 100-500 K temperature range. In particular, low-temperature CO2 separation even for CO2 contents below 0.5 % are found for [1],[2]{2}-grazyne when controlling the membrane exposure contact to the gas mixture, paving the way for exploring and using grazynes for air CO2 remediation.

Understanding the Chemical Bond in Semiconductor/MXene Composites: TiO2 Clusters Anchored on the Ti2C MXene Surface.

First-principles calculations on titania clusters (TiO2)n (n=5 and 10) supported on the pristine Ti2C (0001) surface were carried out to understand the properties of semiconductor/MXene composites with implications in (photo)-catalysis. The reported results reveal a high exothermic interaction accompanied by a substantial charge transfer with a concomitant, notorious, deformation of the titania nanoclusters. The analysis of the density of states analysis of the composite systems evidences a metallic character with titania related states crossing the Fermi level. The picture of the chemical bonds is completed by the analysis of X-Ray Photoelectron Spectra (XPS) features, evidencing clear shifts of the C(1s) and O(1s) related peaks relative to the isolated systems that have a quite complex origin. This detailed analysis provides insights to experimentalists interested in the design and synthesis of these systems with possible applications in catalysis.

The nature of the electronic ground state of M2C (M = Ti, V, Cr, Zr, Nb, Mo, Hf, Ta, and W) MXenes.

A systematic computational study is presented aimed at accurately describing the electronic ground state nature and properties of M2C (M = Ti, V, Cr, Zr, Nb, Mo, Hf, Ta, and W) MXenes. Electronic band structure calculations in the framework of density functional theory (DFT), carried out with different types of basis sets and employing the generalized gradient approach (GGA) and hybrid functionals, provide strong evidence that Ti2C, Zr2C, Hf2C, and Cr2C MXenes exhibit an open-shell conducting ground state with localized spins on the metal atoms, while V2C, Nb2C, Mo2C, Ta2C, and W2C MXenes exhibit a diamagnetic conducting ground state. For Ti2C, Zr2C, Hf2C, and Cr2C, the analysis of the low-lying spin polarized solutions with different spin orderings indicates that their ground states are antiferromagnetic (AFM), consisting of two ferromagnetic (FM) metal layers coupled antiferromagnetically. For the diamagnetic MXenes, the converged spin polarized solutions are significantly less stable than the closed shell solution except for the case of V2C and Mo2C where those excited open shell solutions can be thermally accessible (less than 300 meV per formula unit). The analysis of charge and spin density distributions of the ground state of the MXenes reveals that, in all cases, the metal atoms have a net charge close to +1 e and C atoms close to -2 e. In the case of diamagnetic MXenes, the electronic structure of V2C, Nb2C, and Ta2C is consistent with metal atoms exhibiting a closed-shell s2d2 configuration whereas for Mo2C, and W2C is consistent with a low-spin s1d4 configuration although the FM solution is close in energy for V2C and Mo2C suggesting that they may play a role in their chemistry at high temperature. For the open shell MXenes, the spin density primarily located at the metal atoms showing one unpaired electron per Ti+, Zr+, and Hf+ magnetic center, consistent with s2d1 configuration of the metal atom, and of ∼3.5 unpaired electrons per Cr+ magnetic center interpreted as a mixture of s2d3 and high-spin s1d4 configuration. Finally, the analysis of the density of states reveals the metallic character of all these bare MXenes, irrespective of the nature of the ground state, with significant covalent contributions for Mo2C and W2C.

Quantifying the Accuracy of Density Functionals on Transition Metal Bulk and Surface Properties.

Density functional theory would be exact when the exact exchange-correlation (xc) functional would be known, but since it is regretfully not known, dozens of xc functionals have been developed in the past decades, with some of them better suited for describing certain systems and/or properties. For transition metals (TMs), recent systematic studies assessing bulk properties─shortest interatomic bond distance, δ, cohesive energy, Ecoh, and bulk modulus, B0─and surface features─surface energy, γ, work function, ϕ, and interlayer distances, δij─of 27 TM bulks and 81 TM surfaces, highlighted that generalized gradient approximation (GGA) based xc functionals are, overall, better suited than other types of xc functionals for the TMs bulk and surfaces description, such as Perdew-Burke-Ernzerhof (PBE) or Vega-Viñes (VV). Still, some basic local density approximation xc functionals were not assessed, such as the Hedin-Lundqvist (HL) and Perdew-Zunger (PZ), or GGAs such as the revised Perdew-Burke-Ernzerhof (revPBE) or the Armiento-Mattsson (AM05). Here, we expand the analysis by not only including them but also the recent meta-GGA strongly constrained appropriately normed (SCAN) xc functional, characterized by fulfilling all 17 mathematical conditions an xc must comply, plus the Bayesian error estimation functional (BEEF) xc, a functional parametrized over a large and diverse set of experimental results using machine learning. The present results reveal that none of the xc studied excel neither PBE nor VV, yet AM05 and SCAN performance is quite acceptable, while BEEF xc probably needs more shells of parametrization to reach competitive accuracy levels.

Machine Learning-Driven Discovery of Key Descriptors for CO2 Activation over Two-Dimensional Transition Metal Carbides and Nitrides.

Fusing high-throughput quantum mechanical screening techniques with modern artificial intelligence strategies is among the most fundamental ─yet revolutionary─ science activities, capable of opening new horizons in catalyst discovery. Here, we apply this strategy to the process of finding appropriate key descriptors for CO2 activation over two-dimensional transition metal (TM) carbides/nitrides (MXenes). Various machine learning (ML) models are developed to screen over 114 pure and defective MXenes, where the random forest regressor (RFR) ML scheme exhibits the best predictive performance for the CO2 adsorption energy, with a mean absolute error ± standard deviation of 0.16 ± 0.01 and 0.42 ± 0.06 eV for training and test data sets, respectively. Feature importance analysis revealed d-band center (εd), surface metal electronegativity (χM), and valence electron number of metal atoms (MV) as key descriptors for CO2 activation. These findings furnish a fundamental basis for designing novel MXene-based catalysts through the prediction of potential indicators for CO2 activation and their posterior usage.

Importance of broken geometric symmetry of single-atom Pt sites for efficient electrocatalysis.

Platinum single-atom catalysts hold promise as a new frontier in heterogeneous electrocatalysis. However, the exact chemical nature of active Pt sites is highly elusive, arousing many hypotheses to compensate for the significant discrepancies between experiments and theories. Here, we identify the stabilization of low-coordinated PtII species on carbon-based Pt single-atom catalysts, which have rarely been found as reaction intermediates of homogeneous PtII catalysts but have often been proposed as catalytic sites for Pt single-atom catalysts from theory. Advanced online spectroscopic studies reveal multiple identities of PtII moieties on the single-atom catalysts beyond ideally four-coordinated PtII-N4. Notably, decreasing Pt content to 0.15 wt.% enables the differentiation of low-coordinated PtII species from the four-coordinated ones, demonstrating their critical role in the chlorine evolution reaction. This study may afford general guidelines for achieving a high electrocatalytic performance of carbon-based single-atom catalysts based on other d8 metal ions.

How does thickness affect magnetic coupling in Ti-based MXenes.

The magnetic nature of Ti2C, Ti3C2, and Ti4C3 MXenes is determined from periodic calculations within density functional theory and using the generalized gradient approximation based PBE functional, the PBE0 and HSE06 hybrids, and the on-site Hubbard corrected PBE+U one, in all cases using a very tight numerical setup. The results show that all functionals consistently predict a magnetic ground state for all MXenes, with spin densities mainly located at the Ti surface atoms. The analysis of solutions corresponding to different spin orderings consistently show that all functionals predict an antiferromagnetic conducting ground state with the two ferromagnetic outer (surface) Ti layers being antiferromagnetically coupled. A physically meaningful spin model is proposed, consistent with the analysis of the chemical bond, with closed shell, diamagnetic, Ti2+ like ions in inner layers and surface paramagnetic Ti+ like centers with one unpaired electron per magnetic center. From a Heisenberg spin model, the relevant isotropic magnetic coupling constants are extracted from an appropriate mapping of total energy differences per formula unit to the expected energy values of the spin Hamiltonian. While the numerical values of the magnetic coupling constants largely depend on the used functional, the nearest neighbor intralayer coupling is found to be always ferromagnetic, and constitutes the dominant interaction, although two other non-negligible interlayer antiferromagnetic terms are involved, implying that the spin description cannot be reduced to NN interaction only. The influence of the MXene thickness is noticeable for the dominant ferromagnetic interaction, increasing its value with the MXene width. However, the interlayer interactions are essentially due to the covalency effects observed in all metallic solutions which, as expected, decay with distance. Within the PBE+U approach, a U value of 5 eV is found to closely simulate the results from hybrid functionals for Ti2C and less accurately for Ti3C2 and Ti4C3.

Toward a Rigorous Theoretical Description of Photocatalysis Using Realistic Models.

This Perspective aims at providing a road map to computational heterogeneous photocatalysis highlighting the knowledge needed to boost the design of efficient photocatalysts. A plausible computational framework is suggested focusing on static and dynamic properties of the relevant excited states as well of the involved chemistry for the reactions of interest. This road map calls for explicitly exploring the nature of the charge carriers, the excited-state potential energy surface, and its time evolution. Excited-state descriptors are introduced to locate and characterize the electrons and holes generated upon excitation. Nonadiabatic molecular dynamics simulations are proposed as a convenient tool to describe the time evolution of the photogenerated species and their propagation through the crystalline structure of photoactive material, ultimately providing information about the charge carrier lifetime. Finally, it is claimed that a detailed understanding of the mechanisms of heterogeneously photocatalyzed reactions demands the analysis of the excited-state potential energy surface.

Theoretical Analysis of Magnetic Coupling in the Ti2C Bare MXene.

The nature of the electronic ground state of the Ti2C MXene is unambiguously determined by making use of density functional theory-based calculations including hybrid functionals together with a stringent computational setup providing numerically converged results up to 1 meV. All the explored density functionals (i.e., PBE, PBE0, and HSE06) consistently predict that the Ti2C MXene has a magnetic ground state corresponding to antiferromagnetic (AFM)-coupled ferromagnetic (FM) layers. A spin model, with one unpaired electron per Ti center, consistent with the nature of the chemical bond emerging from the calculations, is presented in which the relevant magnetic coupling constants are extracted from total energy differences of the involved magnetic solutions using an appropriate mapping approach. The use of different density functionals enables us to define a realistic range for the magnitude of each of the magnetic coupling constants. The intralayer FM interaction is the dominant term, but the other two AFM interlayer couplings are noticeable and cannot be neglected. Thus, the spin model cannot be reduced to include nearest-neighbor interactions only. The Néel temperature is roughly estimated to be in the 220 ± 30 K, suggesting that this material can be used in practical applications in spintronics and related fields.

MXenes à la Carte: Tailoring the Epitaxial Growth Alternating Nitrogen and Transition Metal Layers.

A high-throughput analysis based on density functional simulations underscores the viable epitaxial growth of MXenes by alternating nitrogen and metal adlayers. This is supported by an exhaustive analysis of a number of thermodynamic and kinetic thresholds belonging to different critical key steps in the course of the epitaxial growth. The results on 18 pristine N- and C-based MXenes with M2X stoichiometry reveal an easy initial N2 fixation and dissociation, where N2 adsorption is controlled by the MXene surface charge and metal d-band center and its dissociation controlled by the reaction energy change. Furthermore, formation energies indicate the plausible formation of N-terminated M2XN2 MXenes. Moreover, the further covering with metal adlayers is found to be thermodynamically driven and stable, especially when using early transition metal atoms. The most restrictive analyzed criterion is the N2 adsorption and dissociation at nearly full N-covered adlayers, which is yet achievable for almost half of the explored M2X seeds. The present results unfold the possibility of expanding, controlling, and tuning the composition, width, and structure of the MXene family.

Effect of nanostructuring on the interaction of CO2 with molybdenum carbide nanoparticles.

Transition metal carbides are increasingly used as catalysts for the transformation of CO2 into useful chemicals. Recently, the effect of nanostructuring of such carbides has started to gain relevance in tailoring their catalytic capabilities. Catalytic materials based on molybdenum carbide nanoparticles (MoCy) have shown a remarkable ability to bind CO2 at room temperature and to hydrogenate it into oxygenates or light alkanes. However, the involved chemistry is largely unknown. In the present work, a systematic computational study is presented aiming to elucidate the chemistry behind the bonding of CO2 with a representative set of MoCy nanoparticles of increasing size, including stoichiometric and non-stoichiometric cases. The obtained results provide clear trends to tune the catalytic activity of these systems and to move towards more efficient CO2 transformation processes.

Catalytic Reduction of Carbon Dioxide on the (001), (011), and (111) Surfaces of TiC and ZrC: A Computational Study.

We present a computational study of the activity and selectivity of early transition-metal carbides as carbon dioxide reduction catalysts. We analyze the effects of the adsorption of CO2 and H2 on the (001), (011), and metal-terminated (111) surfaces of TiC and ZrC, as carbon dioxide undergoes either dissociation to CO or hydrogenation to COOH or HCOO. The relative stabilities of the three reduction intermediates and the activation energies for their formation allow the identification of favored pathways on each surface, which are examined as they lead to the release of CO, HCOOH, CH3OH, and CH4, thereby also characterizing the activity and selectivity of the two materials. Reaction energetics implicate HCO as the key common intermediate on all surfaces studied and rule out the release of formaldehyde. Surface hydroxylation is shown to be highly selective toward methane production as the formation of methanol is hindered on all surfaces by its barrierless conversion to CO.

Artificial-intelligence-driven discovery of catalyst genes with application to CO2 activation on semiconductor oxides.

Catalytic-materials design requires predictive modeling of the interaction between catalyst and reactants. This is challenging due to the complexity and diversity of structure-property relationships across the chemical space. Here, we report a strategy for a rational design of catalytic materials using the artificial intelligence approach (AI) subgroup discovery. We identify catalyst genes (features) that correlate with mechanisms that trigger, facilitate, or hinder the activation of carbon dioxide (CO2) towards a chemical conversion. The AI model is trained on first-principles data for a broad family of oxides. We demonstrate that surfaces of experimentally identified good catalysts consistently exhibit combinations of genes resulting in a strong elongation of a C-O bond. The same combinations of genes also minimize the OCO-angle, the previously proposed indicator of activation, albeit under the constraint that the Sabatier principle is satisfied. Based on these findings, we propose a set of new promising catalyst materials for CO2 conversion.

Can calculated harmonic vibrational spectra rationalize the structure of TiC-based nanoparticles?

Nanoscale titanium carbide (TiC) is widely used in composites and energy applications. In order to design and optimize these systems and to gain a fundamental understanding of these nanomaterials, it is important to understand the atomistic structure of nano-TiC. Cluster beam experiments have provided detailed infrared vibrational spectra of numerous TixCy nanoparticles with well defined masses. However, these spectra have yet to be convincingly assigned to TixCy nanoparticle structures. Herein, using accurate density functional theory based calculations, we perform a systematic survey of likely candidate nanoparticle structures with masses corresponding to those in experiment. We calculate harmonic infrared vibrational spectra for a range of nanoparticles up to 100 atoms in size, with a focus on systems based on removing either four carbon atoms or a single titanium atom from rocksalt-structured stoichiometric TiC nanoparticles. Our calculations clearly show that Ti-deficient nanoparticles are unlikely candidates to explain the experimental spectra as such structures are highly susceptible to C-C bonding, whose characteristic frequencies are not observed in experiment. However, our calculated infrared spectra for C-deficient nanoparticles have some matching features with the experimental spectra but tend to have more complex infrared spectra with more peaks than those obtained from experiment. We suggest that the discrepancy between experiment and theory may be largely due to thermally induced anharmonicities and broadening in the latter nanoparticles, which are not be accounted for in harmonic vibrational calculations.

Charting the Atomic C Interaction with Transition Metal Surfaces.

Carbon interaction with transition metal (TM) surfaces is a relevant topic in heterogeneous catalysis, either for its poisoning capability, for the recently attributed promoter role when incorporated in the subsurface, or for the formation of early TM carbides, which are increasingly used in catalysis. Herein, we present a high-throughput systematic study, adjoining thermodynamic plus kinetic evidence obtained by extensive density functional calculations on surface models (324 diffusion barriers located on 81 TM surfaces in total), which provides a navigation map of these interactions in a holistic fashion. Correlation between previously proposed electronic descriptors and ad/absorption energies has been tested, with the d-band center being found the most suitable one, although machine learning protocols also underscore the importance of the surface energy and the site coordination number. Descriptors have also been tested for diffusion barriers, with ad/absorption energies and the difference in energy between minima being the most appropriate ones. Furthermore, multivariable, polynomial, and random forest regressions show that both thermodynamic and kinetic data are better described when using a combination of different descriptors. Therefore, looking for a single perfect descriptor may not be the best quest, while combining different ones may be a better path to follow.

Identifying the Atomic Layer Stacking of Mo2C MXene by Probe Molecule Adsorption.

A density functional theory study is presented here aimed at investigating whether the atomic stacking on the new family of two-dimensional MXene materials has an influence on their adsorption properties and whether these properties can provide information about this structural feature. To this end, the Mo2C MXene, exhibiting two nearly degenerate crystal structures with either ABC or ABA atomic stacking, is chosen as a case study. The study of the adsorption of CO, CO2, and H2O on both polymorphs of Mo2C reveals substantial differences that could be used in experiments to provide information about the atomic stacking of a given sample. Particularly, we show that the asymmetric and symmetric stretching modes of the adsorbed CO2 and the CO stretching mode are clear features that allow one to identify the stacking of atomic layers of the Mo2C MXene. The present finding is likely to apply to other MXenes as well.

Acetylene-Mediated Electron Transport in Nanostructured Graphene and Hexagonal Boron Nitride.

The discovery of graphene has catalyzed the search for other 2D carbon allotropes, such as graphynes, graphdiynes, and 2D π-conjugated polymers, which have been theoretically predicted or experimentally synthesized during the past decade. These materials exhibit a conductive nature bound to their π-conjugated sp2 electronic system. Some cases include sp-hybridized moieties in their nanostructure, such as acetylenes in graphynes; however, these act merely as electronic couplers between the conducting π-orbitals of sp2 centers. Herein, via first-principles calculations and quantum transport simulations, we demonstrate the existence of an acetylene-meditated transport mechanism entirely hosted by sp-hybridized orbitals. For that we propose a series of nanostructured 2D materials featuring linear arrangements of closely packed acetylene units which function as sp-nanowires. Because of the very distinct nature of this unique transport mechanism, it appears to be highly complementary with π-conjugation, thus potentially becoming a key tool for future carbon nanoelectronics.