A Simulation Study for Trimetallic Nanosized Alloy ( Ni , Cu , and Ag ) in Hydrogenation of Organic Compounds : A Case Study of ( Nitrophenols )

Trimetallic system (Ni, Cu, and Ag) supported on alumina was utilized for hydrogenation of nitrophenols. The catalytic active centers for hydrogenation were attributed only to the presence of Ni. However, the presence of bior trimetallic systems improves the catalytic activity via extra synergism.The catalytic activity wasmeasured as the time for reaching 100% conversion.The function of synergism was fitted for both bimetallic systems (Ni:Ag; Ni:Cu) individually. Subsequently, three-dimensional function was fitted for trimetallic system (Ni:Cu:Ag) based on the linear combination of data for individual bimetallic system. After a complex calculation areal function was evaluated. An Excel program was written to simply evaluate the catalytic activity of trimetallic system with high accuracy. Characterization of catalysts was performed using EPR and pulsed chemisorption by hydrogen. These characterizations of samples enable us to evaluate particle size, metallic surface area, and degree of dispersion. These values were successfully correlated with the synergism function. The program written then could be capable of predicting these values for any trimetallic system.


Introduction
Hydrogenation reactions are one of the most important reactions for preparation of many fine chemicals and important intermediates [1][2][3][4][5].The selectivity and activity of such reactions are an issue of most recent researches [1,2].Most of catalysts used in such important reactions are metal catalysts [6,7].The basic idea of hydrogenation reactions is to adsorb both metal and target organic moiety subject to hydrogenation.The challenge in this subject is to find that metal has high activity and selectivity towards the desired compound.In order to attain both activity and selectivity researches tend to use bimetallic system.Thus, a combination of Ni-Ag, Ni-Cu, and Ni-Au was used in [1] for active and selective hydrogenation of cinnamaldehyde.In addition, Ni-Pd bimetallic system was used in hydrogenation of nitrobenzyl ethers to amino benzyl [2].Moreover, Ni-Ag bimetallic system was used in hydrogenation of dimethyl oxalate [3].
Many others researches also used many combinations of bimetallic systems in order to enhance the selectivity and activity [4][5][6][7][8].Searching the literature data, only few papers deal with the trimetallic system.Thus, in [9] homogenous catalyst containing trimetallic system of Sn, Ru, Pt in selective reduction of cyclodecatriene.However heterogeneous catalyst containing clusters of PtRu 5 Sn was evaluated as hydrogenation catalyst [10,11].Just as rare sole example we found a paper that deals with hetero quarter metallic nanoparticle hydrogenation catalyst of Ni/Ru/Pt/Au [12].In all above examples of hydrogenation catalysts containing bi-or tri-or even quaternary metallic systems a simple evaluation of one sample or more is performed just to bring the selectivity and compare it to bi-or monometallic catalyst.In this paper we aim to investigate and write a mathematical simulation program containing trimetallic heterogeneous catalyst containing Ni, Ag, and Cu.
In order to fulfill this purpose, we choose a simple effective hydrogenation reaction.Thus hydrogenation of nitrophenol into amino phenol using supported metal and using hydrazine hydrate as a hydrogen source was proven to be 100% selective and finishes with 100% conversion in only few minutes [13,14].The program will evaluate the effect of this trimetallic system and will build three-dimensional model for the activity function of these catalysts and also will be capable of evaluating some physical properties of these metals.

Preparation of the Catalyst.
Typically 10 g of alumina was impregnated with 5 wt% metals.Different mole ratios of mono-, di-, and trimetal are varied keeping 5 wt% constant.After drying at 100 ∘ C the catalyst is subjected to chemical reduction using hydrazine hydrate in alkaline solution of NaOH.After completion of reduction, the catalyst was rapidly filtered and dried and kept in a closed bottle.

Catalytic Reaction.
Typically, a solution containing 0.125 g of nitrophenols with 20 mL of hydrazine hydrate was heated to 80-100 ∘ C in a three-necked flask connected to a condenser.0.5 g of the catalyst was added to the above solution.The time was recorded just upon the addition.After reaction completion the color of the system is changed to grayish white indicating 100% conversion [14].The filtrate was then evaporated under reduced pressure.The residual solids were recrystallized from hot water to give the pure product of corresponding amino phenols in almost 100% yield.

