English

• 单原子催化剂(SAC)在催化领域是相对较新的概念，在过去几年中受到人们的广泛关注^[1-8]。由于它具有很高的化学势，金属单原子很容易迁移并聚集形成大的纳米颗粒。如果一个金属单原子想在高温时稳定存在，那么就要求单原子和金属基底之间具有强相互作用。较强的相互作用可以阻碍单个原子凝聚烧结，从而保持稳定。如果这个相互作用过于强，催化剂会失去它的催化活性。因此平衡催化剂活性和稳定性这两者的关系是很有必要的。单原子催化剂烧结过程包括以下几步：首先原子从基底表面脱离成为小的纳米颗粒，之后发生迁移从而形成较大的纳米颗粒，也就会发生凝聚烧结现象^[9]。一般而言，SAC催化剂往往是过渡金属负载在过渡金属氧化物上，金属-载体的强相互作用有利于SAC的稳定，如Datye课题组^[6]发现氧化铈与铂的相互作用可以稳定氧化铈部分晶面上孤立的Pt原子，可以在高温时保持分散的Pt原子不聚集。另一方面，对于高度配位不饱和的过渡金属单原子催化剂(如ad-atom)则研究的还不够。

  Kentaro等通过场离子显微镜研究了Pt(111)上Pt团簇的形成以及迁移能垒与温度的关系，发现在350 K的低温下，包括五聚体在内的Pt团簇扩散和解离现象都是比较显著的^[10]。Kokalj及其合作者^[11]发现Rh₁/Rh(111)-ad-tom催化剂可以增强甲烷第一次脱氢形成甲基的反应，同时可以抑制甲基的进一步脱氢反应，从而既保证了甲烷的充分活化，同时也避免了积碳现象的发生。那么，该类ad-atom催化剂在反应温度下是否可以稳定存在呢?为探讨该问题，在本工作中，我们通过反应力场(ReaxFF)结合LAMMPS(large-scale atomic/molecular massively parallel simulator)软件包进行了大尺度分子动力学模拟，研究了一系列ad-atom过渡金属催化剂(Cu，Ag，Au，Fe，Ni，Pd，Pt，Ru)在高温的稳定性。同时，通过与密度泛函(DFT)计算结果进行对比，验证Ni/C/H/O力场参数，并探究了金属单原子催化剂模型Ni₁/Ni(111)在H₂和O₂气氛下的稳定性行为。该研究工作为人们制备高活性和选择性ad-atom催化剂提供了一定的理论指导依据。

  1.   计算方法

  采用基于密度泛函理论的从头量子力学的程序包VASP(Vienna ab-initio simulation package)^[12-13]进行计算。基组采用平面波赝势方法^[14-15]，利用PBE关联梯度修正泛函对模型进行充分优化^[16]，采用广义梯度近似(GGA)描述体系的交换关联能。平面波的截断能设置为400 eV，k点^[17]设置为3×3×1。对于含有Fe和Ni的体系，加入了自旋极化。晶胞厚度为4层，扩胞3×3，单原子催化剂模型(M₁/M(111))为添加一个额外的原子在fcc或hcp位点。在优化过程中，上3层驰豫，下2层原子保持固定。

  大尺度分子动力学模拟软件包LAMMPS^[18]和反应性力场(ReaxFF)用于分子动力学的计算。动力学模拟所用模型为Cu单胞切(111)表面，厚度为4层，扩9×9超胞的周期性slab模型。这4层中下2层固定，上面2层放开，在表层上方有1.5 nm的真空层。时间步长设为0.25 fs，模拟时间设为100 ps，系综选择NVT系综。首先设置初始温度100 K模拟10 ps稳定，然后用50 ps将系统升温至所需温度(最终温度为300，400，500和600 K)，温度阻尼为25.0，最后模拟100 ps来稳定此状态。SAC体系(Ag，Au，Fe，Ni，Pd，Pt，Ru和Cu)反应性力场参数来源于Duin等的工作^[19-28]。Ni/C/H/O的力场参数基于Ni/C/H和Pt/O的力场参数^{[23, 25]}。

  2.   结果与讨论

  2.1   验证反应力场参数的可靠性

  我们首先对力场参数的可信度进行了验证。在训练集中加入了Ni块体(bulk)的状态方程，金属Ni的内聚能(E_cohesive, 见式1，E_M^bulk表示金属M块体的能量，N_M^bulk表示块体中包含有金属M原子的数目，E_M^atom代表气相中金属M原子的能量)，CH₃，CH₂，CH，H₂和O₂等吸附物种在Ni(111)表面上的结合能(E_binding, 见式2，E_sp/Ni(111)和E_Ni(111)分别代表吸附物种吸附后的体系和吸附前Ni(111)slab的能量，E_sp^atom表示吸附物种在气相中的能量)，Ni(111)面的表面能(E_surf，见式3，E_slab，E_M^bulk和N_M分别为Ni(111)slab模型的能量，块体中金属Ni原子的能量和slab模型中Ni原子数目，A是切面的面积)，以及H₂和O₂在载体上分解的反应热。采用单参数搜索法^[29]，优化了Ni和H，Ni和O的键参数，以及Ni和H，O的交叉项。

  $ {E_{{\rm{cohesive}}}} = {E_{\rm{M}}}^{{\rm{bulk}}} - {N_{\rm{M}}}^{{\rm{bulk}}} \times {E_{\rm{M}}}^{{\rm{atom}}} $

  (1)

  $ {E_{{\rm{ads}}}} = {E_{x{\rm{/N}}i\left( {111} \right)}} - {E_{{\rm{Ni(111)}}}} - {E_x}^{{\rm{atom}}} $

  (2)

  $ {E_{{\rm{surf}}}} = \frac{{{E_{{\rm{slab}}}} - {N_{\rm{M}}} \times E_{\rm{M}}^{{\rm{bulk}}}}}{{2A}} $

  (3)

  图 1和表 1结果表明力场参数对所研究体系吸附能和反应热进行了拟合，说明本文所用力场参数是可靠的。

  图 1

  

  图 1.  不同吸附物种吸附在Ni(111)不同位点以及H₂和O₂吸附在Ni₁/Ni(111) top位点上DFT和反应力场拟合的吸附能结果

  Figure 1.  ReaxFF fits for the binding energy of different adsorbed species (CH₃, CH₂, CH, O and H) on the different site of Ni(111) surface and H₂, O₂ molecular on the Ni₁/Ni(111)

  下载: 全尺寸图片 幻灯片

  表 1

  

  表 1  DFT计算结果和ReaxFF拟合的H₂和O₂分别在Ni(111)和Ni₁/Ni(111)上分解的反应热的比较

  Table 1.  ReaxFF fits for the reaction energy of the dissociation of H₂ and O₂ on Ni(111) and Ni/Ni(111) respectively

  下载: 导出CSV

  kJ·mol-1
  	H2 on Ni(111)	H2 on Ni1/Ni(111)	O2 on Ni(111)	O2 on Ni1/Ni(111)
  ReaxFF	-115.31	-29.92	-169.08	-258.44
  QM	-83.89	-32.38	-162.84	-313.80

  2.2   单原子金属催化剂(M₁/M(111), M=Cu, Ag, Au, Pd, Pt, Ni, Ru, Fe和Fe₁/Fe(100))的稳定性

  对于单原子金属催化剂模型，我们使用分子动力学方法研究了具有分散的单原子表面在不同温度下的热稳定性。通过DFT对四层扩胞3×3的Cu₁/Cu(111)进行结构优化，然后将优化后的构型再次扩胞3×3最终形成一个四层扩胞9×9的结构，在nvt系综100 K的低温下模拟25 ps达到平衡。模拟温度从100 K逐步升温至300 K，并在300 K温度下模拟100 ps，对模拟三次后的结果进行统计。发现单金属原子Cu催化剂会发生聚集(图 2)，并且随着温度增加聚集的原子数增加。例如在300 K时有4~5个原子聚集形成五聚体，而400~500 K时9个Cu原子发生聚集，甚至在600 K时表层原子排列变得更加混乱，Cu原子显著驰豫。这表明单原子Cu在高于300 K时不能稳定存在。

  图 2

  

  图 2.  从上至下分别为Cu, Ag, Au, Ni, Pd, Pt, Ru, Fe单金属M₁/M(111)模型和Fe₁/Fe(100)模型在300, 400, 500和600 K时动力学构型(100 ps)的侧视图

  Figure 2.  Side view of dynamics configurations (100 ps) Cu, Ag, Au, Ni, Pd, Pt, Ru, Fe single metal model from top to bottom respectively and Fe₁/Fe(100) model at 300, 400, 500 and 600 K

  下载: 全尺寸图片 幻灯片

  对于大部分过渡金属(如Ag，Au，Pd，Pt和Ni)，随着温度升高，单金属原子模型变得不稳定并聚集在一起形成金属团簇。有趣的是对于Fe₁/Fe(100)模型，它可以稳定地存在不产生团簇(图 2)，但Fe₁/Fe(111)模型的结构不能保持稳定。这说明金属表面结构和其原子的排列也是影响单原子金属催化剂稳定性的重要因素。如图 3所示，通过DFT计算可知大多数过渡金属(除了Fe)均具有内聚能(E_cohesive)高于结合能(E_binding)的性质(计算方法见公式(1)和公式(2))，说明单原子的聚集是热力学有利的过程，因此倾向于聚集。然而对于铁，它的内聚能小于结合能，单金属铁原子可以被有效地束缚住，意味着聚集过程不是热力学有利过程，因此Fe单原子可以稳定地存在。

  图 3

  

  图 3.  金属单原子和载体的吸附强度以及对应的金属原子的聚合能两者之间的关系

  Figure 3.  Relationship between binding strength of single atom and the corresponding cohesive energy

  下载: 全尺寸图片 幻灯片

  高温时大部分ad-atom催化剂不稳定，因此研究低温时M₁/M(111)表面单原子团聚成簇现象与温度的关系(如图 4所示)，同样温度下团簇中原子数越少越稳定。可以看出低于200 K稳定性Au(Ru)>Pt>Cu>Pd>Ag>Ni，300 K时Au>Ru(Cu)>Pt(Pd)>Ag(Ni)。而Pt₁/Pt(111)在实验中160 K时形成二聚体，210 K时形成三聚体^[10]，与本工作中模拟结果(165 K时形成二聚体，200 K形成三聚体)比较接近。

  图 4

  

  图 4.  M₁/M(111)模型低温时表面单原子发生团聚现象, 温度与团簇原子数的关系

  Figure 4.  Relationship between temperature and atom number of cluster when the aggregation of M₁/M(111) model accurs at low temperature

  下载: 全尺寸图片 幻灯片

  2.3   Ni₁/Ni(111)单原子金属模型催化剂在H₂(O₂)气氛下的稳定性研究

  在实验上，催化剂一般是在各种气氛下参与反应，如Ni可以高效地使甲烷部分氧化。因此我们也研究了H₂和O₂气氛是如何影响模型催化剂(Ni₁/Ni(111))的稳定性。100 K时(图 5)，Ni₁/Ni(111)催化剂在H₂气氛和O₂气氛下保持稳定，在真空环境中Ni单原子则会聚集。在H₂气氛下，对Ni₁/Ni(111)进行动力学研究。发现100 K的模拟温度下，时间进行至9 ps时出现第1个H₂分子分解(图S1)，至24 ps时出现第2个H₂分子分解，之后则无明显变化，并且Ni₁/Ni(111)模型催化剂保持稳定。300 K条件时，有更多的H₂分解，同时催化剂无法保持稳定，表面孤立的Ni单原子团聚成簇。在O₂气氛下，负载的Ni单原子移动受到O₂分子阻碍，在100 K时保持稳定。当模拟温度设定为300 K时，升温过程中出现O₂分子向邻近Ni原子移动的情况(图S1d)，持续升温会出现O₂分解现象，随后Ni单原子发生聚集。体系升至300 K，全部Ni单原子团聚成簇。而随着温度继续升高至400 K以上，表面Ni原子会脱出表面，分解的少部分O原子则会嵌入表面空穴中。与真空状态模拟结果相比较(图S2)，H₂(尤其是O₂)的存在，一定程度上提高了Ni₁/Ni(111)的热稳定性。其原因可能是Ni与O₂(或O)之间强的相互作用阻碍了聚合现象的发生。

  图 5

  

  图 5.  在H₂ (a)和O₂气氛(b)下，Ni₁/Ni(111)催化模型在100、300、400、500和600 K时的动力学模拟结果

  Figure 5.  Dynamics results of Ni₁/Ni(111) catalyst model under H₂ (a) and O₂ (b) atomosphere at 100, 300, 400, 500 and 600 K

  下载: 全尺寸图片 幻灯片

  3.   结论

  在本工作中，通过反应性力场分子动力学模拟研究了金属载体的单原子催化剂(M₁/M(111))在300到600 K的热稳定性。模拟结果表明，对于SAC模型，只有Fe₁/Fe(100)催化剂体系可以在高的反应温度(约600 K)中稳定存在，而别的单原子金属催化剂会在载体表面形成金属团簇。原因总结如下：Fe₁/Fe(100)的吸附有较强的结合能，高于它的内聚能，而别的过渡金属相反，与载体的结合能力要弱于内聚能，因此Fe₁/Fe(100)可以稳定存在。对于Ni₁/Ni(111)催化剂，H₂和O₂气氛在一定程度上提高了它的稳定性，气氛会阻碍Ni单原子在表面发生团聚的过程。

  Supporting information is available at http://www.wjhxxb.cn

  
  
