这是用户在 2025-8-7 3:32 为 https://pubs.acs.org/doi/10.1021/acs.jpcc.3c02483?fig=agr1&ref=pdf 保存的双语快照页面,由 沉浸式翻译 提供双语支持。了解如何保存?

Pair your account to your Institution

Click the pair button to affiliate your institution with your personal account

ACS Publications. Most Trusted. Most Cited. Most Read
Unveiling Static and Dynamic Structures of Pd Clusters Influenced by Al2O3 Surfaces: DFT and AIMD Studies
  • Subscribed
C: Chemical and Catalytic Reactivity at Interfaces

Unveiling Static and Dynamic Structures of Pd Clusters Influenced by Al2O3 Surfaces: DFT and AIMD Studies
揭示受 Al2O3 表面影响的 Pd 簇的静态和动态结构:DFT 和 AIMD 研究
Click to copy article link
点击复制文章链接
Article link copied!

  • Kazumasa Murata
    Kazumasa Murata
    Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
  • Akira Oda
    Akira Oda
    Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
    Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
    More by Akira Oda
  • Atsushi Satsuma
    Atsushi Satsuma
    Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
    Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
  • Kyoichi Sawabe*
    Kyoichi Sawabe
    Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
    Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
    *Email: sawabe@chembio.nagoya-u.ac.jp
Open PDFSupporting Information (19)open URL

The Journal of Physical Chemistry C

Cite this: J. Phys. Chem. C 2023, 127, 41, 20267–20275
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcc.3c02483
Published October 5, 2023
Copyright © 2023 American Chemical Society

Article Views

1510

Altmetric

-

Citations

-
Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

Abstract  抽象

Click to copy section link
点击复制部分链接
Section link copied!

Metal nanoparticles (MNPs) supported on oxides are used in various catalytic reactions, such as petrochemical processes and exhaust gas purification. Metal–support interactions (MSIs) between MNPs and oxides affect the morphology and thermal stability of the MNPs, which are known to be related to the catalytic activity of the MNPs. However, the effect of MSIs is still a matter of debate. We conducted density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) to investigate the geometric structure and migration behavior of Pd clusters on various Al2O3 surfaces. The MSIs of the Pd clusters with different Al2O3 surfaces resulted in different morphologies of the clusters. The shape of Pd clusters on α-Al2O3(0001) were deformed to oblate and on γ-Al2O3(110) to prolate. On the other hand, the shape of Pd clusters on γ-Al2O3(100) did not change from that in a vacuum. The migration rate of Pd clusters was also affected by the MSI, which was higher on γ-Al2O3(110) than on α-Al2O3(0001) and γ-Al2O3(100). To evaluate the interactions between Pd clusters and Al2O3 surfaces, we calculated the interaction energies for Pd clusters with various sizes on Al2O3 surfaces. The interaction energy of Pd/Al2O3 systems depended on the coordination number (CN, number of bonds per Pd atom) of Pd–Al and Pd–O. Despite the smaller atomic density of Al and O of the γ-Al2O3(110) surface than that of α-Al2O3(0001) and γ-Al2O3(100) surfaces, the number of bonds formed between Pd clusters and γ-Al2O3(110) was greater than that of α-Al2O3(0001) and γ-Al2O3(100). According to AIMD simulations and perpendicular displacements of Al and O atoms bound to the Pd cluster on the Al2O3 surface from the bare Al2O3 surface, coordinatively unsaturated Al3+ sites on γ-Al2O3(110) were more flexible than those on α-Al2O3(0001) and γ-Al2O3(100). The flexibility of γ-Al2O3(110) facilitated the formation of bonds with the Pd cluster, resulting in stronger MSIs. When the coordinatively unsaturated Al3+ sites were hydrated, the MSIs were found to be weakened.
负载在氧化物上的金属纳米颗粒 (MNP) 用于各种催化反应,例如石化过程和废气净化。MNPs 和氧化物之间的金属-载体相互作用(MSI)影响 MNPs 的形貌和热稳定性,已知这与 MNPs 的催化活性有关。然而,MSI 的效果仍然是一个有争议的问题。我们通过密度泛函理论(DFT)计算和从头分子动力学(AIMD)研究了 Pd 簇在各种 Al2O3 表面上的几何结构和迁移行为。具有不同 Al2O3 表面的 Pd 簇的 MSIs 导致了簇的不同形貌。α-Al2O3(0001)上的 Pd 簇形状变形为扁圆形,γ-Al2O3(110)上的 Pd 簇形状变形为扁圆形。另一方面,γ-Al2O3(100)上的 Pd 簇形状与真空中的 Pd 簇形状没有变化。Pd 簇的迁移速率也受到 MSI 的影响,MSI 在γ-Al2O3(110)上高于α-Al2O3(0001)和γ-Al2O3(100)。为了评估 Pd 簇与 Al2O3 表面之间的相互作用,我们计算了 Al2O3 表面上各种尺寸的 Pd 簇的相互作用能。Pd/Al2O3 体系的相互作用能取决于 Pd-Al 和 Pd-O 的配位数(CN,每个 Pd 原子的键数)。 尽管γ-Al2O3(110)表面的 Al 和 O 原子密度小于α-Al2O3(0001)和γ-Al2O3(100)表面,但 Pd 簇与γ-Al2O3(110)之间形成的键数大于α-Al2O3(0001)和γ-Al2O3(100)。根据 AIMD 模拟和 Al2O3 表面与 Al2O 3 表面 Pd 簇结合的 Al 和 O 原子从裸露的 Al2O3 表面的垂直位移,γ-Al2O3(110)上的配位不饱和 Al3+ 位点比α-Al2O3(0001)和γ-Al2O3(100)上的位点更灵活。γ-Al2O3(110)的柔韧性促进了与 Pd 簇的键形成,从而产生了更强的 MSI。当配位不饱和的 Al3+ 位点水合时,发现 MSIs 被削弱。

This publication is licensed under the terms of your institutional subscription. Request reuse permissions.
本出版物根据您的条款获得许可 机构订阅。 请求重用权限。

Copyright © 2023 American Chemical Society
版权所有 © 2023 美国化学会

1. Introduction  1. 简介

Click to copy section link
点击复制部分链接
Section link copied!

The morphology of metal nanoparticles (MNPs) is closely related to their catalytic activity; therefore, their control is significant. The relation between the morphology of MNPs and the catalytic activity has been reported for various catalytic reactions. (1) In general, small-sized MNPs have been supported on oxide and carbon supports with a large surface area in order to efficiently utilize the active metal species. The metal–support interactions (MSIs), which are induced at the interface between the MNPs and the support, affect not only the stability of the small MNPs but also the geometric structure of the MNPs. (2) Even for small MNPs or metal clusters, the number of unsaturated coordination sites on the surface varies depending on their morphology. In addition, the electronic structure of the MNPs is significantly affected by their morphology. (3) Thus, understanding MSIs that control the size and morphology of MNPs on the supports leads to catalytic design.
金属纳米颗粒(MNPs)的形貌与其催化活性密切相关;因此,他们的控制很重要。MNPs 的形貌与催化活性之间的关系已被报道用于各种催化反应。(1) 一般来说,为了有效利用活性金属物质,小型 MNPs 被支撑在具有大表面积的氧化物和碳载体上。在 MNPs 和载体之间的界面处诱导的金属-载体相互作用(MSIs)不仅影响小 MNPs 的稳定性,还影响 MNPs 的几何结构。(2) 即使对于小 MNPs 或金属团簇,表面不饱和配位位点的数量也因形貌而异。此外,MNP 的电子结构受其形貌的显着影响。(3) 因此,了解控制载体上 MNP 尺寸和形态的 MSI 可以进行催化设计。
γ-Al2O3 is a nonreducible oxide and commonly used as the support due to its high surface area. The surface of γ-Al2O3 has coordinatively unsaturated Al3+ sites derived from structural defects and shows unique acid–base properties. (4−6) It has been reported that the coordinatively unsaturated Al3+ sites, such as pentacoordinated Al3+ (AlV) sites on the γ-Al2O3(100) surface, are strongly anchoring metal species, and the difference in the strength of MSIs affects the morphology of MNPs. (7−15) At low Pt loadings, isolated Pt atoms were present on the γ-Al2O3(100) surface, and at high Pt loading, raft-like Pt particles were formed. On the other hand, three-dimensional Pt particles were formed on the α-Al2O3 surface without AlV sites. (11) In addition, the structural disorder of a Pt crystal was induced by strong MSIs in Pt/γ-Al2O3, and amorphous Pt particles were generated on γ-Al2O3. (16,17) Pt nanoparticles were unevenly distributed at defect sites such as edges, steps, or kinks in γ-Al2O3. (18) Moreover, Pt particles on γ-Al2O3 were deformed for three-dimensional to a raft-like structure during H2 reduction at high temperature. (19) Hydroxyl groups formed by the hydration of Al2O3 surfaces also contributed to the stabilization of metal species. (20−22)
γ-Al2O3 是一种不可还原的氧化物,由于其高表面积,通常用作载体。γ-Al2O3 的表面具有源自结构缺陷的配位不饱和 Al3+ 位点,并表现出独特的酸碱特性。(4−6) 据报道,配位不饱和的 Al3+ 位点,如γ-Al2O3(100)表面的五配位 Al3+ (AlV)位点,是强锚定金属物质,MSI 的强度差异影响 MNP 的形貌。(7−15) 在低 Pt 负载下,γ-Al2O3 上存在孤立的 Pt 原子(100)表面,在高 Pt 负载下,形成筏状 Pt 颗粒。另一方面,在没有 AlV 位点的α-Al2O3 表面上形成了三维 Pt 颗粒。(11) 此外,Pt/γ-Al2O3 中的强 MSI 诱导了 Pt 晶体的结构紊乱,并在γ-Al2O3 上生成了无定形的 Pt 颗粒。(16,17) 铂纳米颗粒在γ-Al2O3 的边缘、台阶或扭结等缺陷部位分布不均匀。(18) 此外,在高温下 H2 还原过程中,γ-Al2O3 上的 Pt 颗粒被三维变形为筏状结构。(19) 由 Al2O3 表面水合形成的羟基也有助于金属物质的稳定。(20−22)
In the case of Pd supported on Al2O3, the morphology of the Pd particles and their surface structure were changed by the crystalline phase of Al2O3 used as supports. Spherical Pd particles were formed on θ-Al2O3 and α-Al2O3, and they showed a high activity in complete methane oxidation and hydrogenation of unsaturated aldehyde. (23,24) On the other hand, amorphous-like Pd particles with a high fraction of corner sites preferentially were formed on γ-Al2O3, and they were active in the oxidation of carbon monoxide. (25) In addition, the interaction between Pd and Al2O3 stabilized the active Pd phase in catalysts. (26−28) It is inferred that the different Pd–Al2O3 interactions affected the geometric structure and stability of the Pd particles. However, only a few theoretical reports deal with the geometric structures and the stabilization mechanism of Pd particles on different Al2O3 surfaces, and the details are not well understood. (29,30)
在 Al2O3 负载的 Pd 的情况下,Pd 颗粒的形貌及其表面结构因用作载体的 Al2O3 的晶相而改变。在 θ-Al2O3 和α-Al2O3 上形成球形 Pd 颗粒,在甲烷完全氧化和不饱和醛加氢中表现出较高的活性。(23,24) 另一方面,在γ-Al2O3 上优先形成具有高比例角位点的无定形状钯颗粒,它们在一氧化碳的氧化中具有活性。(25) 此外,Pd 与 Al2O3 之间的相互作用稳定了催化剂中的活性 Pd 相。(26−28) 推断不同的 Pd-Al2O3 相互作用影响了 Pd 颗粒的几何结构和稳定性。然而,只有少数理论报告涉及 Pd 颗粒在不同 Al2O3 表面上的几何结构和稳定机理,细节尚不清楚。(29,30)
DFT calculations have been used to understand the interactions between metals and Al2O3. DFT studies show that isolated metal species on the Al2O3 surface were stabilized by coordinatively unsaturated Al3+ sites or terminal hydroxyl groups formed by hydration. (12,20−22) An optimized structure for metal clusters on Al2O3 has been proposed based on screening adsorption sites for metal clusters. (18,29−32) However, as the size of metal clusters increases, metal clusters interact with multiple adsorption sites on the Al2O3. Therefore, it is difficult to determine the thermodynamically stable adsorption sites and orientations for the metal clusters on Al2O3. In the case of Pd/γ-Al2O3, the adsorption structures of Pd clusters consisting of 1–5 atoms and Pd(111) monolayers on nonhydrated γ-Al2O3(100) and hydrated γ-Al2O3(110) were investigated. (30) In addition, the morphology of Pd particles was estimated from the adhesion of the Pd(111) monolayer on the Al2O3 surface using the Young–Dupré model. However, the Young–Dupré model can only reproduce the morphology of MNPs on the oxide on a macroscopic scale. It was reported that the structure of Pd13 clusters optimized in a vacuum was used to perform the geometry optimization of Pd13 clusters on γ-Al2O3. (29) However, these are the optimized structures at 0 K. In real catalysis systems, Pd clusters deform due to the finite temperature. To understand the relation between the morphology of metal clusters and the Al2O3 surfaces, ab initio molecular dynamics (AIMD) is needed to reveal the MSIs at the atomic scale.
DFT 计算已用于了解金属与 Al2O3 之间的相互作用。DFT 研究表明,Al2O3 表面上分离的金属物质被配位不饱和的 Al3+ 位点或水合形成的末端羟基稳定。(12,20−22) 基于金属团簇的筛选吸附位点,提出了一种在 Al2O3 上金属团簇的优化结构。(18,29−32) 然而,随着金属团簇尺寸的增加,金属团簇与 Al2O3 上的多个吸附位点相互作用。因此,很难确定 Al2O3 上金属团簇的热力学稳定吸附位点和取向。在 Pd/γ-Al2O3 中,研究了由 1-5 个原子和 Pd(111)单层组成的 Pd 簇在非水合γ-Al2O3(100)和水合γ-Al2O3(110)上的吸附结构。(30) 此外,使用 Young-Dupré模型从 Pd(111)单层在 Al2O3 表面的粘附度中估计了 Pd 颗粒的形貌。然而,Young-Dupré 模型只能在宏观尺度上再现氧化物上 MNP 的形态。据报道,利用在真空中优化的 Pd13 簇结构在γ-Al2O3 上对 Pd13 簇进行几何优化。(29) 然而,这些是 0 K 下的优化结构。在实际催化系统中,Pd 簇由于有限的温度而变形。 为了了解金属团簇形貌与 Al2O3 表面之间的关系,需要从头开始分子动力学(AIMD)来揭示原子尺度的 MSI。
In this study, we investigate the static and dynamic structures of Pd clusters on various Al2O3 surfaces by a theoretical approach combining DFT calculations and AIMD simulations. During AIMD annealing, Pd clusters migrate on the Al2O3 surfaces with a structural change. By optimizing their structures obtained by AIMD annealing, the morphology of Pd particles on Al2O3 can be proposed with minimal artificial effects. Based on a systematic study of the MSIs for Pd clusters with different sizes on Al2O3, we elucidate the MSIs that depend on the Al2O3 surfaces. The objective of this work is to enhance our understanding of the influence of exposed atomic sites of Al2O3 on the morphology of Pd nanoparticles, which offers significant insights into the strategic development of active metal nanoparticles.
本研究采用 DFT 计算和 AIMD 模拟相结合的理论方法,研究了各种 Al2O3 表面上 Pd 簇的静态和动态结构。在 AIMD 退火过程中,Pd 簇在 Al2O3 表面上迁移并发生结构变化。通过优化 AIMD 退火获得的结构,可以提出 Pd 颗粒在 Al2O3 上的形貌,而人工效应最小。基于对 Al2O3 上不同尺寸的 Pd 簇的 MSI 的系统研究,我们阐明了依赖于 Al2O3 表面的 MSI。这项工作的目的是加强我们对 Al2O3 暴露原子位点对 Pd 纳米颗粒形态影响的理解,这为活性金属纳米颗粒的战略发展提供了重要的见解。

