Rules Describing CO2 Activation on Single-Atom Alloys from DFT-meta-GGA Calculations and Artificial Intelligence

Herzain I. Rivera-Arrieta,Lucas Foppa
DOI: https://doi.org/10.26434/chemrxiv-2024-1dr10
2024-11-19
Abstract:Single-atom alloys (SAAs) arise as a promising concept for the design of improved CO2 hydrogenation catalysts. However, from the immense number of possible SAA compositions and structures, only a few might display the properties required to be useful catalysts. Thus, the direct, high-throughput screening of materials is inefficient. Here, we use artificial intelligence to derive rules describing surface sites of SAAs that provide an effective CO2 activation, a crucial initial step to convert the molecule into valuable products. We start by modeling the CO2 interaction with 780 sites of flat and stepped surfaces of SAAs composed by Cu, Zn, and Pd hosts via high-quality DFT-mBEEF calculations. Then, we apply subgroup discovery to determine constraints on key physical parameters, out of 24 offered candidate descriptive parameters, characterizing subgroups (SGs) of surface sites where chemisorbed CO2 displays large elongations of its C-O bonds. The key identified parameters are free-atom properties of the elements constituting the surface sites, such as their electron affinity, electronegativity, and radii of the $d$-orbitals. Additionally, the generalized coordination number is selected as a key geometrical parameter. The SG rules are used to identify promising alloys among more than 1,500 possible single-atom and dual-atom alloys. Some of the promising alloys predicted by the SG rules were explicitly evaluated by additional DFT-mBEEF calculations and confirmed to provide a significant CO2 activation.
Chemistry
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to address the challenges in the carbon dioxide (CO₂) activation process, especially to design more effective CO₂ hydrogenation catalysts. Specifically, the goals of the paper are: 1. **Improve CO₂ activation efficiency**: CO₂ molecules are relatively inert and difficult to directly convert into valuable chemical products (such as olefins and alcohols). Therefore, efficient catalysts are required to activate CO₂, making it more likely to participate in chemical reactions. 2. **Screen effective catalysts**: Single - Atom Alloys (SAAs), as a promising catalyst material, have unique electronic properties and can adjust the adsorption strength of reaction species on the surface. However, due to the large compositional and structural space of SAAs, direct high - throughput screening is inefficient. Therefore, an efficient method is needed to screen potentially effective SAAs. 3. **Use artificial intelligence to accelerate discovery**: By combining Density Functional Theory (DFT) calculations and artificial intelligence (AI), especially Subgroup Discovery (SGD) techniques, to identify the key features of SAA surface sites that can effectively activate CO₂. This method can not only accelerate the discovery of new materials but also provide chemical insights into CO₂ activation. 4. **Verify prediction results**: Based on the potential alloys predicted by the SGD rules, further verify their CO₂ activation ability through additional DFT - mBEEF calculations to ensure the validity and reliability of the predictions. ### Main methods of the paper - **DFT - mBEEF calculations**: Use meta - Bayesian - error - estimation functional (mBEEF) for DFT simulations to study the interactions of CO₂ on 780 different SAA surface sites. - **Subgroup Discovery (SGD)**: Through SGD technology, determine the key physical parameters from 24 candidate description parameters, which can characterize the significant C−O bond elongation of CO₂ on SAA surface sites. - **Screening and verification**: Screen more than 1,500 possible SAAs and Diatomic Alloys (DAAs) according to SGD rules, and confirm their CO₂ activation ability through additional DFT - mBEEF calculations. ### Key formulas - CO₂ adsorption energy: \[ E_{\text{ads}}^{\text{CO}_2} = E_{\text{slab}\cdot\cdot\cdot\text{CO}_2} - (E_{\text{slab}} + E_{\text{CO}_2}) \] where \( E_{\text{slab}\cdot\cdot\cdot\text{CO}_2} \) is the surface energy after adsorbing CO₂, \( E_{\text{slab}} \) is the energy of the clean surface, and \( E_{\text{CO}_2} \) is the energy of an isolated CO₂ molecule. - C−O bond elongation: \[ \Delta d_{\text{C−O}}^{\text{max}} = d_{\text{C−O}}^{\text{chem}} - d_{\text{C−O}}^{\text{equil}} \] where \( d_{\text{C−O}}^{\text{chem}} \) is the maximum distance between the two C−O bonds in the adsorbed CO₂ molecule, and \( d_{\text{C−O}}^{\text{equil}} \) is the distance of the C−O bond in the optimized gas - phase CO₂ molecule (1.155 Å). Through the above methods, the paper has successfully identified SAAs that can significantly activate CO₂ and provided important guidance for the future design of more efficient CO₂ hydrogenation catalysts.