• 1. [1]

     Yang X F, Wang A Q, Qiao B T, et al. Acc. Chem. Res., 2013, 46(8):1740-1748 doi: 10.1021/ar300361m
     

  2. [2]

     Kyriakou G, Boucher M B, Jewell, et al. Science, 2012, 335(6073):1209-1212 doi: 10.1126/science.1215864
     

  3. [3]

     Lucci F R, Liu J L, Marcinkowski M D, et al. Nature Com-munications., 2015, 6:8550 doi: 10.1038/ncomms9550
     

  4. [4]

     Shan J J, Li M W, Allard L F, et al. Nature, 2017, 551:605-608 doi: 10.1038/nature24640
     

  5. [5]

     Kwon Y W, Kim T Y, Kwon G, et al. J. Am. Chem. Soc., 2017, 139(48):17694-17699 doi: 10.1021/jacs.7b11010
     

  6. [6]

     Jones J, Xiong H F, DeLaRiva A T, et al. Science, 2016, 353:150-154 doi: 10.1126/science.aaf8800
     

  7. [7]

     Nie L, Mei D H, Xiong H F, et al. Science, 2017, 358(6369):1419-1423 doi: 10.1126/science.aao2109
     

  8. [8]

     Wang A Q, Li J, Zhang T. Nat. Rev. Chem., 2018, 2:65-81 doi: 10.1038/s41570-018-0010-1
     

  9. [9]

     Liu J C, Wang Y G, Li J. J. Am. Chem. Soc., 2017, 139(17):6190-6199 doi: 10.1021/jacs.7b01602
     

  10. [10]

      Kyuno K, Gert E. Surf. Sci., 1999, 437(1):29-37
      

  11. [11]

      Kokalj A, Bonini N, Sbraccia C, et al. J. Am. Chem. Soc., 2004, 126(51):16732-16733 doi: 10.1021/ja045169h
      

  12. [12]

      Kresse G, Furthmüller J. Comput. Mater. Sci., 1996, 6(1):15-50
      

  13. [13]

      Kresse G, Furthmüller J. Phys. Rev. B, 1996, 54(16):11169-11186 doi: 10.1103/PhysRevB.54.11169
      

  14. [14]

      Blöchl P E. Phys. Rev. B, 1994, 50(24):17953-17979 doi: 10.1103/PhysRevB.50.17953
      

  15. [15]

      Kresse G J, Daniel J. Phys. Rev. B, 1999, 59(3):1758-1775 doi: 10.1103/PhysRevB.59.1758
      

  16. [16]

      Perdew J P, Burke K, Ernzerhof M. Phys. Rev. Lett., 1996, 77(18):3865-3868 doi: 10.1103/PhysRevLett.77.3865
      

  17. [17]

      Monkhorst H J, Pack J D. Phys. Rev. B, 1976, 13(12):5188-5192 doi: 10.1103/PhysRevB.13.5188
      

  18. [18]

      Plimpton S. J. Comput. Phys., 1995, 117(1):1-19
      

  19. [19]

      Van Duin A C T, Dasgupta S, Lorant F, et al. J. Phys. Chem. A, 2001, 105(41):9396-9409 doi: 10.1021/jp004368u
      

  20. [20]

      Adam L, David C, Van Duin A C T, et al. Surf. Sci., 2016, 645:67-73 doi: 10.1016/j.susc.2015.11.009
      

  21. [21]

      Bae G T, Aikens C M. J. Phys. Chem. A, 2013, 117(40):10438-10446 doi: 10.1021/jp405992m
      

  22. [22]

      Masoud A, Van Duin A C T, James D K. J. Phys. Chem. A, 2010, 114(21):6298-6307 doi: 10.1021/jp101332k
      

  23. [23]

      Jonathan E M, Van Duin A C T, William A G. J. Phys. Chem. C, 2010, 114(11):4939-4949 doi: 10.1021/jp9035056
      

  24. [24]

      Thomas P S, Michael J J, Van Duin A C T. J. Phys. Chem. C, 2014, 118(9):4967-4981 doi: 10.1021/jp411015a
      

  25. [25]

      Fantauzzi D, Bandlow J, Van Duin A C T, et al. Phys. Chem. Chem. Phys., 2014, 16:23118-23133 doi: 10.1039/C4CP03111C
      

  26. [26]

      Kim S Y, Lee H W, Pai S J, et al. ACS Appl. Mater. Interfaces, 2018, 10(31):26188-26194 doi: 10.1021/acsami.8b05070
      

  27. [27]

      Yeon J, Adams H L, Junkermeier C E, et al. J. Phys. Chem. B, 2018, 122(2):888-896
      

  28. [28]

      Shin Y K, Gai L, Raman S, et al. J. Phys. Chem. A, 2016, 120(41):8044-8055 doi: 10.1021/acs.jpca.6b06770
      

  29. [29]

      Van Duin A C T, Baas J M A, van de Graaf B. J. Chem. Soc. Faraday Trans., 1994, 90(19):2881-2895 doi: 10.1039/ft9949002881
      