2. Computational Methods  2. 计算方法

Click to copy section link
点击复制部分链接
Section link copied!

Density functional theory (DFT) calculations were performed using the Vienna ab initio simulation package (VASP) code. (33−35) The electron–ion interaction with the frozen-core approximation was described using the projector-augmented wave (PAW) method. (36) Wave functions for valence electrons were represented using a periodic-wave basis set with a cutoff energy of 420 eV. The RPBE functional with the D3(BJ) dispersion correction of Grimme was adopted. (37,38) Structural optimization was performed until the energy and force converged within criteria of 1.0 × 10–5 eV and 0.02 eV/Å, respectively. The vacuum spacing for the isolated Pdn clusters (n = 1–4, 6, 9, and 13) was set to 13 × 14 × 15 Å3. All calculations for Al2O3 surfaces were performed on slab models with eight atomic layers. The dehydrated α-Al2O3(0001) (39) and dehydrated γ-Al2O3(100), dehydrated/hydroxylated γ-Al2O3(110) (40) were modeled with 3 × 3 supercells (14.3 × 12.4 Å2) and 2 × 2 supercells (11.1 × 16.7 Å2), 2 × 2 supercells (16.7 × 16.1 Å2), respectively. In the case of Pd13/γ-Al2O3(100), 3 × 2 supercells of dehydrated γ-Al2O3(100) (16.6 × 16.7 Å) were used to avoid the artificial lateral interaction. Based on the previous thermodynamic analysis, the hydroxyl coverage of hydroxylated γ-Al2O3(110) at 773 K in a vacuum was 5.9 OH nm–2. (40) The vacuum spacing of the Al2O3 slab and Pdn cluster (n = 1–4, 6, 9, and 13) on the Al2O3 slab was set to 14 Å. In the structure optimization of Pdn/Al2O3 (n = 1–4, 6, 9, and 13), an ab initio molecular dynamics (AIMD) annealing procedure described below was employed to reduce artifacts manually placing Pdn/Al2O3 (n = 1–4, 6, 9, and 13). (41) AIMD simulations in the NVT ensemble implemented in VASP were performed by applying the PBE exchange-correlation functional coupled with the D3(BJ) dispersion correction method. The Nosé mass control was set to a value of SMASS = 0. The temperature in the simulation was controlled with a Nosé–Hoover thermostat. First, the temperature was raised to 773 K from 0 K within 2 ps with a time step of 1.0 fs and then held at 773 K for up to 8 ps. Upon simulated annealing, the temperature and potential energy do not exhibit any overshoot (Figure S1). However, this annealing process results in the Pd atoms gaining sufficient energy to significantly alter their structure, and they also undergo bond breaking at the interface with Al and O atoms during the annealing (Movies S1–S3). The Pd clusters migrate on the Al2O3 surface, and their adsorbed positions shift from one symmetrically unique site to another within the surface unit. Regardless of the Al2O3 surface used, the Pd clusters did not remain in their initial arbitrary positions. The changes in the structure and position of the Pd clusters before and after simulated annealing suggest that AIMD annealing eliminates the artificial placement of Pd clusters on the Al2O3 surface. After holding for 3 ps, the system was expected to reach equilibrium based on the variation of the average atomic distances (vide infra). To verify that the temperature distribution is monomodal in equilibrium, the histogram distributions of the temperature by holding from 3 to 8 ps are plotted (Figure S2). All temperature distributions are centered around T = 773 K. In equilibrium, the variance of temperature fluctuations in a canonical ensemble is represented by the equation (42)
使用维也纳从头开始模拟包 (VASP) 代码进行密度泛函理论 (DFT) 计算。(33−35) 使用投影仪增强波 (PAW) 方法描述了电子-离子与冻结核心近似的相互作用。(36) 价电子的波函数使用截止能量为 420 eV 的周期波基集表示。采用 RPBE 函数与 Grimme 的 D3(BJ)色散校正。(37,38) 进行结构优化,直到能量和力分别收敛在 1.0 × 10–5 eV 和 0.02 eV/Å 的标准范围内。将分离的 Pdn 簇(n = 1–4、6、9 和 13)的真空间距设置为 13 × 14 × 15 Å3。Al2O3 表面的所有计算都是在具有 8 个原子层的板坯模型上进行的。分别用 3 个×3 个超级细胞(14.3 × 12.4 Å2)和 2 个×2 超级细胞(11.1 × 16.7 Å 2)、2 × 2 个超级细胞(16.7 × 16.1 Å 2)、2 2 个超级细胞(16.7 16.1 Å2)对脱水γ α-Al2O3(100)、脱水/羟基化γ-Al2O3(100 进行建模。 在 Pd13/γ-Al2O3(100)的情况下,使用了 3 个×2 个脱水γ-Al2O3(100)(16.6 × 16.7 Å)的超级细胞来避免人工横向相互作用。根据之前的热力学分析,羟基化γ-Al2O3(110)在真空中 773 K 的羟基覆盖率为 5.9 OH nm-2(40)Al2O3 板和 Pdn 簇(n = 1-4、6、9 和 13)在 Al2O3 板上的真空间距设置为 14 Å。在 Pdn/Al2O3n = 1–4、6、9 和 13)的结构优化中,采用了下述从头开始分子动力学(AIMD)退火程序来减少手动放置 Pdn/Al2O3n = 1–4、6、9 和 13)的伪影。(41) 采用 PBE 交换相关泛函耦合 D3(BJ)色散校正方法,对 VASP 中实现的 NVT 集合进行 AIMD 模拟。将 Nosé 质量控制设置为 SMASS = 0 的值。模拟中的温度由 Nosé-Hoover 恒温器控制。首先,在 2 ps 内将温度从 0 K 提高到 773 K,时间步长为 1.0 fs,然后保持在 773 K,持续 8 ps。在模拟退火时,温度和势能没有表现出任何过冲( 图 S1)。然而,这种退火过程导致 Pd 原子获得足够的能量来显着改变其结构,并且在退火过程中它们也会在与 Al 和 O 原子的界面处发生键断裂( 电影 S1-S3)。Pd 簇在 Al2O3 表面上迁移,它们的吸附位置从表面单元内的一个对称独特的位点移动到另一个位点。无论使用何种 Al2O3 表面,Pd 簇都不会保持在其初始任意位置。 模拟退火前后 Pd 簇结构和位置的变化表明,AIMD 退火消除了 Pd 簇在 Al2O3 表面的人工放置。在保持 3 ps 后,系统有望根据平均原子距离的变化达到平衡( 见下图 )。为了验证温度分布在平衡状态下是单峰的,绘制了保持 3 至 8 ps 的温度直方图分布( 图 S2)。所有温度分布都以 T = 773 K 为中心。在平衡状态下,规范集合中温度波动的方差由方程 (42) 表示
σ2=2T23N
where T, σ, and N denote temperature, the standard deviation, and the number of free atoms, respectively. This equation yields the following standard deviations at 773 K: 53.5 for Pd13/α-Al2O3(0001), 53.5 for Pd13/γ-Al2O3(100), and 54.7 for Pd13/γ-Al2O3(110). These values are comparable to the standard deviations of the temperatures derived from the 5 ps molecular dynamics (MD) simulation conducted at 773 K, which are 52.6 for Pd13/α-Al2O3(0001), 51.7 for Pd13/γ-Al2O3(100), and 51.3 for Pd13/γ-Al2O3(110). The structure optimization of Pdn/Al2O3 (n = 1–4, 6, 9, and 13) was performed, starting from the initial structure obtained through AIMD annealing. The top four atomic layers of the Al2O3 slab and Pdn cluster (n = 1–4, 6, 9, and 13) were allowed to relax, and the rest of the bottom atoms were fixed. For the sampling of the reciprocal space, Monkhorst–Pack k-point sets of (2 × 2 × 1) and (4 × 3 × 1), (2 × 3 × 1) were sampled for dehydrated α-Al2O3(0001) and dehydrated γ-Al2O3(100), dehydrated/hydroxylated γ-Al2O3(110), respectively. Structure models are illustrated using the open visualization tool (OVITO) package. (43) The threshold values for the Pd–Pd and Pd–Al bonds are 3.30 Å, and those for the Pd–O and Al–O bonds are 2.44 Å, respectively. The number of bonds and coordination numbers were calculated for adjacent atoms at distances below the threshold values. It should be noted that our Al2O3 surface model for DFT-AIMD calculations may not correspond to the true state of the Al2O3 surface under experimental conditions.
其中 T、σ 和 N 分别表示温度、标准差和自由原子数。该方程在 773 K 时产生以下标准差:Pd13/α-Al2O3(0001) 为 53.5,Pd13/γ-Al2O3(100) 为 53.5,Pd13/γ-Al2O3(110) 为 54.7。这些值与在 773 K 下进行的 5 ps 分子动力学(MD)模拟得出的温度标准差相当,Pd13/α-Al2O3(0001)为 52.6,Pd13/γ-Al2O3(100)为 51.7,Pd13/γ-Al2O3(110)为 51.3。从 AIMD 退火获得的初始结构开始,对 Pdn/Al2O3n = 1–4、6、9 和 13)进行结构优化。让 Al2O3 板和 Pdn 簇的顶部四个原子层(n = 1–4、6、9 和 13)松弛,其余底部原子固定。对于倒数空间的采样,分别对脱水α-Al2O3(0001)和脱γ水/羟基化γ-Al 2 O 3(110)和脱水/羟基化 -Al2O3(110)的(2 × 2 × 1)和(4 × 3 × ××1)的 Monkhorst-Pack k 点集进行采样。 结构模型使用开放可视化工具 (OVITO) 包进行说明。(43)Pd-Pd 和 Pd-Al 键的阈值分别为 3.30 Å,Pd-O 和 Al-O 键的阈值分别为 2.44 Å。计算距离低于阈值的相邻原子的键数和配位数。 需要注意的是,我们用于 DFT-AIMD 计算的 Al2O3 表面模型可能与实验条件下 Al2O3 表面的真实状态不符。

