Research Article  Open Access
Molecular Dynamics Simulation of Solidification of PdNi Clusters with Different Nickel Content
Abstract
Molecular dynamics simulation has been performed for investigating the glass transition of PdNi alloy nanoparticles in the solidification process. The results showed that the PdNi nanoparticles with composition far from pure metal should form amorphous structure more easily, which is in accordance with the results of the thermodynamic calculation. There are some regular and distorted fivefold symmetry in the amorphous PdNi alloy nanoparticles. The nanoclusters with bigger difference value between formation enthalpies of solutions and glasses will transform to glass more easily than the other PdNi alloy nanoclusters.
1. Introduction
In the past decades, because of the high special surface area and its excellent catalytic performance, pure metal or its alloy nanoparticles have been studied and used widely in catalytic chemistry [1–4]. It is well known that the structure of nanocatalyst particle affects the catalytic activity and selectivity evidently. However, in most experimental cases some restricted conditions were required normally, which blocks the studies on nanostructure of catalyst profoundly. As an economical research method, the computer simulation of molecular dynamics (MD) has been applied to study the metal cluster structure or properties by many researchers successfully [5–11].
Considering the structure stability of catalyst nanoparticles, many MD simulations have been done for researching the melting process of the metal clusters [8–11]. Nevertheless, there are few studies about structure transformation in solidification of the melted clusters. The atomic configuration, which affects the catalyst efficiency, should be different with altering the atomic ratio in the clusters or the solidification process.
PdNi alloy clusters, synthesized by many methods, have been used as catalyst in many reactions, especially in hydrogenation [1]. The melting and solidifying process of PdNi nanoparticles will occur in some preparation or application, such as supporting the cluster on the carbon nanotubes or other catalytic reactions. Therefore, the solidification of melted PdNi nanoparticles has been simulated by molecular dynamics in the present work, and the tendency to forming noncrystalline structure of PdNi alloy with different element content has been also studied here by thermodynamic calculation method.
2. Method of Simulation
The parallel code LAMMPS was used for all MD simulations in this work [12]. The embeddedatom method (EAM), a set of body potentials, has been proved to accurately describe various dynamic properties of transition and noble metals. The EAM potential data used here was obtained from the literature [13]. The Newtonian equations of motion are integrated using the Verlet method with a time step of 5 × 10^{−4} ps.
The calculation is performed in two steps. Firstly, the melting models for simulating solidification process were obtained after the initial models were heated to 2300 K and relaxed enough time (50 ps). Then, the melting models were subjected to cooling process consisting of a series of dropping and keeping temperature MD simulations with temperature decrements of 100 K. The cooling rate of 4 K/ps and the equilibration simulation time of 10 ps were used in the present study. The cooling process of nanoclusters was simulated with constant temperature and constant volume (NVT) molecular dynamic in enough large box with a periodic boundary condition.
The initial models with 1372 atoms (diameter about 2 nm) started with geometries constructed from a FCC block of Pd. Through replacing the Pd atom with Ni atom randomly, the alloy models with different nickel atom content can be created. The exponent of all models can be found in Table 1.

The radial distribution function (RDF) shown as in (1), being regarded as one of the most important parameters, is used to describe the structure characterization of solid, amorphous and liquid states [14]. Consider where is the probability of finding an atom in a distance ranging from to ( is the step of calculation), is the simulated volume of unit cell, is the number of atoms in the system, and is the average number of atoms around the th atom in the sphere shell ranging from to .
Besides MD simulation, Miedema’s thermodynamic theory was used for calculating the forming energies of binary transition metal alloy system, which can explain the glassforming ability of PdNi alloy with different Ni content. On the basis of Miedema’s thermodynamic theory [15–19], the formation enthalpy of PdNi solid solution is expressed as where results from the electronic redistribution occurring when the alloy forms, is the sizemismatch contribution to the formation enthalpy in a binary system, and accounts for the difference in valence and crystal structure of the component metals. However, the structure of palladium is close to that of nickel. That is to say, the can be ignored here.
Then the Gibbs free energy of PdNi solid can be expressed as where is the ideal mixing entropy. Consider where , , , and are the mole ratio and atomic volume of Pd and Ni, respectively. The factor takes the values 8, 5, and 0 for intermetallics, metallic glasses, and solid solutions, respectively; is given by
The parameters, such as , , and , were shown in Table 2. The can also be calculated by this kind of formula.

