Abstract
In this work, we build the model and derive the theory of nonlinear deformation for substitutional alloy AB with interstitial atom C and face-centered cubic structure under pressure from the statistical moment method. The calculation results for FCC-AuCuSi are presented. We obtain the values of density of deformation energy, maximum real stress, limit of elastic deformation, and the stress-strain curve and compare the calculated results with experiments and other calculations.
1. Introduction
The study of elastic and thermodynamic properties of perfect ternary and binary interstitial alloys has also been attractive such as the theory of interstitial alloy [1], calculations from first principles, many-body potentials and dynamical dynamics for defects in metals, alloys and solid solutions [2–4] and elastic and thermodynamic properties of perfect ternary and binary interstitial alloys. We have studied the elastic deformation for body-entered cubic (BCC) and face-centered cubic (FCC) ternary and binary interstitial alloys under pressure by the statistical moment method (SMM) in [5–10].
The dependence of elastic and nonlinear deformations of materials on temperature, pressure, and concentration of components has very important role in predicting and understanding their interatomic interactions, strength, mechanical stability, phase transition mechanisms, and dynamical response. Silicides such as AuCuSi have attracted a lot of attention in recent years because of their functional applications and unusual physical properties. Gold silicide or gold silicon is one of the numerous metal alloys sold by American Elements under the trade name AE AlloysTM.
The experimental data on the real stress and the limit of elastic deformation in the nonlinear deformation of pure metal Au are presented in [11–14].
Thermodynamic and mechanical properties of metals and interstitial alloys are studied by some theoretical methods and simulations. For example, Mehl and Papaconstantopoulos [15] applied a tight-binding (TB) scheme to extend first-principle calculations (ab initio) to regimes containing 102–103 atoms in a unit cell and used two-center, nonorthogonal tight-binding parameters and on-site terms. Numerical calculations for many metals are compared with ab initio calculations and experiments. Deformation mechanisms in the mechanical response of nanoporous gold are investigated by molecular dynamics simulations [16]. In addition, in recent years, some researchers have considered factors affecting the structure, the phase transformation, and the crystallization process of alloys AuCu [17], NiCu [18, 19], and AgCu [20].
In the present paper, we will study nonlinear deformation of FCC ternary alloy (substitutional alloy AB with interstitial atoms C) under pressure by the statistical moment method (SMM) [21–24]. In Section 2, we build the model and theoretical calculations, and in Section 3, we carry out numerical calculations for alloy AuCuSi.
2. Model and Theoretical Calculations
In our model, for interstitial alloy AC with FCC structure and concentration condition cC << cA ( is the concentration of atoms A, is the number of atoms A, is the concentration of atoms C, is the number of atoms C, and is the total number of atoms of the alloy AC), the cohesive energy and the alloy parameters (k is called the harmonic parameter and are called anharmonic parameters) for the interstitial atom C in face centers of cubic unit cell in the approximation of two coordination spheres have the form [5–10]
The cohesive energy and the alloy parameters for main metal atom A1 in body center of cubic unit cell in the approximation of three coordination spheres have the form [5–10]
The cohesive energy and the alloy parameters for the main metal atom A2 in corners of cubic unit cell in the approximation of three coordination spheres have the form [5–10]where is the interaction potential between atoms A and C, is the nearest neighbor distance between the atom X (X = A, A1, A2, C) (A in clean metal and A1, A2, and C in interstitial alloy AC) and other atoms at temperature T, is the nearest neighbor distance between the atom X and other atoms at T = 0 K and is determined from the minimum condition of the cohesive energy is the displacement of atom X from equilibrium position at temperature T. are the corresponding quantities in the clean metal A with FCC structure in the approximation of two coordination spheres and have the form [21, 22]
The equations of state for FCC interstitial alloy at temperature T and pressure P and at 0 K and pressure P are written in the form [5, 6, 8–10, 23]
From that, we can calculate the nearest neighbor distance , the parameters , the displacement of atom X from equilibrium position as in [21], the nearest neighbor distance , and the mean nearest neighbor distance between two atoms A in alloy as follows [5, 6, 8–10]:
The Helmholtz free energy of FCC interstitial alloy AC with the condition cC << cA is determined by [5–10,21]where is the Helmholtz free energy of one atom X, U0X is the cohesive energy, and is the configurational entropy of FCC interstitial alloy AC.
The nearest neighbor distances between two atoms in alloy AC after deformation have the form [9]
The mean nearest neighbor distances between two atoms A in interstitial alloy AC at pressure P and temperature T after deformation have the form [9]
The mean nearest neighbor distances between two atoms A in substitutional alloy AB with interstitial atoms C at pressure P and temperature T and at pressure P and temperature T = 0 K after deformation have the form [5, 6, 9]
The Helmholtz free energy of alloy ABC before and after deformation with the condition cC << cB << cA is determined by [5, 6, 9]
The relationship between the stress and the deformation in nonlinear deformation is given by [9]where σ0ABC and αABC are constants for every alloy.
The density of deformation energy can be written in the form [9]
When the deformation rate is constant, the density of deformation energy of alloy is determined by [9]where CABC is the proportional factor. At the maximum value of the density of deformation energy, we have
The maximum value of stress σABCmax and the maximum real stress σ1ABCmax are [9]
From the maximum condition of stress we derive the deformation corresponding to the maximum value of real stress as follows [9]:
CABC is determined from the experimental condition of stress σ0,2ABC in alloy in the form [9]
From the obtained value of , we can calculate σ0ABC and αABC. The limit of elastic deformation of alloy ABC is determined by [9]
Here is Young’s modulus of alloy ABC and has the form [5, 6]
Here are Young’s moduli of alloy AC and metals A, B, respectively.
3. Numerical Results for Alloy AuCuSi
To describe the interaction between atoms Au and Si, we apply the Mie–Lennard-Jones pair interaction potential in the form [25]where D is the depth of potential well corresponding to the equilibrium distance r0 and m and n are determined empirically. Then, the potential parameters for the interaction Au-Si are determined by [26]
We find mAu-Si and nAu-Si by fitting the theoretical result with the experimental data for Young’s modulus of interstitial alloy AuSi3% at room temperature. The Mie–Lennard-Jones potential parameters for the interactions Au-Au, Si-Si, and Au-Si are given in Table 1. The Poisson ratio is 0.42 for Au [28], 0.34 for Cu [29], and 0.21 for Si [30]. Our investigated range of temperature up to 900 K is below the melting temperature of Au, Cu, and Si in the range of pressure from zero to 6 GPa as studies on the melting curve for these materials [31–33].
Fromthe experimental value of Young’s modulus E = 89.1.109 Pa for Au at T = 300 K and P = 0 [34]. From that, we obtain Pa. Figure 1 shows the density of deformation energy and the real stress of AuCuSi at T = 300 K, , cSi = 1%, and cCu = 1, 3, and 5% calculated by the SMM. For AuCuSi at the same temperature, pressure, and concentration of interstitial atoms, when the concentration of substitutional atoms increases, the maximum real stress decreases by 2% and the elastic limit also decreases by 2%. The density of deformation energy has the maximum value when the strain , and from that, we can find the maximum real stress the elastic limit , and the elastic strain for AuCuSi as shown in Table 2.