3. Results and Discussion

Click to copy section linkSection link copied!

3.1. Al2O3 Surface Models

Initial structures of dehydrated α-Al2O3(0001), γ-Al2O3(100), and γ-Al2O3(110) surfaces for DFT calculations and AIMD simulations were based on a previous report. (40) Previous studies have shown that the (110) surfaces are predominantly present on γ-Al2O3 (about 70% of all surfaces), followed by (100) surfaces (about 20%). (44,45) In the bulk of α-Al2O3, Al is present only at octahedral sites, whereas Al is present at octahedral (75%) and tetrahedral (25%) sites in the bulk of γ-Al2O3. In this study, Al2O3 surface models were obtained by AIMD simulation using the initial structures described above at 773 K (Movies S4–S6), followed by structural optimization at 0 K. The α-Al2O3(0001), γ-Al2O3(100), and γ-Al2O3(110) surface models are shown in Figure 1. The α-Al2O3(0001) surface had one coordinately unsaturated tricoordinated Al atom (AlIII) per unit cell, generated from the octahedral site in bulk. The γ-Al2O3(100) surface had four AlV sites per unit cell, generated from octahedral sites in bulk. The α-Al2O3(0001) and γ-Al2O3(100) surface models were consistent with previously reported surface structures. (40) The γ-Al2O3(110) surface model was more complex than γ-Al2O3(110) surface model without AIMD simulations (shown in Figure S3) due to the surface reconstruction during AIMD annealing at 773 K. The periodicity of (2 × 2) supercells for γ-Al2O3(110) surface was reduced to (2 × 1) supercells after the procedure of the AIMD annealing and structural optimization. The γ-Al2O3(110) surface displayed an uneven structure with six tetracoordinated Al3+ (AlIV) sites, a hexacoordinated Al3+ (AlVI) site, an AlV site, and two AlIII sites per surface unit cell. Such coordinatively unsaturated Al3+ sites function as Lewis acid sites. (4,5,22) Table 1 shows the density of the O and Al atoms on the topmost layers of Al2O3 (see Figure 1). The atom density of O on the Al2O3 surface increased in the order γ-Al2O3(110) < γ-Al2O3(100) < α-Al2O3(0001). In contrast, the atom density of Al was larger in the order α-Al2O3(0001) < γ-Al2O3(110) < γ-Al2O3(100). Despite the same crystalline phase, the O and Al atom density of the γ-Al2O3(110) surface was smaller than that of γ-Al2O3(100).

Figure 1

Figure 1. Top and side views of (a, b) α-Al2O3(0001), (c, d) γ-Al2O3(100), and (e, f) γ-Al2O3(110) surface. Only the topmost layers of the periodic slab are represented using color balls. A dashed line indicates the unit cell; the O and Al atoms are shown as red and brown balls.

Table 1. Density of the O and Al Atoms on the Topmost Layers of Al2O3
modeldensity of O atom (nm–2)density of Al atom (nm–2)
α-Al2O3(0001)15.25.1
γ-Al2O3(100)13.010.8
γ-Al2O3(110)10.47.5

3.2. Static Structure of Pd Clusters on Al2O3

For structural optimization of Pdn clusters (n = 1–4, 6, 9, 13) on the Al2O3 surfaces, AIMD simulations were employed to reduce artifacts manually placing Pdn clusters on the Al2O3. (41) AIMD simulations showed that the Pd clusters on the Al2O3 surfaces migrate to stable adsorption sites with structural changes. At least three structures of Pdn/Al2O3 with low energy in AIMD annealing held at 773 K for 4 ps were sampled as initial structures for the structural optimization at 0 K (Figure S4), (46) which resulted in different optimized structures. Among them, the optimized geometry with the lowest energy was adopted as the most stable structure. Figure 2 and Tables S1–S4 show the optimized structures of the Pdn clusters on the Al2O3 surface and those in a vacuum. The adsorption sites and geometric structure of the Pdn clusters (n = 1–4) on γ-Al2O3(100) in this work agreed with previous reports by Valero et al. (Table S2). (30) At cluster sizes of n ≥ 4, three-dimensional Pdn clusters were formed on γ-Al2O3(100). The point group can classify the shape of clusters with high symmetry. However, such a classification cannot be applied to clusters having no symmetry. Here, we adopt the molecular shape classification used in rotational spectra of polyatomic molecules. The diagonalized matrix of moments of inertia was used to classify the morphology of clusters. Three elements of the diagonal component are called the principal axes of inertia, which are generally assigned to IA, IB, and IC in descending order. Comparing the values of IA, IB, and IC leads to the systematic classification of molecular shapes. Thus, the shape of Pd13 clusters was quantitatively evaluated by the moment of inertia (Table 2). (47) Because the moment of inertia of a Pd13 cluster in a vacuum was IA < IB < IC, the Pd13 cluster displayed an asymmetric geometry. (46,48) A Pd13 cluster on γ-Al2O3(100) also showed an asymmetric geometry and consisted of three layers (Figures 2c,d). On the other hand, because the moment of inertia of a Pd13 cluster on α-Al2O3(0001) was IAIB < IC, the morphology of the Pd13 cluster was oblate, and its polar axis was perpendicular to the surface (Figures 2a,b). In contrast, the moment of inertia of a Pd13 cluster on γ-Al2O3(110) was IA < IBIC. Thus, the morphology of the Pd13 cluster on γ-Al2O3(110) was prolate and consisted of two layers, and its polar axis was parallel to the surface (Figures 2e,f).

Figure 2

Figure 2. Top and side views of (a, b) Pd13/α-Al2O3(0001), (c, d) Pd13/γ-Al2O3(100), and (e, f) Pd13/γ-Al2O3(110). O, Al, and Pd atoms are shown as red, brown, and navy-blue balls.

Table 2. Moment of Inertia of Pd13 Clusters in a Vacuum and on Al2O3 Surfaces
 moment of inertia 
modelIA (amu Å2)IB (amu Å2)IC (amu Å2)IA/IAIB/IAIC/IA morphology
Pd13 cluster5202653681751.01.31.6IA < IB < ICasymmetric top molecule
Pd13/α-Al2O3(0001)5695621582551.01.11.4IAIB < ICoblate
Pd13/γ-Al2O3(100)5220666285141.01.31.6IA < IB < ICasymmetric top molecule
Pd13/γ-Al2O3(110)46568348104851.01.82.3IA < IBICprolate

3.3. Interactions between Pd Clusters and Al2O3 Surfaces

The static geometry of the Pd clusters in the Pd/Al2O3 systems was strongly affected by the Pd–Al2O3 interactions originating from the surface structures. The Pd–Al2O3 interactions can be discussed based on the binding energy per Pd atom (Eb) and the interaction energy per Pd atom (Ei) shown in the following equations: (29,32)
Eb(for isolated Pdncluster)=1n(EPdnclusternEPd atom)
Eb(for Pdn/Al2O3)=1n(EPdn/Al2O3nEPd atomEAl2O3)
Ei=1n(EPdn/Al2O3EPdnclusterEAl2O3)
where EPdn/Al2O3, EPd atom, EAl2O3, and EPdn cluster are the total energies of Pdn/Al2O3, Pd atom, Al2O3, and Pdn clusters, respectively. Thus, the larger the negative values of these energetic parameters, the more stable the Pd clusters are on the Al2O3 surfaces. Figure 3a shows a plot of the binding energies, Eb, against the number of constituent atoms of Pd clusters on Al2O3. Eb of any Pdn clusters on the Al2O3 surfaces was more negative than those of the isolated Pdn clusters, indicating that the Pd clusters are stabilized by Al2O3 supports. As the size of the Pd clusters increased, Eb generally decreased. This means that the formation of larger Pd clusters is energetically more stable. When the size of the Pd clusters was small (n ≤ 4), Eb changed significantly with increasing the size of Pd clusters to strongly reflect the interactions of Pd with the local structure of the Al2O3. The Eb value of Pd1/γ-Al2O3(100) was more negative than those of Pd2/γ-Al2O3(100) and Pd3/γ-Al2O3(100). In the case of Pdn/γ-Al2O3(110), Eb had a local minimum when n = 2. These results are similar to the previous report for Pt/Al2O3, suggesting that the Pd clusters consisting of a few atoms were energetically stabilized on γ-Al2O3. (12) In the size region of n ≥ 4, Eb of Pdn/Al2O3 decreased gradually with increasing the size. As the size increased, Eb was more strongly affected by the surface structures of Al2O3 than by the cluster size. Comparing Eb of the Pd13 clusters on the different Al2O3 surfaces, they were in the order γ-Al2O3(110) < α-Al2O3(0001) < γ-Al2O3(100). This means that the Pd13 cluster was most stabilized on γ-Al2O3(110). The interactions between Pd clusters and Al2O3 surfaces reflected the broad region of the surface corresponding to the entire unit cell, rather than the local sites on Al2O3 surfaces. The interaction energy Ei for the Pd/Al2O3 systems contains the deformation of Pd clusters and Al2O3 surfaces by the adsorption of Pd clusters on the Al2O3 surfaces. Figure 3b shows a plot of the Ei of Pdn/Al2O3 versus Pd cluster size. In the size region of n ≤ 4, Ei increased rapidly as the size increased. In sizes of n > 4, all Ei values were slightly changed with respect to the cluster size, whereas the differences in Ei for the Al2O3 surfaces were varied. Ei of Pd/Al2O3 was in the order γ-Al2O3(110) < α-Al2O3(0001) < γ-Al2O3(100). Thus, the interactions between the Pd clusters and Al2O3 were stronger for Pd/γ-Al2O3(110) than those for Pd/γ-Al2O3(100) and Pd/α-Al2O3(0001). Furthermore, it is suggested that the suppression of the Pd13 shape on γ-Al2O3(110) comes from the increase in the interactions between Pd atoms and the Al2O3 surface. Such suppression on the other Al2O3 surfaces was less than that on the γ-Al2O3(110) surface.

