IntroductionHistorically metals have been used as construction materials due to their strength. Some metals and alloys are also used as high temperature materials in industries as gas turbines and heat exchangers. Nickel based alloys are used as gas turbine engines, space vehicles, nuclear reactors and many more [1]. Nickel is a corrosion resistant metal [2]. The Nickel based alloys are frequently used as in condensers and heat exchangers due to their resistance to corrosion and stress-corrosion cracking [3]. An attempt to enhance the specific properties has been done by alloying metals for many years. Nickel could be alloyed with different metals to increase their resistance to corrosion.The purpose of this work is to study mechanical properties of Nickel alloys of the type  and  where A = Au, Ag. The alloys and metals were simulated using general utility lattice program (GULP). In this paper, Sutton-Chen (SC) formalism was used to model the energy of interaction among the atoms of the metals and their alloys. This empirical method is frequently used to model metals and alloys of face centered cubic (fcc) crystal lattice structure. This empirical method is applied for determining the mechanical properties such as elastic constants, bulk and shear moduli, poisson’s ratio, intrinsic hardness and ductility coefficient.1. Structure of Metals and AlloysThe metals Ni, Ag and Au have face centered cubic (fcc) while the alloys  and  may have variety of structures under different set of conditions [4]. The study focuses on the alloys having crystal structures. Here the letters A represent the pure metals Ag and Au. To conduct the simulation, an input file was prepared. An input file contains all the necessary information to conduct simulation as a single text file. The coordinates of the atoms used for simulation are constructed from the basic elementary cell known as unit cell. For example, the unit cell of the ordered alloy has the following basis vectors: Ag: a(0:0; 0:0; 0:0); Ni: a(0:5; 0:5; 0:0); Ni: a(0:0; 0:5; 0:5); Ni: a(0:5; 0:0; 0:5) where a is the edge of the elementary cube. From the unit cell, a supercell was constructed by translating the unit cell in three directions to create a cube of 27 unit cells () with a total of 108 atoms. The same procedure was followed in constructing the supercell for all studied metals and alloys. To simulate the influence of chemical disorder atoms were arranged using python random number program designed solely constructing supercells with chemical disorder. It worth noting that the arrangement of atoms was done in such away that the stoichiometric ratio is not violated. For instance, elementary cell of will have always three atoms of nickel and one atom of silver. It is worth noting that in the ordered structures the elementary cells are equivalent to each other. This may not be the case for the chemically disordered structures due to the random arrangement of atom in the unit cell.2. Estimation of Elastic ModuliElastic constants can provide information on the stability, stiffness, brittleness, ductility, and anisotropy of a material [5, 6]. They are essential parameters that correlate the microscopic properties to macroscopic mechanical properties. Precise calculation of the elastic constants is essential for gaining insight into the mechanical strength of solids, verifying their stability, and designing material applications [7].For a general 3-D materials there are six components of stress and a corresponding six components of strain. Applying Hooke’s law, the relation between stress and strain is expressed as [8]:,                                                    (1) where C is elastic constant,  is stress and   is strain. The equilibrium crystal structure and elastic moduli of the crystals were calculated using the GULP software package [9]. To describe interatomic interactions in metals and their alloys, the Sutton-Chen mathematical model was used [10]. The total potential energy of a system of n interacting atoms in the framework of this model can be found in the article published by Desta O. G., Bykova M. I. and Timoshenko Yu. K. [11, 12]. The Sutton-Chen potential parameters for the metals Au, Ag and Ni are given in table 1.Table 1 – Parameters of the Sutton-Chen potential for metals Au, Ag and NiMetalca(Å)ReferenceNi6939.4323.52[13]Ag612144.414.0900[14]Au81034.4084.0800[14] The calculations were performed at 0 K. To find the equilibrium crystal structure at absolute zero temperature, the total potential energy was minimized by the Newton-Raphson method. Spatial configurations in the zeroth approximation were constructed using data on the crystal structure of metals Ni, Au and Ag, which under normal conditions have a face centered cubic lattice [15].In the GULP formalism, these moduli are estimated from elastic constants. The values of these properties are approximated from semi-empirical potential that represents the interaction energy among the atoms of the metals and their alloys. The second derivative of potential energy function with respect to strain gives elastic constants [10]: ,                                                          (2)where  is a component of the stiffness matrix C, U is the energy expression, V is the volume of the unit cell,  and  are strain. For fcc cubic crystals the only unique elements are the elastic constants,  and . In terms of the interatomic potentials, the equation used to compute the values for these independent elastic constants can be expressed as follows [16]: ; ;   .                          (3) Mechanical stability (structural stiffness) of a material determines how much a material deforms under load. For the cubic crystal, the mechanical stability criteria is [6]:  ;   ; ; .                    (4) Bulk and shear moduli are important material properties. The bulk modulus (B) describes the elastic response to volume change. Conceptually, shear modulus it is the ratio of shear stress to shear strain. For isotropic polycrystalline materials, the bulk modulus (B) and shear modulus (G) can be estimated from elastic constants  [7]. One of the most widely used methods to estimate the elastic characteristics of polycrystalline materials is the averaging of the single crystal elastic moduli, known as the Voigt (V), Reuss (R) and Hill (H)[7]. Hill showed that the Voigt approximation leads to overestimated values of the elastic moduli, while the Reuss approximation underestimated the values [17], and recommended to take the arithmetic mean of value of these approximations [18]. The formulas for estimation of bulk and shear moduli due to the Voigt(V) and Ruess (R) approximations are given as follows: ; .      (5) Here we used the Hill(H) approximations to estimate the bulk and shear moduli of the metals and their alloys. For fcc crystal structures  while the value for shear modulus according to the Hill’s approximation is computed using the formula below:.                                                   (6)3. 1 Poisson’s ratio, Intrinsic Hardness and Ductility CoefficientPoisson’s ratio  defined as the ratio of transverse strain to the longitudinal strain is used to reflect the stability of the material against shear and provides information about the nature of the bonding forces [6]. The Poisson’s ratio, , can be computed using the formalism given as follows from the values of bulk and shear moduli [19, 20]:.                                                   (7)Hardness is one of the most important mechanical properties of metals and alloys, which characterizes its own resistance to deformation under the influence of external forces. It determines the wear resistance of surfaces, is an indicator of the strength of the material [21] and, ultimately, is determined at the atomic level [22].The ductility of metals and alloys is closely related to the hardness of the material. Ductility refers to the ability of a material to deform without breaking. Pugh [23] suggested the coefficient of ductility could be estimated taking the ratio between bulk and shear moduli as follows:.                                                          (8)K indicates ductility or fragility of the material. A large value of K indicates that the material is ductile and a small K describes it is brittle. A given substance will be brittle if K < 1.75 [7]. Otherwise the material is characterized as ductile. There are a number of different empirical formulas for estimating intrinsic hardness of metals and alloys. Here we apply one of the formalisms followed Desta and his co-authors given as [11,12]:.                                              (9) 4. Results and Discussion In all our simulations, we used geometry optimization method by keeping the number of atoms (N) and pressure (P) constant. In addition periodic boundary conditions were applied. This condition makes it possible to simulate a small system that is not terminated by a surface. This is due to the fact that the simulation cell is periodically repeated in all directions. This allows us to consider the systems under study to be macroscopically homogeneous. The simulation was conducted at 0 Kelvin and 0.00 GPa.The Newton-Raphson optimizer with Broyden-Fletcher-Goldfarb-Shanno (BFGS) hessian matrix updater was used in the optimization process. For optimized structure the elastic constants and bulk and shear moduli were calculated. From the results of the simulations, intrinsic hardness for each material was calculated using equation (9) while the coefficient of ductility was estimated using equation (8) for both chemically ordered disordered structures. The calculation results are shown in tables 2 and 3 and in figures 2-7.                              a                                                                                    bFigure 1 – Comparison of bulk moduli of the ordered and disordered alloys: a) Ni -Ag; b) Ni-Au.   a                                                                                bFigure 2 – Comparison of shear moduli of the ordered and disordered alloys:a – Ni -Ag; b – Ni-Au. The values of elastic constants  were estimated using equation (3). The bulk and shear moduli were calculated using the Hill approximation. The Poisson’s ratio for each metal or alloy was calculated using equation (7). The values of intrinsic hardness were calculated for all the pure metals Ni, Ag and Au and as well as for the alloys  and  where the letter A represent the pure metals Ag and Au. Similarly the coefficient of ductility was calculated for all the metals and their alloys under investigation. The calculations of the mechanical properties for the alloys, has been done for both chemically ordered and disordered materials to study the influence of chemical disorder on the properties.a                                                                      bFigure 3 – Comparison of elastic moduli of the ordered and disordered alloys: a – Ni -Ag; b – Ni-Au Table 2 – Chemically ordered alloys: , , , B, G, , in GPa; K and   are dimensionless quantitiesMetalNi230.8340177.136579.6187195.035751.542753.30273.78400.38196.6286140.169675.9039158.989251.067524.09603.11330.35129.438888.956353.6246102.450536.290463.59462.82310.34Ag106.645672.938544.880484.1742130.310933.22412.77700.34211.2794157.335161.2112175.316144.056352.79103.97940.38187.8893149.150146.2503162.063232.619451.75284.96830.40Au179.8710147.789942.1339158.483628.610831.41165.53930.41 a                                                                           bFigure 4 – Comparison of elastic moduli of the ordered and disordered alloys: a – Ni -Ag; b – Ni-Au a                                                                          bFigure 5 – Comparison of elastic moduli of the ordered and disordered alloys: a – Ni -Ag; b – Ni-Au Table 3 – Chemically disordered alloys: , , , B, G, , in GPa; K and  are dimensionless quantitiesMetal157.4414128.350823.4546138.301524.91931.27735.55000.42103.776872.519214.150284.2555514.49950.82625.81090.42195.0908168.958550.3417176.355830.40571.39865.80010.42183.7217151.462045.5370162.150429.71371.47475.45710.41 a                                                                  bFigure 6 – Comparison of intrinsic hardness of the ordered and disordered alloys: a) Ni -Ag; b) Ni-Au In all the alloy systems for both the ordered and disordered alloys, the elastic constants and moduli of the type   is greater than that of  . The same can be said about the alloys  and . All the alloy systems investigated here satisfy the stability condition given in equation (4). Further, for all pure metals, ordered and disordered alloys, .The Poisson’s ratio for the metals and ordered alloys falls in the interval 0:34   0:41 while for the disordered alloys it is approximately 0:41 or 0:42. We were able to estimate the intrinsic hardness of the metals and their alloys. For both ordered and disordered alloys, the hardness values of  is greater than . The value of hardness for the ordered  is greater than . For the disordered alloys, the hardness of  is less than .As it can be seen from tables 2 and 3 the metals and their alloys have good ductility coefficients. The value for ordered alloys, the ductility coefficient, K, of alloy AgNi3 is greater than the ductility coefficient of. For the disordered alloys the reverse is true. For ordered alloys, the ductility coefficient of the alloy  is lower than that of . Again the opposite is true with regard the ductility coefficient of disordered alloys of  and .It is worth noting that the hardness values of the pure metal A is less than the hardness values of the alloys of the type ANi3 and NiA3. The same can be said about the ductility coefficient except for the alloys AuNi3 and NiAu3 where the value is lower a                                                                  bFigure 7 – Comparison of ductility coefficient of the ordered and disordered alloys:a – Ni -Ag; b – Ni-Au than the ductility of pure gold.The influence of chemical disorder is much higher in the Ni-Ag alloy system in comparison to Ni -Au. In general terms, the values of the elastic constants and moduli of the disordered alloys are less than the ordered counterparts except for the elastic constant  for the alloy system Ni-Au. In all the materials, intrinsic hardness decreased with disorder while coefficient of ductility increased. The highest influenced mechanical properties by chemically disorder are the shear modulus (G), the elastic constant , intrinsic hardness and coefficient of ductility.ConclusionIn the present work, we were able to estimate bulk and shear moduli, elastic constants, Poisson’s ratio, intrinsic hardness and ductility coefficients of alloys using computer simulation. Further the extent of influence of chemical disorder on the mechanical properties of the alloys was estimated.For binary ordered alloy of nickel with gold or silver, the alloy will have higher values of  the elastic moduli and intrinsic hardness in comparison the pure metals Au or Ag. The higher the ratio of nickel in the alloy, the higher the values of the elastic moduli and intrinsic hardness is.Generally, chemical disorder in the alloys decreases the values of elastic constants and moduli and intrinsic hardness while it increases the values of ductility coefficient and Poisson’s ratio. The influence of chemical disorder is much higher in the alloys involving nickel and silver. Chemical disorder influences higher the shear modulus, elastic constant , intrinsic hardness and coefficient of ductility.