Analysis of Pd-Ni Nanobelts Melting Process Using Molecular Dynamics Simulation

The melting process of Pd-Ni alloy nanobelts with different Ni atom content has been simulated by molecular dynamic (MD) method.The radial distribution function, the Lindemann index, and pair analysis methodwere used to characterize Pd-Ni nanobelt models in simulation. The results indicate that the melting temperature of Pd-Ni nanobelt with composition far from pure metal was lower than that of othermodels, and the breaking point of the nanobelt can be illustrated by the Lindemann index. Pair analysis indicates that the number of FCCpairswill decrease and almost disappear atmelting pointwith increasing temperature.Themelting points of Pd-Ni alloy nanobelts were also calculated by thermodynamic method, and the results were close to that obtained by MD simulation.


Introduction
After carbon nanotube has been reported in 1991 [1], onedimensional (1D) materials have received much attention due to their excellent physical property and have been expected to be used in a wide range of application fields [2][3][4][5][6].Because of their dramatically different behavior from bulk materials, much research work has been done for the structure and properties of nanowires [7][8][9][10].Thermal stability and melting behavior of nanomaterial are another subject that has been studied theoretically and experimentally for a rather long time.The melting process of nanomaterial determined by experiment is very difficult.Nevertheless, computer simulations provide an excellent method to study these melting processes at an atomic level.Up to date, the melting process and thermal properties of many kinds of 1D metal nanomaterial, including Au [8], ZnSe [11], GaN [12], Zr [13], Ag [14], Cu [15], Ni [16,17], Pd [18], and Cu 3 Au [19] nanowires, have been investigated by molecular dynamic (MD) or other computer calculation methods.However, few efforts have been focused on the thermal stability and melting behavior of the 1D binary alloy nanomaterial with different component.
Pd-Ni alloy can be used as catalyst in many reactions [20,21] and hydrogen sensors [22].However, its catalytic properties and performance of sensor relate to Pd or Ni concentration in the Pd-Ni alloy.Therefore, the structural stability and melting characteristics of Pd-Ni 1D alloy with different nickel concentration are studied by MD in the present work.And the relationship of melting points and fusion enthalpy of 1D Pd-Ni alloy models for MD process has also been investigated by using the model introduced in the literature [23].

Method of Simulation
The parallel code lrge-scale atomic/molecular massively parallel simulator (LAMMPS) [24] was used for all MD simulations.And visualization of MD simulation was performed using visual molecular dynamics (VMD) [25] and open visualization tool (Ovito) [26].The embedded-atom method (EAM), a set of -body potentials have been proved that it can accurately describe various dynamic properties of transition and noble metals.The interactions of atoms were modeled using the EAM potential, which is a convenient model to study metal clusters with FCC structure because the local density effects are included in its parameterization.In this work, the EAM potential data was obtained from the literature [27].
The geometry of nanowire is constructed as a regular FCC lattice with initial surface orientations of [100], [010], and [001] in the , , and  directions of simulation box, respectively.An infinitely long nickel nanowire is modeled by applying periodic boundary conditions in the [001] direction.1.And the integer random distribution number (RND  ) was generated by In order to remove the interactions which might lead to local structural distortion and result in the unstable simulation, energy minimization of the structure has been performed before starting molecular dynamics simulation.The MD simulations were carried out in the canonical ensemble with a constant number of atoms  and volume ().The equations of motion were integrated using the Verlet leapfrog algorithm with a time step of 5 × 10 −4 ps.The nanobelt was heated to 1800 K in increments of 50 K, and the heating rate was 0.5 K/ps.Near the melting point (in the region of 900-1300 K), the temperature increments were reduced to 20 K to account for the large temperature fluctuations.At each temperature 10 5 MD trajectory steps were propagated, which is sufficient to ensure the validity of the results presented here.
The radial distribution function (RDF) shown in (2), being regarded as one of the most important parameters, is used to describe the structure characterization of solid, amorphous, and liquid states [11].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 averaged number of atoms around the th atom in the sphere shell ranging from  to  + Δ.
The Lindemann index was recently applied to characterize nanoparticles, and it was proposed that melting of a nanoparticle at a critical value in the range of 0.03-0.05[28].The Lindemann indices of each atom in the models were calculated at each temperature as where   and  are the Lindemann indices of the th atom and the whole model, respectively; ⟨⋅ ⋅ ⋅ ⟩ T in (3) denotes the thermal average at temperature ;   is the distance between the th and th atoms.Pair analysis (PA) technique, which can effectively describe the characteristics of the geometric structural evolvement, has been used for analysing the geometric features of the atomic cluster [29].In the present study, PA method was used to analyse the structural changes accompanying the melting process of Pd-Ni 1D models.Based on the regulation of bond pair, if two atoms are within a specified cutoff distance, they are called a bonded pair of models [29,30].In RDF curves, (  ) is the first minimum value; then   was defined as cutoff distance for PA.The four-index number  is used to express bonded pairs of atomic clusters [29,30].If any atomic pair A-B formed a bond,  = 1 and otherwise  = 2;  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 short-range 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, 1541, and 1431, corresponding to a pentagonal bipyramid, is the characteristic of icosahedral order.