Figure 3

Figure 3. Binding and interaction energies for isolated Pdn clusters and Pdn clusters (n = 1–4, 6, 9, and 13) supported on various Al2O3 surfaces as a function of the number of Pd atoms per Pd cluster.

3.4. Relation between Pd–Al2O3 Interactions and Surface Structure of Al2O3

To understand the relation between Pd–Al2O3 interactions and the surface structure of Al2O3, the interaction energies, Ei, were plotted against the CN (coordination number of bonds per Pd atom) of Pd–O and Pd–Al (Figure 4). Ei generally decreased as the coordination numbers, CN, of Pd–O and Pd–Al decreased. When the CN of Pd–O was less than 0.5 and the CN of Pd–Al was less than 1, Ei was linearly proportional to CN regardless of the Al2O3 surfaces. When CN was greater than the above values, the linear proportionality to CN was lost. The surface-to-surface variation increased because the contact area between Al2O3 and small clusters was restricted. Furthermore, because almost all of the atoms composing the small clusters were bound to the surface atoms, the interfacial interaction of the small clusters is strongly influenced by the local structure of the surface. As the cluster size increased, the Pd atoms on the interface could not be freely arranged due to the influence of other Pd atoms. As a result, CN decreased and the Pd–Al2O3 interactions are approximately determined by the number of Al and O atoms of Al2O3 bound to the Pd clusters. As in Table S5, CNs of Pd–O and Pd–Al at the same size of the Pd cluster were the largest for Pd13/γ-Al2O3(110), followed by Pd13/α-Al2O3(0001). That is to say, Pd/γ-Al2O3(110) could form more bonds between Pd clusters and Al2O3, and the strong Pd–Al2O3 interactions were induced. Here, one might infer that the Pd–Al2O3 interactions depend on the Al and O atom density of the top surface of Al2O3, which is directly involved in the bonding with Pd. However, the Ei values of Pd/γ-Al2O3(110) were smaller than those of Pd/α-Al2O3(0001) and Pd/γ-Al2O3(100) even though the atom density of Al and O in γ-Al2O3(110) was not so large (Table 1). The reason the interactions cannot be explained in terms of the atomic density was involved with the surface flexibility. During the AIMD annealing (Movies S4–S6), the AlIII sites and adjacent O atoms in γ-Al2O3(110) were flexibly oscillating. The γ-Al2O3(110) surface showed that the AlIII sites in γ-Al2O3(110) did not bond with the O atoms in the direction perpendicular to the surface. The α-Al2O3(0001) with AlIII sites was also a flexible structure. On the other hand, the flexibility of the γ-Al2O3(100) surface was less because the AlV sites on the γ-Al2O3(100) surface were bound to O atoms not only in the direction parallel but also in the direction perpendicular to the surface. Furthermore, to evaluate the flexibility of the Al2O3 surface, the displacement of atoms on the surface due to the adsorption of Pd clusters on Al2O3 was investigated (Table 3). In Table 3, the z-axis was the direction perpendicular to the surface, and only the coordinates of Al and O atoms bonded to Pd were used for the statistics. The absolute value of the average z-coordinates of the displacement was largest in the Pd13/γ-Al2O3(110). The z-coordinates of the AlIII sites in close proximity to the Pd13 cluster on the Pd13/γ-Al2O3(110) system were significantly changed upward. On the other hand, the z-coordinates of O atoms on Pd13/γ-Al2O3(110) were changed downward. Thus, the γ-Al2O3(110) surface was flexible to the adsorption of the Pd cluster. In contrast, for Pd13/α-Al2O3(0001), the average z-coordinates of displacement of the Al and O atoms were almost zero. However, they did not mean no change of the z-coordinates of Al atoms on the Pd13/α-Al2O3(0001). According to the maximum and minimum displacement, the z-coordinates of Al atoms in α-Al2O3(0001) fluctuated up (+0.45 Å) and down (−0.54 Å) due to the adsorption of the Pd13 cluster, resulting in an average displacement of zero. Although the α-Al2O3(0001) surface had AlIII sites, the absolute value of the maximum displacement for α-Al2O3(0001) was less than that for γ-Al2O3(110). Therefore, the α-Al2O3(0001) surface was less flexible than the γ-Al2O3(110) surface. The absolute values of average displacement of O and Al atoms on Pd13/γ-Al2O3(100) was the smallest, and the surface structure of γ-Al2O3(100) was stiff. The above results indicate that the Pd–Al2O3 interactions are influenced by not only the atomic density but also the surface flexibility of Al2O3. The increase in the interaction between Pd and Al2O3 surface induced by the surface flexibility also led to the shape suppression of the Pd13 cluster on γ-Al2O3(110). On the other hand, such suppression on the γ-Al2O3(100) was less than that on the other Al2O3 surfaces.

Figure 4

Figure 4. Interaction energies for Pdn clusters (n = 1–4, 6, 9, and 13) supported on various Al2O3 surfaces as a function of the coordination number of (a) Pd–O and (b) Pd–Al.

Table 3. Average, Minimum, and Maximum z-Coordinate Values of the Displacement of Al and O Atoms Bonded to Pd in the Pd13-Supported Al2O3 Surface from the Bare Al2O3 Surface
 Al atom (Å)O atom (Å)
modelavminmaxavminmax
Pd13/α-Al2O3(0001)0.05–0.540.450.01–0.090.06
Pd13/γ-Al2O3(100)0.050.000.12–0.05–0.140.05
Pd13/γ-Al2O3(110)0.21–0.160.72–0.11–0.360.05
In practical terms, coordinatively unsaturated Al3+ sites on the Al2O3 surface are hydrated, forming partially terminal/bridging hydroxyl groups. (49−52) The structures of metal clusters and isolated metal atoms on hydrated Al2O3 surfaces were investigated in previous studies. For example, for the γ-Al2O3(100), Rhn clusters (n ≤ 5) interacted more strongly on the nonhydrated surface than on the hydrated surface. (31) Also, on hydrated γ-Al2O3(110), Rhn clusters bonded strongly to hydroxyl groups. Pt13 and Pd13 clusters interacted more strongly with the nonhydrated γ-Al2O3(100) surface than with the hydrated γ-Al2O3(110) surface. (29) The morphology of Pt22 changed from a biplanar shape on pristine γ-Al2O3(110) to three-dimensional, if hydroxyl coverage on γ-Al2O3(110) is greater than 8.9 OH nm–2. (53) Immobilization of Ag species by terminal hydroxyl groups on the Al2O3 surface was reported. (20) Isolated Ag+ ions were reversibly dispersed by exchange with protons of terminal hydroxyl groups adjacent to strong Lewis acid sites (coordinatively unsaturated Al3+ sites) on Al2O3. (21) With the transfer of terminal hydroxyl groups and protons on the Al2O3 surfaces onto Pt clusters, the interaction between Pt and Al2O3 becomes stronger. (32) In this study, the effect of hydration of Al2O3 on the supported Pd clusters was also investigated. The hydration of γ-Al2O3(110) at 773 K is 5.9 OH nm–2, which was used as a model. (40) Table S6 shows the Pd clusters on the hydrated γ-Al2O3(110) surface obtained by AIMD simulations and structural optimization. During the AIMD annealing at 773 K, the movement of hydroxyl groups on the hydrated Al2O3 surface occurred, but no proton was transferred from the terminal hydroxyl groups onto the Pd cluster. The Ei values of Pdn/hydrated γ-Al2O3(110) (n = 1, 4, 13) were larger than those of Pdn/dehydrated γ-Al2O3(110) (Figure S5). This result was consistent with previous report that the adsorption of Pd clusters (n = 1–5) on the hydrated γ-Al2O3(110) surface is weaker than on the nonhydrated γ-Al2O3(110) surface. (30) The coordinatively unsaturated Al3+ sites on the Al2O3 surface were hydrated, resulting in weakening of Pd–Al2O3 interactions. Although Al2O3 surfaces are very complex in practice, it was unveiled that Al2O3 surfaces with coordinatively unsaturated Al3+ sites have stronger interactions with Pd.

3.5. Dynamic Structure of Pd Clusters on Al2O3

Understanding the dynamics of Pd clusters on Al2O3 surfaces is important to gain more insight into the stability of Pd clusters and the interactions between Pd and Al2O3. AIMD simulations were performed in which the Pd13/Al2O3 systems were held at 773 K for 8 ps (Movies S7–S9 and Figure S6) after a temperature increase from 0 to 773 K during 2 ps. Figure 5 shows the average distances dadj between only adjacent Pd atoms and dall between all Pd atoms. For 0–3 ps, dall of Pd13/Al2O3 decreased significantly with time; after 3 ps, the dall remained almost constant, suggesting that the Pd13 clusters on Al2O3 were almost in equilibrium. On the other hand, dadj showed little change with respect to time. Therefore, the Pd clusters underwent structural changes while keeping the distance between adjacent Pd atoms constant. Figure 6 shows the variation of the x-, y-, and z-coordinates of the centroid of the Pd13 clusters. The Pd13 cluster on γ-Al2O3(100) migrated while stretching perpendicular to the Al2O3 surface (Figures 6b and S3b and Movie S8). On the other hand, the Pd13 cluster on γ-Al2O3(110) was anchored on the surface while maintaining the prolate shape (Figures 6c and S3c and Movie S9). Moreover, the migration rates of the Pd13 clusters on the Al2O3 surfaces were calculated by dividing the variation of the centroid of the Pd13 cluster by time (Table 4). The migration rate of the Pd13 cluster on γ-Al2O3(110) was 0.31 Å/ps, which was smaller than those of the Pd13 clusters on γ-Al2O3(100) and α-Al2O3(0001). The γ-Al2O3(110) surface suppressed the migration of the Pd13 cluster more in comparison with the γ-Al2O3(100) and α-Al2O3(0001) surface. This indicated that a barrier for the migration of the Pd13 cluster on γ-Al2O3(110) was higher than that on γ-Al2O3(100) and α-Al2O3(0001). During the thermal deformation of the Pd13 cluster, the γ-Al2O3(110) surface could keep the Pd–Al and Pd–O bonds without migration of the Pd13 cluster because the Al2O3 surface deforms flexibly in response to the deformation. In contrast, the deformation of the Pd13 cluster on the γ-Al2O3(100) surface was accompanied by the breakage of the Pd–Al or Pd–O bonds as the surface structure was stiff. As a result, the entire Pd cluster had to move across the Al2O3 surface to form new bonds with other surface atoms to stabilize the Pd cluster. The migration rate of the Pd13 cluster was also suppressed on α-Al2O3(0001) compared to γ-Al2O3(100). Therefore, we showed that the design of the Al2O3 surface can control the shape and dispersion of Pd clusters on Al2O3. Because the AIMD simulation time for calculating the migration rate is very short (5 ps), one should be aware that the simulation does not provide a complete sampling of the diffusion path on the surface. Thus, the migration rate may not correspond to the real values. However, the above results, which indicate the dependence of the diffusion of Pd NPs on the Al2O3 surfaces, are supported by the movies and the calculated migration rates.