XRD (X-Ray Diffraction Analysis
). X-Ray diffractograms of various solids were collected using a Bruker D8 advance instrument with CuK1 target with second monochromator 40 kV, 40 mA.

EPR (Electron Paramagnetic Resonance).
The EPR spectra were recorded on EMX Bruker instrument operated at -band frequency.The following parameters are generalized to all samples otherwise mentioned in the text.Microwave frequency is 9.79 GHz.Receiver gain is 20.Sweep width is 6000, center at 3480 Gs.Microwave power is 0.202637 W.

Pulsed Chemisorption.
The H 2 chemisorption experiments were carried out at 373 K using a pulse reactor to determine the Ni metal area and Ni metal particle size.In a typical experiment, approximately 150 mg of the catalyst sample was loaded in a micro quartz reactor (8 mm i.d., and 250 mm long) and the catalyst sample was first subject to hydrogen flow at 420 K for 1 h and kept at the same temperature for 1 h under He flow and then the reactor was cooled to 373 K under helium gas flow.The outlet of the reactor was connected to a microthermal conductivity detector (TCD) of GC-17A (Quantachrome ChemBET 3000, USA) through an automatic six-port valve after cooling of the sample to 303 K, pulses of gas (5% H 2 balance He) (500 mL), until there is no further change in the intensity of TCD peak due to gas pulse.

X-Ray Diffraction (XRD)
. Figure 1 represents the XRD patterns of bare alumina, 5% Ni, and some trimetallic systems.From these figures it could be shown that Ni metal phase was observed only in pure 5% wt Ni.While in the rest of the samples there was no observation of such phase.This could be explained by the lower amount of metal which could not be detected by XRD techniques.Another reason is that the metal phases could be existing in amorphous nanoscale with bi-or trimetallic system.This will be declared further in EPR section.

Electron Spin Resonance (EPR).
EPR spectra of all samples are given in Figure 2. In this figure we can observe an EPR signal characterized for the nanosized nickel [13].