Similar to the form of , the can be expressed as where and are bulk moduli of Ni and Pd, respectively, and are shear moduli of Ni and Pd, respectively.
The Gibbs free energy of formatting amorphous PdNi alloy is
The contains the free energy of the mixing of the pure liquid metals (the first two terms) plus the free energy difference between the amorphous state and the crystalline state of the pure components. The difference can be expressed as where the is the molar heat of fusion, is the melting points of metal.
For PdNi cluster, surface energy should be considered in the forming energy. For the same size of the clusters and simplifying calculation, the surface energies of noncrystalline and crystalline clusters were considered as equal approximately. Then the term of surface energy was neglected in this study.
3. Results and Discussion
When liquid nanocluster translates to crystal cluster, the jump should appear in caloric curve. However, the cluster only has 1372 atoms, which leads to less energy decrease as crystallization. Therefore, the visible jump is not found in the Figure 1.
Figure 2 shows the radial distribution function (RDF) curves of PdNi alloy clusters at 300 K. Crystalline states are visible through the welldefined peaks, which represent positions of the first, second, third, and th nearest neighbor atoms. It proves the presence of an FCC crystal structure in Pd, Ni200, Ni1200, and Ni models. Because of the equilibrium distance of Pd atom bigger than that of Ni, the Pd atom content increasing, the first neighbor distance generally increasing, which is found in Figure 2. For Ni400, Ni600, Ni800, and Ni1000 clusters, the number of the first neighbor atoms is fewer than that of Pd, Ni200, Ni1200, and Ni, while the splitting of the second peak, denoted in Figure 2 by single arrow line, appears in the RDF curves of Ni400, Ni600, Ni800, and Ni1000 clusters. The splitting of the second peak is a wellknown characteristic feature of solid amorphous structure of metallic glass, and the number of the first neighbor atoms is fewer, which can effectively describe the characteristics of the geometric structural evolvement [20–22].
Pair analysis (PA) technique, which can effectively describe the characteristics of the geometric structural evolvement, was used for the analysis of the geometric features of the atomic cluster [20–22]. In this study, PA method was used to analyze the structural changes accompanying the solidification process of melting PdNi clusters. Based on the regulation of bond pair, two atoms are within a specified cutoff distance of each other and they are called a bonded pair of clusters [22, 23]. In RDF curves, if is the first minimum value, then was defined as cutoff distance for PA. Fourindex number is used to express bonded pairs of atomic clusters by Honeycutt and Andersen [20]. If any atomic pair AB forms a bond, and otherwise ; refers to the number of near neighbors which form bonds with both atom A and atom B; stands for the number of pairs among the neighboring atoms forming bonds; is a special distinguished index parameter. Based on the PA technique, the 1201 and 1311 bond pairs represent the rhombus symmetrical features of shortrange order. The FCC structure has the type of 1421 bond pairs, whereas the HCP crystal has the equal number of 1421 and 1422 bond pairs. The difference between 1421 and 1422 bond pairs is the topological arrangement of the two bonds between the four neighbors. The bond pair 1551, corresponding to a pentagonal bipyramid, is the characteristic of icosahedral order.
Figure 3 shows the normalized abundance of selected pairs and polyhedra in the PdNi clusters simulation in the present study. The results indicate that the number of the 1201, 1311, 1321, 1421, 1422, 1551, 1541, and 1431 pairs will change with temperature dropping, especially when the temperature lower than 1000 K. Because all the clusters are in the liquid state at high temperature, the 1201 and 1311 pairs with relative high proportion shown in Figures 3(a) and 3(b), indicating shortrange order structure, exist in all the clusters. According to the case of proportion of every type pair, all the clusters calculated in present work can be assorted into two kinds. The change tendency of the Ni, Ni1200, Ni200, and Pd, labeled with the first type, is different with the second type clusters including Ni400, Ni600, Ni800, and Ni1000. With the temperature dropping, solidification process leads to the proportion of 1201 pairs in all PdNi clusters decreasing gradually. For the first type clusters, the proportion of 1421 and 1422 pairs increased observably to the relative high value shown in Figures 3(c) and 3(d), which indicates that most atoms in Ni, Ni1200, Ni200, and Pd are arranged mainly as FCC or HCP structure. And this result is in accordance with characteristic peaks shown in Figure 3. There are also 15–20% 1311 pairs in the first type clusters due to the surface atoms packing for the small surface energy. It can estimate the crystalline point where the 1421 pair number is evaluated abruptly. And the crystalline temperature of Ni, Ni1200, Ni200, and Pd nanoclusters was about 800 K, 700 K, 600 K, and 800 K, respectively.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
Relative amount of various bonded pairs in the second type PdNi clusters was different from that of the first kind of clusters. In these clusters, the notable increase of 1311, 1321, 1551, 1541, and 1431 pairs was presented as the temperature dropping took place. These results suggested that the main atom structure of the second PdNi clusters was shortrange order at low temperature; namely, Ni400, Ni600, Ni800, and Ni1000 clusters are in amorphous states. In addition, the three types of pairs (1551, 1541, and 1431) corresponding to inherent structures with regular and distorted fivefold symmetries account for about 25% of all pairs, which implies that part of shortrange order is local icosahedral order, meaning the order characteristic of a 13atom icosahedron. There is also small proportion of 1421 and 1422 bond pairs in the second type clusters which indicates that crystalline shortrange order exists in the disorder system [20, 21].
The results of MD simulation, which indicate that the clusters were easier to form amorphous structure when the number of Pd was close to that of Ni atom, can also be explained by thermodynamic calculation. Figure 4(a) shows the formation enthalpies of glasses () and solutions () of PdNi alloy clusters calculated by the Miedema’s model [16]. It is commonly regarded that the alloys could be vitrified into glass more easily only when the enthalpies of glasses are bigger than those of solution. So the compositions of a possible amorphous structure former should be ranging from about Pd75Ni25 to Pd28Ni72, being marked by dash line in Figure 4(a). In this composition region, the difference value between formation enthalpies of solutions and glasses has been used for expressing the glass formation ability (GFA) of PdNi alloy clusters.
(a)
(b)
Evidently, in the curve of versus nickel atom content, shown in Figure 4(b), the cluster with about 50 percent of nickel can form glass more easily than other clusters. More leaving away from this composition point the value of the difference is much smaller, which suggested that the GFA increases firstly and then decreases with more nickel atoms in the cluster. For the verification of this law, the slower cooling process, in which the cooling rate was 0.8 K/ps and the equilibration simulation time was 100 ps, has also been used to simulate solidification of PdNi alloy clusters. Figure 5 shows that the Pd, Ni200, Ni1200, and Ni clusters were presenting more perfect crystal structure and Ni1000 cluster showed the crystalline structure after the solidification process with slower cooling rate.
4. Conclusion
In summary, the solidification of PdNi clusters with different nickel content has been studied in the present work. The results of molecule dynamic simulation and thermodynamic study certified that PdNi clusters with composition close to 50 at.% will form noncrystal structure more easier than other PdNi clusters.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This study was financially supported by the National Natural Science Foundation of China (nos. 51207031 and 51177022), China Postdoctoral Science Foundation (Grant no. 2013M541368), the Promotive Research Fund for Excellent Young and MiddleAged Scientists of Shandong Province (no. BS2011NJ002), the Science Foundation of HIT at Weihai (Grant no. IMJQ10080026), and the Natural Scientific Research Innovation Foundation of HIT (Grant no. IMXK57080026).
References
 Z. J. Wu, M. H. Zhang, Z. F. Zhao, W. Li, and K. Tao, “Synthesis of a Pd on NiB nanoparticle catalyst by the replacement reaction method for hydrodechlorination,” Journal of Catalysis, vol. 256, no. 2, pp. 323–330, 2008. View at: Publisher Site  Google Scholar
 C.H. Lu and F.C. Chang, “Polyhedral oligomeric silsesquioxaneencapsulating amorphous palladium nanoclusters as catalysts for heck reactions,” ACS Catalysis, vol. 1, no. 5, pp. 481–488, 2011. View at: Publisher Site  Google Scholar
 G. Z. Hua, T. Sharifi, F. Nitze, H. R. Barzegar, C.W. Tai, and T. Wagberg, “Phasetransfer synthesis of amorphous palladium nanoparticlefunctionalized 3D helical carbon nanofibers and its highly catalytic performance towards hydrazine oxidation,” Chemical Physics Letters, vol. 543, pp. 96–100, 2012. View at: Publisher Site  Google Scholar
 S. Y. Shen, T. S. Zhao, J. B. Xu, and Y. S. Li, “Synthesis of PdNi catalysts for the oxidation of ethanol in alkaline direct ethanol fuel cells,” Journal of Power Sources, vol. 195, no. 4, pp. 1001–1006, 2010. View at: Publisher Site  Google Scholar
 D. H. Kim, H. Y. Kim, H. G. Kim, J. H. Ryu, and H. M. Lee, “The solidtoliquid transition region of an AgPd bimetallic nanocluster,” Journal of Physics Condensed Matter, vol. 20, no. 3, Article ID 035208, 2008. View at: Publisher Site  Google Scholar
 L. Miao, V. R. Bhethanabotla, and B. Joseph, “Melting of Pd clusters and nanowires: a comparison study using molecular dynamics simulation,” Physical Review B, vol. 72, no. 13, Article ID 134109, 12 pages, 2005. View at: Publisher Site  Google Scholar
 L. Zhu, K. S. Liang, B. Zhang, J. S. Bradley, and A. E. DePristo, “Bimetallic PdCu catalysts: Xray diffraction and theoretical modeling studies,” Journal of Catalysis, vol. 167, no. 2, pp. 412–416, 1997. View at: Publisher Site  Google Scholar
 S. Valkealahti and M. Manninen, “Melting of copper clusters,” Computational Materials Science, vol. 1, no. 2, pp. 123–134, 1993. View at: Publisher Site  Google Scholar
 Y. Qi, T. Çain, W. L. Johnson, and W. A. Goddard III, “Melting and crystallization in Ni nanoclusters: the mesoscale regime,” Journal of Chemical Physics, vol. 115, no. 1, pp. 385–394, 2001. View at: Publisher Site  Google Scholar
 Y. N. Zhang, L. Wang, W. M. Wang, X. F. Liu, X. L. Tian, and P. Zhang, “Pressure effect on the structural transition of liquid Au,” Physics Letters A, vol. 320, no. 56, pp. 452–458, 2004. View at: Publisher Site  Google Scholar
 L. Wang, Y. N. Zhang, X. F. Bian, and Y. Chen, “Melting of Cu nanoclusters by molecular dynamics simulation,” Physics Letters A, vol. 310, no. 23, pp. 197–202, 2003. View at: Publisher Site  Google Scholar
 S. J. Plimpton, “Fast parallel algorithms for shortrange molecular dynamics,” Journal of Computational Physics, vol. 117, no. 1, pp. 1–19, 1995. View at: Publisher Site  Google Scholar
 X. W. Zhou, R. A. Johnson, and H. N. G. Wadley, “Misfitenergyincreasing dislocations in vapordeposited CoFe/NiFe multilayers,” Physical Review B, vol. 69, no. 14, Article ID 144113, 10 pages, 2004. View at: Publisher Site  Google Scholar
 S. Sengul and S. S. Dalgic, “Size dependence of melting process of ZnSe nanowires: molecular dynamics simulations,” Journal of Optoelectronics and Advanced Materials, vol. 13, no. 1112, pp. 1542–1547, 2011. View at: Google Scholar
 D. Kim, B.J. Lee, and N. J. Kim, “Thermodynamic approach for predicting the glass forming ability of amorphous alloys,” Intermetallics, vol. 12, no. 1011, pp. 1103–1107, 2004. View at: Publisher Site  Google Scholar
 F. R. de Boer, R. Boom, W. C. M. Mattens, A. R. Miedema, and A. K. Niessen, Cohesion in Metals: Transition Metals Alloys, NorthHolland, Amsterdam, The Netherlands, 1988.
 J. A. V. Butler, “The thermodynamics of the surfaces of solutions,” Proceedings of the Royal Society A, vol. 135, no. 827, pp. 348–375, 1932. View at: Publisher Site  Google Scholar
 Z. Y. Qiao, L. J. Yan, Z. M. Cao, and G. R. Liu, “Surface tension calculation of undercooled alloys,” Rare Metals, vol. 20, no. 1, pp. 58–61, 2001. View at: Google Scholar
 J. Brillo, Y. Plevachuk, and I. Egry, “Surface tension of liquid AlCuAg ternary alloys,” Journal of Materials Science, vol. 45, no. 19, pp. 5150–5157, 2010. View at: Publisher Site  Google Scholar
 J. D. Honeycutt and H. C. Andersen, “Molecular dynamics study of melting and freezing of small LennardJones clusters,” The Journal of Physical Chemistry, vol. 91, no. 19, pp. 4950–4963, 1987. View at: Publisher Site  Google Scholar
 L. Hui, F. Pederiva, G. H. Wang, and B. L. Wang, “Local clusters and defects in onedimensional gold wires,” The Journal of Chemical Physics, vol. 119, no. 18, pp. 9771–9776, 2003. View at: Publisher Site  Google Scholar
 H. Li, F. Ding, J. L. Wang, X. F. Bian, and G. H. Wang, “Structural studies of clusters in melt of FeAl compound,” The Journal of Chemical Physics, vol. 114, no. 14, pp. 6413–6416, 2001. View at: Publisher Site  Google Scholar
 H. Li, G. Wang, F. Ding, J. L. Wang, and W. F. Shen, “Molecular dynamics computation of clusters in liquid FeAl alloy,” Physics Letters A, vol. 280, no. 56, pp. 325–332, 2001. View at: Publisher Site  Google Scholar
 J. A. Dean, Lange's Handbook of Chemistry, R. R. Donnelley & Sons, New York, NY, USA, 15th edition, 1999.
 A. K. Niessen, F. R. de Boer, R. Boom, P. F. de Châtel, W. C. M. Mattens, and A. R. Miedema, “Model predictions for the enthalpy of formation of transition metal alloys II,” Calphad, vol. 7, no. 1, pp. 51–70, 1983. View at: Publisher Site  Google Scholar
 W. Martienssen and H. Warlimont, Springer Handbook of Condensed Matter and Materials Data, Springer, Berlin, Germany, 2005.
Copyright
Copyright © 2014 Chen Gang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.