Figure 5

Figure 5. Variation of the dall and dadj of (a) Pd13/α-Al2O3(0001), (b) Pd13/γ-Al2O3(100), and (c) Pd13/γ-Al2O3(110) for the MD cycle.

Figure 6

Figure 6. Variation of the average position of Pd13 clusters in (a) Pd13/α-Al2O3(0001), (b) Pd13/γ-Al2O3(100), and (c) Pd13/γ-Al2O3(110) in the x-, y-, and z-axis direction for MD cycle (center). Top and side views of Pd13/Al2O3 at 3 and 8 ps (left and right). O, Al, and Pd atoms are shown as red, brown, and navy-blue balls.

Table 4. Average Distance between Adjacent Pd Atoms dadj, Average Distance between All Pd Atoms dall, and Migration Rate of Pd13 Cluster for Pd13/α-Al2O3(0001), Pd13/γ-Al2O3(100), and Pd13/γ-Al2O3(110)
modeldadj (Å)dall (Å)migration rate of Pd particle (Å/ps)
Pd13/α-Al2O3(0001)2.77 ± 0.034.12 ± 0.100.66
Pd13/γ-Al2O3(100)2.77 ± 0.044.24 ± 0.171.34
Pd13/γ-Al2O3(110)2.77 ± 0.034.00 ± 0.060.31

4. Conclusions

Click to copy section linkSection link copied!

The geometric structure and dynamics of Pd clusters on various Al2O3 surfaces were investigated by DFT calculations combined with AIMD simulations. The structure of a Pd13 cluster with low symmetry on the γ-Al2O3(100) was similar to that of a Pd13 cluster in a vacuum. On the other hand, Pd13 clusters on α-Al2O3(0001) and γ-Al2O3(110) were oblate and prolate, respectively. These Pd clusters deformed to bond more with the Al2O3 surface. Comparing the reconstruction of the Al2O3 surface when Pd clusters were adsorbed, the γ-Al2O3(110) surface showed a greater displacement of the Al and O atoms than did the other surfaces. The γ-Al2O3(110) surface with coordinatively unsaturated 3-coordinated Al3+ sites was found to be more flexible than the other surfaces. From the point of view of dynamics, Pd13 clusters on γ-Al2O3(100) diffused during AIMD annealing at 773 K, stretching and shrinking perpendicularly to the Al2O3 surface. In contrast, Pd13 clusters on γ-Al2O3(110) were anchored on Al2O3 while maintaining a prolate shape. Focusing on Al2O3, the γ-Al2O3(110) surface deformed flexibly as the Pd13 cluster deformed, so that surface Al and O atoms moved closer to the Pd cluster to form more bonds. Because the γ-Al2O3(100) surface was stiff, the deforming of the Pd cluster led to the migration of the entire Pd cluster across the Al2O3 surface to maintain the number of Pd–Al and Pd–O bonds per Pd atom. The interactions between Pd and Al2O3 were discussed in terms of binding energy Eb and interaction energy Ei per Pd atom. In the case of Pd13 on Al2O3, both Eb and Ei were stronger in the order γ-Al2O3(100) < α-Al2O3(0001) < γ-Al2O3(110). The interactions between Pd clusters and Al2O3 were larger as the number of Al and O atoms bonded to the Pd cluster increased at the interface. However, the Al and O atom density of the γ-Al2O3(110) surface was smaller than that of α-Al2O3(0001) and γ-Al2O3(100). The γ-Al2O3(110) surface bonded more with Pd clusters due to the structural flexibility of the 3-coordinated Al3+ sites and the adjacent O atoms, causing stronger Pd–Al2O3 interactions. Hydration of the Al2O3 surface weakened the Pd–Al2O3 interaction because hydroxyl adsorbates covered the coordinatively unsaturated Al3+ sites.

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcc.3c02483.

  • Figure S1: transition of the temperature and the potential energy of Pd13 cluster on Al2O3 surfaces for MD cycle; Figure S2: histograms of the temperature distribution for 5 ps after 3 ps of holding at 773 K; Figure S3: initial structure of γ-Al2O3(110) surface before AIMD simulation; Figure S4: sampled initial structures for the structural optimization; Table S1: molecular structures of Pd clusters supported on α-Al2O3(0001); Table S2: molecular structures of Pd clusters supported on γ-Al2O3(100); Table S3: molecular structures of Pd clusters supported on γ-Al2O3(110); Table S4: molecular structures of Pd clusters in a vacuum; Table S5: coordination number, average interatomic distance, binding energies, and interaction energies, for Pd clusters supported on Al2O3; Table S6: molecular structures of Pd clusters supported on hydrated γ-Al2O3(110); Figure S5: binding energies for Pd cluster supported on dehydrated and hydrated γ-Al2O3(110) surface as a function of the number of Pd atoms per Pd cluster; Figure S6: snapshots of MD simulations at 773 K for Pd13/Al2O3 (PDF)

  • Movie S1a: top view of Pd13/α-Al2O3(0001) (MP4)

  • Movie S1b: side view of Pd13/α-Al2O3(0001) (MP4)

  • Movie S2a: top view of Pd13/γ-Al2O3(100) (MP4)

  • Movie S2b: side view of Pd13/γ-Al2O3(100) (MP4)

  • Movie S3a: top view of Pd13/γ-Al2O3(110) (MP4)

  • Movie S3b: side view of Pd13/γ-Al2O3(110) (MP4)

  • Movie S4a: top view of α-Al2O3(0001) (MP4)

  • Movie S4b: side view of α-Al2O3(0001) (MP4)

  • Movie S5a: top view of γ-Al2O3(100) (MP4)

  • Movie S5b: side view of γ-Al2O3(100) (MP4)

  • Movie S6a: top view of γ-Al2O3(110) (MP4)

  • Movie S6b: side view of γ-Al2O3(110) (MP4)

  • Movie S7a: top view of Pd13/α-Al2O3(0001) (MP4)

  • Movie S7b: side view of Pd13/α-Al2O3(0001) (MP4)

  • Movie S8a: top view of Pd13/γ-Al2O3(100) (MP4)

  • Movie S8b: side view of Pd13/γ-Al2O3(100) (MP4)

  • Movie S9a: top view of Pd13/γ-Al2O3(110) (MP4)

  • Movie S9b: side view of Pd13/γ-Al2O3(110) (MP4)

Unveiling Static and Dynamic Structures of Pd Clusters Influenced by Al2O3 Surfaces: DFT and AIMD Studies

24 views

0 shares

0 downloads

S1
Supporting information
Unveiling Static and Dynamic Structures of Pd Clusters
Influenced by Al
2
O
3
Surfaces: DFT and AIMD Studies
Kazumasa Murata
, Akira Oda
, ‡
, Atsushi Satsuma
, ‡
, Kyoichi Sawabe*
, ‡
Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-
8520, Japan
Corresponding Author
* E-mail:
sawabe@chembio.nagoya-u.ac.jp
S2
Figure S1.
Transition of the temperature and potential energy of Pd
13
cluster on Al
2
O
3
surfaces during the simulated annealing
and holding by AIMD. (a) Pd
13
/
-Al
2
O
3
(0001), (b) Pd
13
/
-Al
2
O
3
(100), and (c) Pd
13
/
-Al
2
O
3
(110).
(
a
)
(
b
)
(
c
)
Figure S2.
Histograms of the temperature distribution for 5 ps after 3 ps of holding at 773K. (a) Pd
13
/
-Al
2
O
3
(0001), (b)
Pd
13
/
-Al
2
O
3
(100), and (c) Pd
13
/
-Al
2
O
3
(110).

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Author
  • Authors
    • Kazumasa Murata - Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
    • Akira Oda - Graduate School of Engineering, Nagoya University, Nagoya 464-8603, JapanElements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, JapanOrcidhttps://orcid.org/0000-0002-0217-4202
    • Atsushi Satsuma - Graduate School of Engineering, Nagoya University, Nagoya 464-8603, JapanElements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, JapanOrcidhttps://orcid.org/0000-0001-6946-7678
  • Notes
    The authors declare no competing financial interest.

Acknowledgments

Click to copy section linkSection link copied!

This work was partly supported by the JSPS KAKENHI Grant-in-Aids for Scientific Research (B) (Grant No. 18H01787) and JSPS Fellows (Grant No. 19J15440) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. A portion of this work was performed under the management of the Elements Strategy Initiative for Catalysts & Batteries (ESICB, JPMXP0112101003), which is also supported by MEXT. We used the supercomputer of ACCMS, Kyoto University.

Abbreviations

Click to copy section linkSection link copied!

MNPs

metal nanoparticles

MSIs

metal–support interactions

DFT

density functional theory

AIMD

ab initio molecular dynamics

CN

coordination number.

References

Click to copy section linkSection link copied!

