THE STUDY OF THE INFLUENCE OF CHEMICAL DISORDER ON MECHANICAL PROPERTIES OF BINARY ALLOYS OF THE TYPE AND (A = Au; Ag) USING COMPUTER MODELING

In this work, computer simulations of binary alloys of nickel with gold and silver were performed using the Sutton-Chen mathematical model. The simulation was conducted using the general lattice program (GULP) software package [9]. For the alloys, a simulation of ordered and chemically disordered structures were carried out. Elastic moduli, Poisson’s ratio, intrinsic hardness and ductility coefficient were estimated for both ordered and disordered structures. The influence of chemical disorder on the mechanical properties of the alloys was estimated. An attempt was made to assess the significance of the ratio of nickel in the alloys to their mechanical properties. The study found that the influence of chemical disorder on all studied mechanical properties of the metals is significant. In particular, the influence is much higher in shear modulus, on the elastic constant , intrinsic hardness and coefficient of ductility.


Introduction
Historically 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.
Воронежский научно-технический вестник № 3(33) сентябрь 2020 г. 5 The purpose of this work is to study mechanical properties of Nickel alloys of the type 3 ANi and 3 NiA 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.

Structure of Metals and Alloys
The metals Ni, Ag and Au have face centered cubic (fcc) while the alloys 3 ANi and 3 NiA may have variety of structures under different set of conditions [4]. The study focuses on the alloys having 2 1 L 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 3 AgNi 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 3 AgNi 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.

Estimation of Elastic Moduli
Elastic 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 Воронежский научно-технический вестник № 3(33) сентябрь 2020 г.
6 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]: 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. 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]: Воронежский научно-технический вестник № 3(33) сентябрь 2020 г.
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]: 11 12 0 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 ij C [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: Воронежский научно-технический вестник № 3(33) сентябрь 2020 г.

1 Poisson's ratio, Intrinsic Hardness and Ductility Coefficient
Poisson'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]: 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: 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]

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 Воронежский научно-технический вестник № 3(33) сентябрь 2020 г. 9 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)  Воронежский научно-технический вестник № 3(33) сентябрь 2020 г. 10 The values of elastic constants ij C 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 3 ANi and 3 NiA 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.     In all the alloy systems for both the ordered and disordered alloys, the elastic constants and moduli of the type 3 AgNi is greater than that of 3 NiAg . The same can be said about the alloys 3 AuNi and 3 NiAu . All the alloy systems investigated here satisfy the stability condition given in equation (4). Further, for all pure metals, ordered and disordered alloys, 12 11 C B C  .
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 3 AgNi is greater than 3 NiAg . The value of hardness for the ordered 3 AuNi is greater than 3 NiAu . For the disordered alloys, the hardness of 3 AuNi is less than 3 NiAu .
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 3 NiAg . For the disordered alloys the reverse is true. For ordered alloys, the ductility coefficient of the alloy 3 AuNi is lower than that of 3 NiAu . Again the opposite is true with regard the ductility coefficient of disordered alloys of 3 AuNi and 3 NiAu .
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 Воронежский научно-технический вестник № 3(33) сентябрь 2020 г. 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 12 C 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 44 C , intrinsic hardness and coefficient of ductility.

Conclusion
In 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 al-