• 

  图 1  不同吸附物种吸附在Ni(111)不同位点以及H₂和O₂吸附在Ni₁/Ni(111) top位点上DFT和反应力场拟合的吸附能结果

  Figure 1  ReaxFF fits for the binding energy of different adsorbed species (CH₃, CH₂, CH, O and H) on the different site of Ni(111) surface and H₂, O₂ molecular on the Ni₁/Ni(111)

  下载: 全尺寸图片 幻灯片
  

  图 2  从上至下分别为Cu, Ag, Au, Ni, Pd, Pt, Ru, Fe单金属M₁/M(111)模型和Fe₁/Fe(100)模型在300, 400, 500和600 K时动力学构型(100 ps)的侧视图

  Figure 2  Side view of dynamics configurations (100 ps) Cu, Ag, Au, Ni, Pd, Pt, Ru, Fe single metal model from top to bottom respectively and Fe₁/Fe(100) model at 300, 400, 500 and 600 K

  下载: 全尺寸图片 幻灯片
  

  图 3  金属单原子和载体的吸附强度以及对应的金属原子的聚合能两者之间的关系

  Figure 3  Relationship between binding strength of single atom and the corresponding cohesive energy

  下载: 全尺寸图片 幻灯片
  

  图 4  M₁/M(111)模型低温时表面单原子发生团聚现象, 温度与团簇原子数的关系

  Figure 4  Relationship between temperature and atom number of cluster when the aggregation of M₁/M(111) model accurs at low temperature

  下载: 全尺寸图片 幻灯片
  

  图 5  在H₂ (a)和O₂气氛(b)下，Ni₁/Ni(111)催化模型在100、300、400、500和600 K时的动力学模拟结果

  Figure 5  Dynamics results of Ni₁/Ni(111) catalyst model under H₂ (a) and O₂ (b) atomosphere at 100, 300, 400, 500 and 600 K

  下载: 全尺寸图片 幻灯片

  表 1  DFT计算结果和ReaxFF拟合的H₂和O₂分别在Ni(111)和Ni₁/Ni(111)上分解的反应热的比较

  Table 1.  ReaxFF fits for the reaction energy of the dissociation of H₂ and O₂ on Ni(111) and Ni/Ni(111) respectively

  kJ·mol-1
  	H2 on Ni(111)	H2 on Ni1/Ni(111)	O2 on Ni(111)	O2 on Ni1/Ni(111)
  ReaxFF	-115.31	-29.92	-169.08	-258.44
  QM	-83.89	-32.38	-162.84	-313.80

  下载: 导出CSV
•