This article references 53 other publications.

  1. 1
    Liu, L.; Corma, A. Metal Catalysts for Heterogeneous Catalysis: From Single Atoms to Nanoclusters and Nanoparticles. Chem. Rev. 2018, 118, 49815079,  DOI: 10.1021/acs.chemrev.7b00776
  2. 2
    van Deelen, T. W.; Hernández Mejía, C.; de Jong, K. P. Control of Metal-Support Interactions in Heterogeneous Catalysts to Enhance Activity and Selectivity. Nat. Catal. 2019, 2, 955970,  DOI: 10.1038/s41929-019-0364-x
  3. 3
    Omoda, T.; Takano, S.; Tsukuda, T. Toward Controlling the Electronic Structures of Chemically Modified Superatoms of Gold and Silver. Small 2021, 17, 118,  DOI: 10.1002/smll.202170136
  4. 4
    Wischert, R.; Laurent, P.; Copéret, C.; Delbecq, F.; Sautet, P. γ-Alumina: The Essential and Unexpected Role of Water for the Structure, Stability, and Reactivity of “Defect” Sites. J. Am. Chem. Soc. 2012, 134, 1443014449,  DOI: 10.1021/ja3042383
  5. 5
    Wischert, R.; Copéret, C.; Delbecq, F.; Sautet, P. Optimal Water Coverage on Alumina: A Key to Generate Lewis Acid-Base Pairs That Are Reactive towards the C-H Bond Activation of Methane. Angew. Chemie - Int. Ed. 2011, 50, 32023205,  DOI: 10.1002/anie.201006794
  6. 6
    Zhao, Z.; Xiao, D.; Chen, K.; Wang, R.; Liang, L.; Liu, Z.; Hung, I.; Gan, Z.; Hou, G. Nature of Five-Coordinated Al in Γ-Al2O3 Revealed by Ultra-High-Field Solid-State NMR. ACS Cent. Sci. 2022, 8, 795803,  DOI: 10.1021/acscentsci.1c01497
  7. 7
    Ayoola, H. O.; Bonifacio, C. S.; Zhu, Q.; Li, C. H.; House, S. D.; Kas, J. J.; Jinschek, J.; Rehr, J. J.; Saidi, W. A.; Yang, J. C. Probing the Local Bonding at the Pt/γ-Al2O3 Interface. J. Phys. Chem. C 2020, 124, 98769885,  DOI: 10.1021/acs.jpcc.9b12029
  8. 8
    Dessal, C.; Sangnier, A.; Chizallet, C.; Dujardin, C.; Morfin, F.; Rousset, J. L.; Aouine, M.; Bugnet, M.; Afanasiev, P.; Piccolo, L. Atmosphere-Dependent Stability and Mobility of Catalytic Pt Single Atoms and Clusters on γ-Al2O3. Nanoscale 2019, 11, 68976904,  DOI: 10.1039/C9NR01641D
  9. 9
    Wan, C.; Hu, M. Y.; Jaegers, N. R.; Shi, D.; Wang, H.; Gao, F.; Qin, Z.; Wang, Y.; Hu, J. Z. Investigating the Surface Structure of γ-Al2O3 Supported WOX Catalysts by High Field 27Al MAS NMR and Electronic Structure Calculations. J. Phys. Chem. C 2016, 120, 2309323103,  DOI: 10.1021/acs.jpcc.6b09060
  10. 10
    Kwak, J. H.; Hu, J. Z.; Kim, D. H.; Szanyi, J.; Peden, C. H. F. Penta-Coordinated Al3+ Ions as Preferential Nucleation Sites for BaO on Al2O3: An Ultra-High-Magnetic Field 27Al MAS NMR Study. J. Catal. 2007, 251, 189194,  DOI: 10.1016/j.jcat.2007.06.029
  11. 11
    Kwak, J. H.; Hu, J.; Mei, D.; Yi, C.-W.; Kim, D. H.; Peden, C. H. F.; Allard, L. F.; Szanyi, J. Coordinatively Unsaturated Al3+ Centers as Binding Sites for Active Catalyst Phases of Platinum on γ-Al2O3. Science 2009, 325, 16701673,  DOI: 10.1126/science.1176745
  12. 12
    Mei, D.; Kwak, J. H.; Hu, J.; Cho, S. J.; Szanyi, J.; Allard, L. F.; Peden, C. H. F. Unique Role of Anchoring Penta-Coordinated Al3+ Sites in the Sintering of γ-Al2O3-Supported Pt Catalysts. J. Phys. Chem. Lett. 2010, 1, 26882691,  DOI: 10.1021/jz101073p
  13. 13
    Zhang, Z.; Zhu, Y.; Asakura, H.; Zhang, B.; Zhang, J.; Zhou, M.; Han, Y.; Tanaka, T.; Wang, A.; Zhang, T.; Yan, N. Thermally Stable Single Atom Pt/m-Al2O3 for Selective Hydrogenation and CO Oxidation. Nat. Commun. 2017, 8, 16100,  DOI: 10.1038/ncomms16100
  14. 14
    Hackett, S. F. J.; Brydson, R. M.; Gass, M. H.; Harvey, I.; Newman, A. D.; Wilson, K.; Lee, A. F. High-Activity, Single-Site Mesoporous Pd/Al2O3 Catalysts for Selective Aerobic Oxidation of Allylic Alcohols. Angew. Chemie - Int. Ed. 2007, 46, 85938596,  DOI: 10.1002/anie.200702534
  15. 15
    Tang, N.; Cong, Y.; Shang, Q.; Wu, C.; Xu, G.; Wang, X. Coordinatively Unsaturated Al3+ Sites Anchored Subnanometric Ruthenium Catalyst for Hydrogenation of Aromatics. ACS Catal. 2017, 7, 5987,  DOI: 10.1021/acscatal.7b01816
  16. 16
    Li, L.; Wang, L. L.; Johnson, D. D.; Zhang, Z.; Sanchez, S. I.; Kang, J. H.; Nuzzo, R. G.; Wang, Q.; Frenkel, A. I.; Li, J. Noncrystalline-to-Crystalline Transformations in Pt Nanoparticles. J. Am. Chem. Soc. 2013, 135, 1306213072,  DOI: 10.1021/ja405497p
  17. 17
    Murata, K.; Onoda, J.; Yamamoto, Y.; Oda, A.; Ohyama, J.; Satsuma, A. Enhancement of Toluene Hydrogenation Activity of Supported Pt Nanoparticles with Increasing the Crystallinity of Pt. Appl. Catal. A, Gen. 2022, 629, 118425,  DOI: 10.1016/j.apcata.2021.118425
  18. 18
    Batista, A. T. F.; Baaziz, W.; Taleb, A.; Chaniot, J.; Moreaud, M.; Aguilar-tapia, A.; Proux, O.; Hazemann, J.; Diehl, F.; Gay, A. Atomic Scale Insight into the Formation, Size, and Location of Platinum Nanoparticles Supported on γ-Alumina. ACS Catal. 2020, 10, 41934204,  DOI: 10.1021/acscatal.0c00042
  19. 19
    Vaarkamp, M.; Miller, J. T.; Modica, F. S.; Koningsberger, D. C. On the Relation between Particle Morphology, Structure of the Metal-Support Interface, and Catalytic Properties of Pt/γ-Al2O3. J. Catal. 1996, 163, 294305,  DOI: 10.1006/jcat.1996.0330
  20. 20
    Wang, F.; Ma, J.; Xin, S.; Wang, Q.; Xu, J.; Zhang, C.; He, H.; Cheng Zeng, X. Resolving the Puzzle of Single-Atom Silver Dispersion on Nanosized γ-Al2O3 Surface for High Catalytic Performance. Nat. Commun. 2020, 11, 529,  DOI: 10.1038/s41467-019-13937-1
  21. 21
    Kubota, H.; Mine, S.; Toyao, T.; Maeno, Z.; Shimizu, K. Redox-Driven Reversible Structural Evolution of Isolated Silver Atoms Anchored to Specific Sites on γ-Al2O3. ACS Catal. 2022, 12, 544599,  DOI: 10.1021/acscatal.1c04924
  22. 22
    Valero, M. C.; Raybaud, P.; Sautet, P. Influence of the Hydroxylation of γ-Al2O3 Surfaces on the Stability and Diffusion of Single Pd Atoms: A DFT Study. J. Phys. Chem. B 2006, 110, 17591767,  DOI: 10.1021/jp0554240
  23. 23
    Murata, K.; Mahara, Y.; Ohyama, J.; Yamamoto, Y.; Arai, S.; Satsuma, A. The Metal-Support Interaction Concerning the Particle Size Effect of Pd/Al2O3 on Methane Combustion. Angew. Chemie Int. Ed. 2017, 56, 1599315997,  DOI: 10.1002/anie.201709124
  24. 24
    Murata, K.; Shiotani, T.; Ohyama, J.; Satsuma, A. Selective Hydrogenation of C═C Bond in Cinnamaldehyde on Pd Step Sites of Pd/Al2O3. Chem. Lett. 2021, 50, 599602,  DOI: 10.1246/cl.200856
  25. 25
    Murata, K.; Eleeda, E.; Ohyama, J.; Yamamoto, Y.; Arai, S.; Satsuma, A. Identification of Active Sites in CO Oxidation over a Pd/Al2O3 Catalyst. Phys. Chem. Chem. Phys. 2019, 21, 1812818137,  DOI: 10.1039/C9CP03943K
  26. 26
    Duan, H.; You, R.; Xu, S.; Li, Z.; Qian, K.; Cao, T.; Huang, W.; Bao, X. Penta-Coordinated Al3+ Stabilized Active Pd Structures on Al2O3 Coated Palladium Catalysts for Methane Combustion. Angew. Chemie Int. Ed. 2019, 58, 1204312048,  DOI: 10.1002/anie.201904883
  27. 27
    Yang, X.; Li, Q.; Lu, E.; Wang, Z.; Gong, X.; Yu, Z.; Guo, Y.; Wang, L.; Guo, Y.; Zhan, W.; Zhang, J.; Dai, S. Taming the Stability of Pd Active Phases through a Compartmentalizing Strategy toward Nanostructured Catalyst Supports. Nat. Commun. 2019, 10, 1611,  DOI: 10.1038/s41467-019-09662-4
  28. 28
    Murata, K.; Shiotani, T.; Ohyama, J.; Wakabayashi, R.; Maruoka, H.; Kimura, T.; Satsuma, A. Relationship between Penta-Coordinated Al3+ Sites in the Al2O3 Supports and CH4 Combustion Activity of Pd/Al2O3 Catalysts. Catal. Sci. Technol. 2021, 11, 23742378,  DOI: 10.1039/D1CY00098E
  29. 29
    Hu, C. H.; Chizallet, C.; Mager-Maury, C.; Corral-Valero, M.; Sautet, P.; Toulhoat, H.; Raybaud, P. Modulation of Catalyst Particle Structure upon Support Hydroxylation: Ab Initio Insights into Pd13 and Pt13/γ-Al2O3. J. Catal. 2010, 274, 99110,  DOI: 10.1016/j.jcat.2010.06.009
  30. 30
    Valero, M. C.; Raybaud, P.; Sautet, P. Nucleation of Pdn (N = 1–5) Clusters and Wetting of Pd Particles on γ-Al2O3 Surfaces: A Density Functional Theory Study. Phys. Rev. B 2007, 75, 112,  DOI: 10.1103/PhysRevB.75.045427
  31. 31
    Shi, X. R.; Sholl, D. S. Nucleation of Rhn (n = 1–5) Clusters on γ-Al2O3 Surfaces: A Density Functional Theory Study. J. Phys. Chem. C 2012, 116, 1062310631,  DOI: 10.1021/jp301114n
  32. 32
    Mager-Maury, C.; Chizallet, C.; Sautet, P.; Raybaud, P. Platinum Nanoclusters Stabilized on γ-Alumina by Chlorine Used As a Capping Surface Ligand: A Density Functional Theory Study. ACS Catal. 2012, 2, 13461357,  DOI: 10.1021/cs300178y
  33. 33
    Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 1550,  DOI: 10.1016/0927-0256(96)00008-0
  34. 34
    Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 1116911186,  DOI: 10.1103/PhysRevB.54.11169
  35. 35
    Kresse, G.; Hafner, J. Ab. Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558561,  DOI: 10.1103/PhysRevB.47.558
  36. 36
    Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 1795317979,  DOI: 10.1103/PhysRevB.50.17953
  37. 37
    Hammer, B.; Hansen, L. B.; Nørskov, J. K. Improved Adsorption Energetics within Density-Functional Theory Using Revised Perdew-Burke-Ernzerhof Functionals. Phys. Rev. B - Condens. Matter Mater. Phys. 1999, 59, 74137421,  DOI: 10.1103/PhysRevB.59.7413
  38. 38
    Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104,  DOI: 10.1063/1.3382344
  39. 39
    Eng, P. J.; Trainor, T. P.; Brown, G. E., Jr.; Waychunas, G. A.; Newville, M.; Sutton, S. R.; Rivers, M. L. Structure of the Hydrated α-Al2O3(0001) Surface. Science 2000, 288, 10291033,  DOI: 10.1126/science.288.5468.1029
  40. 40
    Digne, M.; Sautet, P.; Raybaud, P.; Euzen, P.; Toulhoat, H. Use of DFT to Achieve a Rational Understanding of Acid-Basic Properties of γ-Alumina Surfaces. J. Catal. 2004, 226, 5468,  DOI: 10.1016/j.jcat.2004.04.020
  41. 41
    Sawabe, K.; Satsuma, A. Theoretical Study on Carbon Monoxide Adsorption on Unsupported and γ-Al2O3-Supported Silver Nanoparticles: Size, Shape, and Support Effects. ACS Omega 2022, 7, 44054412,  DOI: 10.1021/acsomega.1c06208
  42. 42
    Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications; Academic press: 2002, p 141
  43. 43
    Stukowski, A. Visualization and Analysis of Atomistic Simulation Data with OVITO-the Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012,  DOI: 10.1088/0965-0393/18/1/015012
  44. 44
    Digne, M.; Sautet, P.; Raybaud, P.; Euzen, P.; Toulhoat, H. Hydroxyl Groups on γ-Alumina Surfaces: A DFT Study. J. Catal. 2002, 211, 15,  DOI: 10.1006/jcat.2002.3741
  45. 45
    Nortier, P.; Fourre, P.; Saad, A.B. M.; Saur, O.; Lavalley, J.C. Effects of Crystallinity and Morphology on the Surface Properties of Alumina. Appl. Catal. 1990, 61, 141160,  DOI: 10.1016/S0166-9834(00)82140-5
  46. 46
    Wang, L. L.; Johnson, D. D. Density Functional Study of Structural Trends for Late-Transition-Metal 13-Atom Clusters. Phys. Rev. B 2007, 75, 110,  DOI: 10.1103/PhysRevB.75.235405
  47. 47
    Atkins, P.; de Paula, J. Atkins’ Physical Chemistry, 8th ed.; Oxford University Press: 2006; pp 441442.
  48. 48
    Zhang, W.; Xiao, L.; Hirata, Y.; Pawluk, T.; Wang, L. The Simple Cubic Structure of Ir Clusters and the Element Effect on Cluster Structures. Chem. Phys. Lett. 2004, 383, 6771,  DOI: 10.1016/j.cplett.2003.11.005
  49. 49
    Hendriksen, B. A.; Pearce, D. R.; Rudham, R. Heats of Adsorption of Water on α- and γ-Alumina. J. Catal. 1972, 24, 8287,  DOI: 10.1016/0021-9517(72)90010-3
  50. 50
    Thissen, P.; Grundmeier, G.; Wippermann, S.; Schmidt, W. G. Water Adsorption on the α-Al2O3 (0001) Surface. Phys. Rev. B - Condens. Matter Mater. Phys. 2009, 80, 16,  DOI: 10.1103/PhysRevB.80.245403
  51. 51
    Pigeon, T.; Chizallet, C.; Raybaud, P. Revisiting γ-Alumina Surface Models through the Topotactic Transformation of Boehmite Surfaces. J. Catal. 2022, 405, 140151,  DOI: 10.1016/j.jcat.2021.11.011
  52. 52
    Ranea, V. A.; Schneider, W. F.; Carmichael, I. DFT Characterization of Coverage Dependent Molecular Water Adsorption Modes on α-Al2O3(0001). Surf. Sci. 2008, 602, 268275,  DOI: 10.1016/j.susc.2007.10.029
  53. 53
    Erfani, G. S.; Hong, S.; Rahman, T. S. Effects of γ-Al2O3 Support on the Morphology and Electronic Structure of Pt Nanoparticles. J. Phys. Chem. C 2019, 123, 1689316901,  DOI: 10.1021/acs.jpcc.9b04770