Results and Discussion
Figure 2 shows the total energy of nanobelts with different Ni content versus temperature (TE curve) during heating process.Total energies increase linearly with temperature in the early stage.When close to the melting transition temperature, simple jumps in total energy can be easily observed.The melting transition temperature of model can be estimated according to the TE curves.However, being different from that of cluster and bulk system, in the TE curve of nanobelt besides sharp increasing region, there is also a decreasing part after melting takes place.Therefore, the derivative curve of TE versus temperature shown in Figure 2, which can be considered as heat capacity [11], should have sharp peak and valley region.The melting points of the nanobelts can be determined by the position of the peak because of the heat absorption of melting.For reducing the total energy, especially surface energy, liquid 1D nanomaterial will change to sphericallike material due to the little interaction of atoms in the model at high temperature.That is to say, the nanobelt will break for decreasing the TE after 1D nanomaterial melting.From the derivative curve, the breaking temperature of 1D Pd-Ni can be estimated by the position of the valley.The melting temperature and breaking temperature of all the nanobelt models can be found in Figure 2.There are two sharp increasing and decreasing regions in the TE curve of PdNi1200 model.The reasons of emerging the first increasing and decreasing region were the melting of the model and decreasing of surface energy.In Figure 3  motion but merely vibrate around their original lattice positions.Following the first stage, the Lindemann index increased and decreased rapidly.It is evident that the melting point and breaking temperature of the nanobelt can be estimated from the derivative curve of the Lindemann index versus temperature.At the temperature corresponding to the peak of derivative curve, lots of atoms in the nanobelt model exhibit large amplitude diffusion, which should be considered as melting of nanobelt.In the curve, there is also a sudden collapse of the 1D structure at the temperature corresponding to the valley.The melting temperature and breaking temperature estimated from derivative curves were the same as those shown in Figure 2. It can be found that with Ni concentration increase in the model, the fusing temperature decreases first and then increases.However, the breaking temperature of nanobelt with different Ni concentration did not obey this rule due to relating with the atom position and diffusion process of models.The Lindemann index of every atom in models can be used to determine the melting point and breaking temperature more precisely.From Figure 5, it can be found that the Lindemann index of atoms in PdNi1600 model was bigger than 0.05 at 1080 K which indicates that melting process occurs at this temperature [28].At 1200 K, though there was no breaking, the Lindemann index of some atoms was very high, which suggests that the nanobelt will break at little higher temperature.Structural properties and changes of nanobelt during heating are of interest in understanding mechanical and other properties of materials.Figure 6 shows that there are only 1421, 1311, and 1211 pairs in the 1D Pd-Ni alloy models at very low temperature.It can be found that the number of 1421 pairs will decrease obviously with increasing temperature.Due to the collapse of 1421 pairs, atoms in the FCC structure should form other pairs, which can be identified by the elevation of the 1201, 1301, and 1311 pair numbers.At melting temperature, the number of 1421 and 1311 pairs in the nanobelt drops abruptly due to melting transition of nanobelt, and the 1421 pairs almost disappear after nanobelt melting.It can also be found that the number of 1201 and 1321 bondtypes increased obviously at melting temperature, which was caused by the collapse of 1311 and 1421 pairs.
Through pure element thermodynamic data, the classical thermodynamics method can be used to predict many thermodynamic parameters of bulk alloy, such as enthalpy and melting point.Therefore, by building the relationship of thermodynamic properties of nanomaterial and bulk material, thermodynamic parameters of nanomaterial can also be calculated by thermodynamics method.In this study, the fusion enthalpy of 1D Pd-Ni alloy model was calculated by ( 5) which was reported by Guisbiers and Buchaillot [31].Consider where  *  and Δ  are the size-dependent melting point and fusion enthalpy, respectively;  ,∞ and Δ ,∞ are the melting point and fusion enthalpy of bulk material, respectively.Here,  *  value can use the results of MD simulation.Thermodynamic parameters of bulk Pd-Ni alloy can be calculated from the pure Ni and Pd thermodynamic data [32] which has been shown in the Appendix.For example, Δ ,∞ , Δ ,∞ , and  ,∞ can be calculated by (6) [33] and (7), respectively.The expression of   and   , which were Gibbs energy of liquid and solid phases of bulk material, used in ( 6) and ( 7) can be illustrated by (8).Then, using the thermodynamic data of bulk material calculated previously, Δ  and Δ  will be obtained by (5).Consider All the calculated results have been shown in Table 2.It can be found that the relation between melting point of nanobelt model and Ni concentration is the same as that of bulk phase Pd-Ni alloy.However, because the specific surface area of nanomaterial is very great, the Pd-Ni nanobelt melt temperature is smaller than that of bulk phase.The melting points of nanobelts can be reproduced by equation of  * *  = Δ  / Δ  .The results shown in Table 2 indicate that the melting point obtained by MD simulation is close to  * *  , and the error is below 5%.Particularly, the difference of melt point between MD and thermodynamic method for PdNi800 and PdNi1200 was only 3-4 K.The nanobelt models used in this work were considered as the ideal solution, so the error may be different when the other mix model was used, such as aggregation state model.