Program for Trimetallic System.
In this research we aimed to study and evaluate trimetallic system of Ni, Cu, and Ag as an effective catalyst for hydrogenation of nitrophenols.In order to study this system we make a series of bimetallic systems Ni:Cu and Ni:Ag with different molar ratios (Table 1).The catalytic activity was taken as time to reach 100% conversion [13].However we added a term "catalytic activity response" as a reciprocal of this time.This was taken as real measure of catalytic activity and thus as time decreases catalytic activity response increases.
From experimental test it was found that neither Cu nor Ag alone gives any catalytic activity in these reactions; however in presence of Ni as a bimetallic system they could enhance the catalytic activity.
In order to study the effect of these metals on the activity of Ni, we first determined what is called the linear catalytic activity which is the catalytic activity of catalyst based on its % nickel theoretically calculated from the pure nickel catalyst.
Based on these values we can check any presence of synergism occurring due to addition of either Cu or Ag as bimetallic system with nickel.A function of synergism was evaluated for bimetallic system individually (Ni:Cu and Ni:Ag), was given below, and was fitted with regression with  2 more than 0.99 as indicated below (Figures 3 and 4).where  is wt% of Cu in the sample and  is % syn/Ni.From Figures 3 and 4 curves it was found that the synergism function is increased as the % added metal (Cu or Ag) increases reaching maximum and then decreases again.
In order to study the trimetallic system we make a linear addition of both functions of individual synergism.Figure 5 represents the 3D graph of raw data given.Figures 5 and 6 showed clearly that the 3D model of the linear synergism functions showed more than 3D peaks varied according to the composition of the system.The maximum synergism was found to be of value about 22% syn/Ni.This 3D figure could give us the ability to evaluate the four parameters in one time (% syn, % Ni, % Cu, and % Ag).However, subsequently, analysis of the real system should be compared which will be performed in the next section.
After applying the previous function on real trimetallic systems we observed that extra synergism occurs due to presence of both Cu and Ag together.In order to evaluate this extra synergism we correlate the difference in synergism with the % Ni (Figure 7).Thus, it was found that extra synergism is observed only in % Ni between 5 and 40%.In order to calculate it automatically we simulate this function as follows: where  is difference between simulated linear synergism and real trimetallic synergism and  is % Ni.After applying the above function on the threedimensional linear synergism function we could obtain the real function of synergism (Figure 8).
Obtaining the previous function enables us to deduce the catalytic activity and consequently the time to reach 100% for trimetallic system simply by input of the mole ratio of Ni:Cu:Ag.
Particle Size Simulation.Particle size was evaluated by Kawabata equation correlating the peak-to-peak width of EPR spectra of nanoparticles with corresponding particle size: where  is proportionality constant and Δ  is the peak width.Knowing that "" for Ni is 1.6 [13] we can evaluate the particle size of each catalyst.
We correlate the particle size of different ratio of trimetallic catalysts with the function of real synergism Figure 9. where  is particle size by nm and  is the syn/Ni.
By the above function we can deduce the particle size of any trimetallic system by the syn/Ni function.
In order to evaluate the secret of catalytic activity we make pulsed chemisorption to measure the metallic surface area and degree of dispersion and average crystallite size of Ni.This pulsed chemisorption was performed at reaction temperature (Table 2).
Correlating surface area per unite % Ni with the synergism function we can obtain a good correlation as in Figure 10.
We can simulate the previous function as follows:  the trimetallic system.The program written in Excel could simulate easily with a good precision the whole system for the all isomers of nitrophenol (p-, m-, and o-).
Figure 11 is the image of this Excel program which is attached with this manuscript as Supplementary Material (available online at https://doi.org/10.1155/2017/9464209).
Evaluation of the Program.In order to evaluate the program the catalytic activities of some random samples were measured and compared with data given from the program, Table 3.
From Table 3 we could realize that the program could predict precisely the time to reach 100% with average % error not exceeding 1.83% and max error not increasing by 8%.
Mechanism of Reduction of Nitrophenols.Figure 12 represents the kinetic curve of the three nitrophenols which shows that all reduction processes follow first-order equation with respect to nitrophenols.
Durability Study of Some Samples Given.From Table 4 it was found that most of trimetallic system showed a good durability upon five times of successive use.

Conclusions
The trimetallic system (Ni, Cu, and Ag) was proved to be very effective for reduction of nitrophenols using hydrazine hydrate.Thus, nickel metal was found to be the only active center in this reaction.However, the presence of bimetallic and trimetallic alloy with nickel exhibits extra synergism.A program with high precision was successfully written to simulate the synergism function and the catalytic activity of the trimetallic system.The method of calculation used will open new aspect to deal with similar system in the

Figure 6
Figure 6 represents the 3D model of simulated function which gives us high accuracy.Figures5 and 6showed clearly that the 3D model of the linear synergism functions showed more than 3D peaks varied according to the composition of the system.The maximum synergism was found to be of value about 22% syn/Ni.This 3D figure could give us the ability to evaluate the four parameters in one time (% syn, % Ni, % Cu, and % Ag).However, subsequently, analysis of the real system should be compared which will be performed in the next section.After applying the previous function on real trimetallic systems we observed that extra synergism occurs due to

Figure 9 :
Figure 9: Correlation between particle size and function of synergism over % Ni.
=  +  * exp(− exp(−(( −  ln(ln(2)) − )/)where  is the metallic surface area by % Ni and  is the synergism function per % Ni.By the above correlation we can simulate the metallic surface area to the synergism function easily to fully simulate

Figure 11 :
Figure 11: Image of Excel program.

Table 1 :
Examples of some measured bimetallic and trimetallic systems.

Table 2 :
Experimental data for some selected samples for dispersion and metallic surface area.

Table 3 :
A comparison between real and calculated results of some random samples combined with % error.

Table 4 :
A durability study for some trimetallic samples for reduction of PNP.