Cited By

Click to copy section linkSection link copied!

This article has not yet been cited by other publications.

The Journal of Physical Chemistry C

Cite this: J. Phys. Chem. C 2023, 127, 41, 20267–20275
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcc.3c02483
Published October 5, 2023
Copyright © 2023 American Chemical Society

Article Views

1510

Altmetric

-

Citations

-
Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Top and side views of (a, b) α-Al2O3(0001), (c, d) γ-Al2O3(100), and (e, f) γ-Al2O3(110) surface. Only the topmost layers of the periodic slab are represented using color balls. A dashed line indicates the unit cell; the O and Al atoms are shown as red and brown balls.

    Figure 2

    Figure 2. Top and side views of (a, b) Pd13/α-Al2O3(0001), (c, d) Pd13/γ-Al2O3(100), and (e, f) Pd13/γ-Al2O3(110). O, Al, and Pd atoms are shown as red, brown, and navy-blue balls.

    Figure 3

    Figure 3. Binding and interaction energies for isolated Pdn clusters and Pdn clusters (n = 1–4, 6, 9, and 13) supported on various Al2O3 surfaces as a function of the number of Pd atoms per Pd cluster.

    Figure 4

    Figure 4. Interaction energies for Pdn clusters (n = 1–4, 6, 9, and 13) supported on various Al2O3 surfaces as a function of the coordination number of (a) Pd–O and (b) Pd–Al.

    Figure 5

    Figure 5. Variation of the dall and dadj of (a) Pd13/α-Al2O3(0001), (b) Pd13/γ-Al2O3(100), and (c) Pd13/γ-Al2O3(110) for the MD cycle.

    Figure 6

    Figure 6. Variation of the average position of Pd13 clusters in (a) Pd13/α-Al2O3(0001), (b) Pd13/γ-Al2O3(100), and (c) Pd13/γ-Al2O3(110) in the x-, y-, and z-axis direction for MD cycle (center). Top and side views of Pd13/Al2O3 at 3 and 8 ps (left and right). O, Al, and Pd atoms are shown as red, brown, and navy-blue balls.

  • References


    This article references 53 other publications.

    1. 1
      Liu, L.; Corma, A. Metal Catalysts for Heterogeneous Catalysis: From Single Atoms to Nanoclusters and Nanoparticles. Chem. Rev. 2018, 118, 49815079,  DOI: 10.1021/acs.chemrev.7b00776
    2. 2
      van Deelen, T. W.; Hernández Mejía, C.; de Jong, K. P. Control of Metal-Support Interactions in Heterogeneous Catalysts to Enhance Activity and Selectivity. Nat. Catal. 2019, 2, 955970,  DOI: 10.1038/s41929-019-0364-x
    3. 3
      Omoda, T.; Takano, S.; Tsukuda, T. Toward Controlling the Electronic Structures of Chemically Modified Superatoms of Gold and Silver. Small 2021, 17, 118,  DOI: 10.1002/smll.202170136
    4. 4
      Wischert, R.; Laurent, P.; Copéret, C.; Delbecq, F.; Sautet, P. γ-Alumina: The Essential and Unexpected Role of Water for the Structure, Stability, and Reactivity of “Defect” Sites. J. Am. Chem. Soc. 2012, 134, 1443014449,  DOI: 10.1021/ja3042383
    5. 5
      Wischert, R.; Copéret, C.; Delbecq, F.; Sautet, P. Optimal Water Coverage on Alumina: A Key to Generate Lewis Acid-Base Pairs That Are Reactive towards the C-H Bond Activation of Methane. Angew. Chemie - Int. Ed. 2011, 50, 32023205,  DOI: 10.1002/anie.201006794
    6. 6
      Zhao, Z.; Xiao, D.; Chen, K.; Wang, R.; Liang, L.; Liu, Z.; Hung, I.; Gan, Z.; Hou, G. Nature of Five-Coordinated Al in Γ-Al2O3 Revealed by Ultra-High-Field Solid-State NMR. ACS Cent. Sci. 2022, 8, 795803,  DOI: 10.1021/acscentsci.1c01497
    7. 7
      Ayoola, H. O.; Bonifacio, C. S.; Zhu, Q.; Li, C. H.; House, S. D.; Kas, J. J.; Jinschek, J.; Rehr, J. J.; Saidi, W. A.; Yang, J. C. Probing the Local Bonding at the Pt/γ-Al2O3 Interface. J. Phys. Chem. C 2020, 124, 98769885,  DOI: 10.1021/acs.jpcc.9b12029
    8. 8
      Dessal, C.; Sangnier, A.; Chizallet, C.; Dujardin, C.; Morfin, F.; Rousset, J. L.; Aouine, M.; Bugnet, M.; Afanasiev, P.; Piccolo, L. Atmosphere-Dependent Stability and Mobility of Catalytic Pt Single Atoms and Clusters on γ-Al2O3. Nanoscale 2019, 11, 68976904,  DOI: 10.1039/C9NR01641D
    9. 9
      Wan, C.; Hu, M. Y.; Jaegers, N. R.; Shi, D.; Wang, H.; Gao, F.; Qin, Z.; Wang, Y.; Hu, J. Z. Investigating the Surface Structure of γ-Al2O3 Supported WOX Catalysts by High Field 27Al MAS NMR and Electronic Structure Calculations. J. Phys. Chem. C 2016, 120, 2309323103,  DOI: 10.1021/acs.jpcc.6b09060
    10. 10
      Kwak, J. H.; Hu, J. Z.; Kim, D. H.; Szanyi, J.; Peden, C. H. F. Penta-Coordinated Al3+ Ions as Preferential Nucleation Sites for BaO on Al2O3: An Ultra-High-Magnetic Field 27Al MAS NMR Study. J. Catal. 2007, 251, 189194,  DOI: 10.1016/j.jcat.2007.06.029
    11. 11
      Kwak, J. H.; Hu, J.; Mei, D.; Yi, C.-W.; Kim, D. H.; Peden, C. H. F.; Allard, L. F.; Szanyi, J. Coordinatively Unsaturated Al3+ Centers as Binding Sites for Active Catalyst Phases of Platinum on γ-Al2O3. Science 2009, 325, 16701673,  DOI: 10.1126/science.1176745
    12. 12
      Mei, D.; Kwak, J. H.; Hu, J.; Cho, S. J.; Szanyi, J.; Allard, L. F.; Peden, C. H. F. Unique Role of Anchoring Penta-Coordinated Al3+ Sites in the Sintering of γ-Al2O3-Supported Pt Catalysts. J. Phys. Chem. Lett. 2010, 1, 26882691,  DOI: 10.1021/jz101073p
    13. 13
      Zhang, Z.; Zhu, Y.; Asakura, H.; Zhang, B.; Zhang, J.; Zhou, M.; Han, Y.; Tanaka, T.; Wang, A.; Zhang, T.; Yan, N. Thermally Stable Single Atom Pt/m-Al2O3 for Selective Hydrogenation and CO Oxidation. Nat. Commun. 2017, 8, 16100,  DOI: 10.1038/ncomms16100
    14. 14
      Hackett, S. F. J.; Brydson, R. M.; Gass, M. H.; Harvey, I.; Newman, A. D.; Wilson, K.; Lee, A. F. High-Activity, Single-Site Mesoporous Pd/Al2O3 Catalysts for Selective Aerobic Oxidation of Allylic Alcohols. Angew. Chemie - Int. Ed. 2007, 46, 85938596,  DOI: 10.1002/anie.200702534
    15. 15
      Tang, N.; Cong, Y.; Shang, Q.; Wu, C.; Xu, G.; Wang, X. Coordinatively Unsaturated Al3+ Sites Anchored Subnanometric Ruthenium Catalyst for Hydrogenation of Aromatics. ACS Catal. 2017, 7, 5987,  DOI: 10.1021/acscatal.7b01816
    16. 16
      Li, L.; Wang, L. L.; Johnson, D. D.; Zhang, Z.; Sanchez, S. I.; Kang, J. H.; Nuzzo, R. G.; Wang, Q.; Frenkel, A. I.; Li, J. Noncrystalline-to-Crystalline Transformations in Pt Nanoparticles. J. Am. Chem. Soc. 2013, 135, 1306213072,  DOI: 10.1021/ja405497p
    17. 17
      Murata, K.; Onoda, J.; Yamamoto, Y.; Oda, A.; Ohyama, J.; Satsuma, A. Enhancement of Toluene Hydrogenation Activity of Supported Pt Nanoparticles with Increasing the Crystallinity of Pt. Appl. Catal. A, Gen. 2022, 629, 118425,  DOI: 10.1016/j.apcata.2021.118425
    18. 18
      Batista, A. T. F.; Baaziz, W.; Taleb, A.; Chaniot, J.; Moreaud, M.; Aguilar-tapia, A.; Proux, O.; Hazemann, J.; Diehl, F.; Gay, A. Atomic Scale Insight into the Formation, Size, and Location of Platinum Nanoparticles Supported on γ-Alumina. ACS Catal. 2020, 10, 41934204,  DOI: 10.1021/acscatal.0c00042
    19. 19
      Vaarkamp, M.; Miller, J. T.; Modica, F. S.; Koningsberger, D. C. On the Relation between Particle Morphology, Structure of the Metal-Support Interface, and Catalytic Properties of Pt/γ-Al2O3. J. Catal. 1996, 163, 294305,  DOI: 10.1006/jcat.1996.0330
    20. 20
      Wang, F.; Ma, J.; Xin, S.; Wang, Q.; Xu, J.; Zhang, C.; He, H.; Cheng Zeng, X. Resolving the Puzzle of Single-Atom Silver Dispersion on Nanosized γ-Al2O3 Surface for High Catalytic Performance. Nat. Commun. 2020, 11, 529,  DOI: 10.1038/s41467-019-13937-1
    21. 21
      Kubota, H.; Mine, S.; Toyao, T.; Maeno, Z.; Shimizu, K. Redox-Driven Reversible Structural Evolution of Isolated Silver Atoms Anchored to Specific Sites on γ-Al2O3. ACS Catal. 2022, 12, 544599,  DOI: 10.1021/acscatal.1c04924
    22. 22
      Valero, M. C.; Raybaud, P.; Sautet, P. Influence of the Hydroxylation of γ-Al2O3 Surfaces on the Stability and Diffusion of Single Pd Atoms: A DFT Study. J. Phys. Chem. B 2006, 110, 17591767,  DOI: 10.1021/jp0554240
    23. 23
      Murata, K.; Mahara, Y.; Ohyama, J.; Yamamoto, Y.; Arai, S.; Satsuma, A. The Metal-Support Interaction Concerning the Particle Size Effect of Pd/Al2O3 on Methane Combustion. Angew. Chemie Int. Ed. 2017, 56, 1599315997,  DOI: 10.1002/anie.201709124
    24. 24
      Murata, K.; Shiotani, T.; Ohyama, J.; Satsuma, A. Selective Hydrogenation of C═C Bond in Cinnamaldehyde on Pd Step Sites of Pd/Al2O3. Chem. Lett. 2021, 50, 599602,  DOI: 10.1246/cl.200856
    25. 25
      Murata, K.; Eleeda, E.; Ohyama, J.; Yamamoto, Y.; Arai, S.; Satsuma, A. Identification of Active Sites in CO Oxidation over a Pd/Al2O3 Catalyst. Phys. Chem. Chem. Phys. 2019, 21, 1812818137,  DOI: 10.1039/C9CP03943K
    26. 26
      Duan, H.; You, R.; Xu, S.; Li, Z.; Qian, K.; Cao, T.; Huang, W.; Bao, X. Penta-Coordinated Al3+ Stabilized Active Pd Structures on Al2O3 Coated Palladium Catalysts for Methane Combustion. Angew. Chemie Int. Ed. 2019, 58, 1204312048,  DOI: 10.1002/anie.201904883
    27. 27
      Yang, X.; Li, Q.; Lu, E.; Wang, Z.; Gong, X.; Yu, Z.; Guo, Y.; Wang, L.; Guo, Y.; Zhan, W.; Zhang, J.; Dai, S. Taming the Stability of Pd Active Phases through a Compartmentalizing Strategy toward Nanostructured Catalyst Supports. Nat. Commun. 2019, 10, 1611,  DOI: 10.1038/s41467-019-09662-4
    28. 28
      Murata, K.; Shiotani, T.; Ohyama, J.; Wakabayashi, R.; Maruoka, H.; Kimura, T.; Satsuma, A. Relationship between Penta-Coordinated Al3+ Sites in the Al2O3 Supports and CH4 Combustion Activity of Pd/Al2O3 Catalysts. Catal. Sci. Technol. 2021, 11, 23742378,  DOI: 10.1039/D1CY00098E
    29. 29
      Hu, C. H.; Chizallet, C.; Mager-Maury, C.; Corral-Valero, M.; Sautet, P.; Toulhoat, H.; Raybaud, P. Modulation of Catalyst Particle Structure upon Support Hydroxylation: Ab Initio Insights into Pd13 and Pt13/γ-Al2O3. J. Catal. 2010, 274, 99110,  DOI: 10.1016/j.jcat.2010.06.009
    30. 30
      Valero, M. C.; Raybaud, P.; Sautet, P. Nucleation of Pdn (N = 1–5) Clusters and Wetting of Pd Particles on γ-Al2O3 Surfaces: A Density Functional Theory Study. Phys. Rev. B 2007, 75, 112,  DOI: 10.1103/PhysRevB.75.045427
    31. 31
      Shi, X. R.; Sholl, D. S. Nucleation of Rhn (n = 1–5) Clusters on γ-Al2O3 Surfaces: A Density Functional Theory Study. J. Phys. Chem. C 2012, 116, 1062310631,  DOI: 10.1021/jp301114n
    32. 32
      Mager-Maury, C.; Chizallet, C.; Sautet, P.; Raybaud, P. Platinum Nanoclusters Stabilized on γ-Alumina by Chlorine Used As a Capping Surface Ligand: A Density Functional Theory Study. ACS Catal. 2012, 2, 13461357,  DOI: 10.1021/cs300178y
    33. 33
      Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 1550,  DOI: 10.1016/0927-0256(96)00008-0
    34. 34
      Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 1116911186,  DOI: 10.1103/PhysRevB.54.11169
    35. 35
      Kresse, G.; Hafner, J. Ab. Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558561,  DOI: 10.1103/PhysRevB.47.558
    36. 36
      Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 1795317979,  DOI: 10.1103/PhysRevB.50.17953
    37. 37
      Hammer, B.; Hansen, L. B.; Nørskov, J. K. Improved Adsorption Energetics within Density-Functional Theory Using Revised Perdew-Burke-Ernzerhof Functionals. Phys. Rev. B - Condens. Matter Mater. Phys. 1999, 59, 74137421,  DOI: 10.1103/PhysRevB.59.7413
    38. 38
      Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104,  DOI: 10.1063/1.3382344
    39. 39
      Eng, P. J.; Trainor, T. P.; Brown, G. E., Jr.; Waychunas, G. A.; Newville, M.; Sutton, S. R.; Rivers, M. L. Structure of the Hydrated α-Al2O3(0001) Surface. Science 2000, 288, 10291033,  DOI: 10.1126/science.288.5468.1029
    40. 40
      Digne, M.; Sautet, P.; Raybaud, P.; Euzen, P.; Toulhoat, H. Use of DFT to Achieve a Rational Understanding of Acid-Basic Properties of γ-Alumina Surfaces. J. Catal. 2004, 226, 5468,  DOI: 10.1016/j.jcat.2004.04.020
    41. 41
      Sawabe, K.; Satsuma, A. Theoretical Study on Carbon Monoxide Adsorption on Unsupported and γ-Al2O3-Supported Silver Nanoparticles: Size, Shape, and Support Effects. ACS Omega 2022, 7, 44054412,  DOI: 10.1021/acsomega.1c06208
    42. 42
      Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications; Academic press: 2002, p 141
    43. 43
      Stukowski, A. Visualization and Analysis of Atomistic Simulation Data with OVITO-the Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012,  DOI: 10.1088/0965-0393/18/1/015012
    44. 44
      Digne, M.; Sautet, P.; Raybaud, P.; Euzen, P.; Toulhoat, H. Hydroxyl Groups on γ-Alumina Surfaces: A DFT Study. J. Catal. 2002, 211, 15,  DOI: 10.1006/jcat.2002.3741
    45. 45
      Nortier, P.; Fourre, P.; Saad, A.B. M.; Saur, O.; Lavalley, J.C. Effects of Crystallinity and Morphology on the Surface Properties of Alumina. Appl. Catal. 1990, 61, 141160,  DOI: 10.1016/S0166-9834(00)82140-5
    46. 46
      Wang, L. L.; Johnson, D. D. Density Functional Study of Structural Trends for Late-Transition-Metal 13-Atom Clusters. Phys. Rev. B 2007, 75, 110,  DOI: 10.1103/PhysRevB.75.235405
    47. 47
      Atkins, P.; de Paula, J. Atkins’ Physical Chemistry, 8th ed.; Oxford University Press: 2006; pp 441442.
    48. 48
      Zhang, W.; Xiao, L.; Hirata, Y.; Pawluk, T.; Wang, L. The Simple Cubic Structure of Ir Clusters and the Element Effect on Cluster Structures. Chem. Phys. Lett. 2004, 383, 6771,  DOI: 10.1016/j.cplett.2003.11.005
    49. 49
      Hendriksen, B. A.; Pearce, D. R.; Rudham, R. Heats of Adsorption of Water on α- and γ-Alumina. J. Catal. 1972, 24, 8287,  DOI: 10.1016/0021-9517(72)90010-3
    50. 50
      Thissen, P.; Grundmeier, G.; Wippermann, S.; Schmidt, W. G. Water Adsorption on the α-Al2O3 (0001) Surface. Phys. Rev. B - Condens. Matter Mater. Phys. 2009, 80, 16,  DOI: 10.1103/PhysRevB.80.245403
    51. 51
      Pigeon, T.; Chizallet, C.; Raybaud, P. Revisiting γ-Alumina Surface Models through the Topotactic Transformation of Boehmite Surfaces. J. Catal. 2022, 405, 140151,  DOI: 10.1016/j.jcat.2021.11.011
    52. 52
      Ranea, V. A.; Schneider, W. F.; Carmichael, I. DFT Characterization of Coverage Dependent Molecular Water Adsorption Modes on α-Al2O3(0001). Surf. Sci. 2008, 602, 268275,  DOI: 10.1016/j.susc.2007.10.029
    53. 53
      Erfani, G. S.; Hong, S.; Rahman, T. S. Effects of γ-Al2O3 Support on the Morphology and Electronic Structure of Pt Nanoparticles. J. Phys. Chem. C 2019, 123, 1689316901,  DOI: 10.1021/acs.jpcc.9b04770
  • Supporting Information

    Supporting Information


    The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcc.3c02483.

    • Figure S1: transition of the temperature and the potential energy of Pd13 cluster on Al2O3 surfaces for MD cycle; Figure S2: histograms of the temperature distribution for 5 ps after 3 ps of holding at 773 K; Figure S3: initial structure of γ-Al2O3(110) surface before AIMD simulation; Figure S4: sampled initial structures for the structural optimization; Table S1: molecular structures of Pd clusters supported on α-Al2O3(0001); Table S2: molecular structures of Pd clusters supported on γ-Al2O3(100); Table S3: molecular structures of Pd clusters supported on γ-Al2O3(110); Table S4: molecular structures of Pd clusters in a vacuum; Table S5: coordination number, average interatomic distance, binding energies, and interaction energies, for Pd clusters supported on Al2O3; Table S6: molecular structures of Pd clusters supported on hydrated γ-Al2O3(110); Figure S5: binding energies for Pd cluster supported on dehydrated and hydrated γ-Al2O3(110) surface as a function of the number of Pd atoms per Pd cluster; Figure S6: snapshots of MD simulations at 773 K for Pd13/Al2O3 (PDF)

    • Movie S1a: top view of Pd13/α-Al2O3(0001) (MP4)

    • Movie S1b: side view of Pd13/α-Al2O3(0001) (MP4)

    • Movie S2a: top view of Pd13/γ-Al2O3(100) (MP4)

    • Movie S2b: side view of Pd13/γ-Al2O3(100) (MP4)

    • Movie S3a: top view of Pd13/γ-Al2O3(110) (MP4)

    • Movie S3b: side view of Pd13/γ-Al2O3(110) (MP4)

    • Movie S4a: top view of α-Al2O3(0001) (MP4)

    • Movie S4b: side view of α-Al2O3(0001) (MP4)

    • Movie S5a: top view of γ-Al2O3(100) (MP4)

    • Movie S5b: side view of γ-Al2O3(100) (MP4)

    • Movie S6a: top view of γ-Al2O3(110) (MP4)

    • Movie S6b: side view of γ-Al2O3(110) (MP4)

    • Movie S7a: top view of Pd13/α-Al2O3(0001) (MP4)

    • Movie S7b: side view of Pd13/α-Al2O3(0001) (MP4)

    • Movie S8a: top view of Pd13/γ-Al2O3(100) (MP4)

    • Movie S8b: side view of Pd13/γ-Al2O3(100) (MP4)

    • Movie S9a: top view of Pd13/γ-Al2O3(110) (MP4)

    • Movie S9b: side view of Pd13/γ-Al2O3(110) (MP4)


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.