Conclusion
Summarily, the heating process of Pd-Ni alloy nanobelts has been simulated by MD and the melting temperature was also calculated by thermodynamic method.The melting and breaking of Pd-Ni alloy nanobelts can be indicated by increasing and decreasing regions of total energy relating to temperature curve.With increasing temperature, number of FCC pairs will decrease and almost disappear at melting point.The melting point is lower than bulk phase Pd-Ni alloy, and the melt point obtained by MD is close to that of calculated by thermodynamic method.

Figure 2 :
Figure 2: Total energy versus temperature curves for Pd-Ni nanobelt models during heating process.
Figure2shows the total energy of nanobelts with different Ni content versus temperature (TE curve) during heating process.Total energies increase linearly with temperature in the early stage.When close to the melting transition temperature, simple jumps in total energy can be easily observed.The melting transition temperature of model can be estimated according to the TE curves.However, being different from that of cluster and bulk system, in the TE curve of nanobelt besides sharp increasing region, there is also a decreasing part after melting takes place.Therefore, the derivative curve of TE versus temperature shown in Figure2, which can be considered as heat capacity[11], should have sharp peak and valley region.The melting points of the nanobelts can be determined by the position of the peak because of the heat absorption of melting.For reducing the total energy, especially surface energy, liquid 1D nanomaterial will change to sphericallike material due to the little interaction of atoms in the model at high temperature.That is to say, the nanobelt will break for decreasing the TE after 1D nanomaterial melting.From the derivative curve, the breaking temperature of 1D Pd-Ni can be estimated by the position of the valley.The melting temperature and breaking temperature of all the nanobelt models can be found in Figure2.There are two sharp increasing and decreasing regions in the TE curve of PdNi1200 model.The reasons of emerging the first increasing and decreasing region were the melting of the model and decreasing of surface energy.In Figure3(a), part of crystal structure can be found in the model, which results in PdNi1200 model holding 1D nanomaterial.With the temperature increasing, the melting of whole model (Figure3(b)) and the collapse of nanobelt should induce the second increasing and decreasing region in the derivative curve.Figure4shows the temperature dependence upon the Lindemann index of the nanobelts during heating.It is clear from Figure4that the Lindemann index increases slowly and linearly with temperature at low temperature stage because of the linear increase in atomic kinetic energy with temperature.The value of the Lindemann index during this stage is very small since most of the atoms do not have large amplitude

Figure 4 Figure 3 :Figure 4 :
Figure2shows the total energy of nanobelts with different Ni content versus temperature (TE curve) during heating process.Total energies increase linearly with temperature in the early stage.When close to the melting transition temperature, simple jumps in total energy can be easily observed.The melting transition temperature of model can be estimated according to the TE curves.However, being different from that of cluster and bulk system, in the TE curve of nanobelt besides sharp increasing region, there is also a decreasing part after melting takes place.Therefore, the derivative curve of TE versus temperature shown in Figure2, which can be considered as heat capacity[11], should have sharp peak and valley region.The melting points of the nanobelts can be determined by the position of the peak because of the heat absorption of melting.For reducing the total energy, especially surface energy, liquid 1D nanomaterial will change to sphericallike material due to the little interaction of atoms in the model at high temperature.That is to say, the nanobelt will break for decreasing the TE after 1D nanomaterial melting.From the derivative curve, the breaking temperature of 1D Pd-Ni can be estimated by the position of the valley.The melting temperature and breaking temperature of all the nanobelt models can be found in Figure2.There are two sharp increasing and decreasing regions in the TE curve of PdNi1200 model.The reasons of emerging the first increasing and decreasing region were the melting of the model and decreasing of surface energy.In Figure3(a), part of crystal structure can be found in the model, which results in PdNi1200 model holding 1D nanomaterial.With the temperature increasing, the melting of whole model (Figure3(b)) and the collapse of nanobelt should induce the second increasing and decreasing region in the derivative curve.Figure4shows the temperature dependence upon the Lindemann index of the nanobelts during heating.It is clear from Figure4that the Lindemann index increases slowly and linearly with temperature at low temperature stage because of the linear increase in atomic kinetic energy with temperature.The value of the Lindemann index during this stage is very small since most of the atoms do not have large amplitude

Table 1 :
Content of nickel of every 1D Pd-Ni alloy model.