(a)

(b)
When concentrations we obtain the nonlinear deformation of main metal Au. Figure 2 shows the curve of real stress for Au at and T = 300 K where we have the comparison between the SMM and experiments [11, 12]. The maximum real stress and the elastic limit calculated by the SMM in Table 3 are in good agreement with the molecular dynamics (MD) results [16] and experiments [13, 14]. The error of the maximum real stress between the SMM calculation and the MD result [16] is 0.6%. The error of the elastic limit between the SMM calculation and the experimental data [13] is 5.65%. Thus, the SMM calculations of the maximum real stress and the elastic limit for Au at and T = 300 K are in good agreement with other calculations and experiments.

Figure 3 shows the density of deformation energy and the real stress of AuCuSi at T = 300 K, , cCu = 10%, and cSi = 0, 1, and 2% calculated by the SMM. For AuCuSi at the same temperature, pressure, and concentration of substitutional atoms, when the concentration of interstitial atoms increases, the maximum real stress and the elastic limit decrease strongly. We have the comparison on values of εF, σ1max, εe, and σe for alloys AuSi, AuCu, and AuCuSi in Table 4.

(a)

(b)
Figure 4 shows the density of deformation energy and the real stress of AuCuSi at T = 300, 600, and 900 K, , cCu = 10%, and cSi = 2% calculated by the SMM. For AuCuSi at the same pressure, concentration of substitutional atoms, and concentration of interstitial atoms, when the temperature increases, the maximum real stress and the elastic limit increase. We have the comparison on values of εF, σ1max, εe, and σe for Au88Cu10Si2 at P = 0 and T = 300, 600, and 900 K in Table 5.

(a)

(b)
Figure 5 shows the real stress of AuCuSi at T = 300 K, , 4, and 6 GPa, cCu = 10%, and cSi = 3% calculated by the SMM. For AuCuSi at the same temperature, concentration of substitutional atoms, and concentration of interstitial atoms, when the pressure increases, the maximum real stress increases and the elastic limit decreases. When the pressure increases from 2 to 4 GPa, the maximum real stress increases by 3% and the elastic limit decreases by 7.86%. When the pressure increases from 4 to 6 GPa, the maximum real stress increases by 11% and the elastic limit decreases by 10.56%. We have the comparison on values of εF, σ1max, εe, and σe for Au87Cu10Si3 at , 4, and 6 GPa and T = 300 K in Table 6.

4. Conclusion
In our paper, we derive the analytic expressions of characteristic quantities for the nonlinear deformation such as the density of deformation energy, the maximum real stress, and the limit of elastic deformation depending on temperature, pressure, and concentration of components together with the strain-stress of substitutional alloy AB with interstitial atom C and FCC structure under pressure. We apply theoretical results to alloy AuCuSi. The maximum real stress and the elastic limit calculated by the SMM are in good agreement with MD results [16] and experiments [13, 14]. For AuCuSi at the same temperature, concentration of substitutional atoms, and concentration of interstitial atoms, when the pressure increases, the maximum real stress increases and the elastic limit decreases. When the pressure increases from 2 to 4 GPa, the maximum real stress increases by 3% and the elastic limit decreases by 7.86%. When the pressure increases from 4 to 6 GPa, the maximum real stress increases by 11% and the elastic limit decreases by 10.56%. Our calculated results for alloy AuCuSi are compared with ones for alloys AuCu, AuSi, and metal Au. If we use more coordination spheres and we have exact experimental data for the stress of alloy AuCuSi at different temperatures, pressures, and concentrations of components, we will obtain better calculations.
Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.