State Key Laboratory of Chemical Engineering, East China University of Science and Technology (ECUST), Shanghai, China.
J Chem Phys. 2011 Apr 7;134(13):134704. doi: 10.1063/1.3575193.
First-principles calculations have been performed to investigate CH(4) dissociation and C diffusion during the Ni∕Fe-catalyzed growth of carbon nanofibers (CNFs). Two bulk models with different Ni to Fe molar ratios (1:1 and 2:1) are constructed, and x-ray diffraction (XRD) simulations are conducted to evaluate their reliability. With the comparison between the calculated and experimental XRD patterns, these models are found to be well suited to reproduce the crystalline structures of Ni∕Fe bulk alloys. The calculations indicate the binding of the C(1) derivatives to the Ni∕Fe closest-packed surfaces is strengthened compared to that on Ni(111), arising from the upshift of the weighted d-band centers of catalyst surfaces. Then, the transition states for the four successive dehydrogenation steps in CH(4) dissociation are located using the dimer method. It is found that the energy barriers for the first three steps are rather close on the alloyed Ni∕Fe and Ni surfaces, while the activation energy for CH dissociation is substantially lowered with the introduction of Fe. The dissolution of the generated C from the surface into the bulk of the Ni∕Fe alloys is thermodynamically favorable, and the diffusion of C through catalyst particles is hindered by the Fe component. With the combination of density functional theory calculations and kinetic analysis, the C concentration in catalyst particles is predicted to increase with the Fe content. Meanwhile, other experimental conditions, such as the composition of carbon-containing gases, feedstock partial pressure, and reaction temperature, are also found to play a key role in determining the C concentration in bulk metal, and hence the microstructures of generated CNFs.
运用第一性原理计算研究了 Ni∕Fe 催化生长碳纳米纤维(CNF)过程中 CH(4)的离解和 C 的扩散。构建了两种 Ni 与 Fe 摩尔比不同的体相模型(1:1 和 2:1),并进行了 X 射线衍射(XRD)模拟以评估其可靠性。通过对计算和实验 XRD 图谱的比较,发现这些模型非常适合再现 Ni∕Fe 体相合金的晶体结构。计算表明,与 Ni(111)相比,C(1)衍生物与 Ni∕Fe 最密排表面的结合得到了增强,这是由于催化剂表面加权 d 带中心上移所致。然后,使用二聚体方法定位了 CH(4)离解的四个连续脱氢步骤的过渡态。结果发现,在合金化的 Ni∕Fe 和 Ni 表面上,前三个步骤的能垒相当接近,而 Fe 的引入显著降低了 CH 离解的活化能。从表面生成的 C 溶解到 Ni∕Fe 合金的体相是热力学有利的,C 通过催化剂颗粒的扩散受到 Fe 组分的阻碍。通过密度泛函理论计算和动力学分析的结合,预测催化剂颗粒中的 C 浓度会随 Fe 含量的增加而增加。同时,其他实验条件,如含碳气体的组成、原料分压和反应温度,也被发现对确定金属体相中 C 的浓度,从而对生成的 CNF 的微观结构起着关键作用。