Next Article in Journal
Three-Dimensional Vibration Analysis of Laminated Composite Rectangular Plate with Cutouts
Previous Article in Journal
Europium Doping Impact on the Properties of MBE Grown Bi2Te3 Thin Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First-Principles Study on the Cu/Fe Interface Properties of Ternary Cu-Fe-X Alloys

1
State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai 200240, China
2
Shanghai Key Laboratory of Advanced High-Temperature Materials and Precision Forming, Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(14), 3112; https://doi.org/10.3390/ma13143112
Submission received: 28 May 2020 / Revised: 3 July 2020 / Accepted: 7 July 2020 / Published: 13 July 2020
(This article belongs to the Section Materials Physics)

Abstract

:
The supersaturated Fe in Cu is known to reduce the electrical conductivity of Cu severely. However, the precipitation kinetics of Fe from Cu are sluggish. Alloying is one of the effective ways to accelerate the aging precipitation of Cu-Fe alloys. Nucleation plays an important role in the early stage of aging. The interface property of Cu/γ–Fe is a key parameter in understanding the nucleation mechanism of γ-Fe, which can be obviously affected with the addition of alloying elements. In this paper, first principles calculations were carried out to investigate the influence of alloying elements on the interface properties, including the geometric optimizations, interfacial energy, work of adhesion and electronic structure. Based on the previous research, 14 elements including B, Si, P, Al, Ge, S, Mg, Ag, Cd, Sn, In, Sb, Zr and Bi were selected for investigation. Results showed that all these alloying elements tend to concentrate in the Cu matrix with the specific substitution position of the atoms determined by the binding energy between Fe and alloy element (X). The bonding strength of the Cu/γ-Fe interface will decrease obviously after adding Ag, Mg and Cd, while a drop in interfacial energy of Cu/γ–Fe will happen when alloyed with Al, B, S, P, Si, Ge, Sn, Zr, Bi, Sb and In. Further study of the electronic structure found that Al and Zr were not effective alloying elements.

1. Introduction

Supersaturated Fe in a Cu matrix is known to reduce the electrical conductivity of Cu severely due to the slow precipitation kinetics of Fe [1,2,3]. Adding a third element allows a larger variety of possible kinetic paths during aging precipitation to improve the strength and conductivity of the binary Cu-Fe alloys [4,5]. The aging precipitation process of Cu-Fe alloy includes nucleation, growth and coarsening. However, few works focusing on the effect of the alloying elements on the nucleation of γ-Fe can be found. In the past, experimental investigations [6,7] of γ-Fe precipitates in Cu reveal that the antiferromagnetic behavior of γ-Fe precipitates at low temperatures has been well established. According to the order of magnetic moments of γ-Fe precipitates, the interface orientation relationship between γ-Fe precipitated and Cu matrix was (100)[010]Cu//(100)[010]γ–Fe, while confirming the absence of inter diffusion with surface alloy formation. The interface properties between the Cu matrix and γ-Fe, especially the atomic structure and interfacial energy, are crucial to make clear the nucleation mechanism of the γ-Fe in ternary Cu-Fe-X alloys. However, solid/solid interfacial adhesion and energy are difficult to be obtained by experiment. Recently, the first-principles calculation has been successfully used in the detailed study of metal/ceramic interface adhesion of Al-base [8], Mg [9] Fe [10] and Ni [11]. In the present work, the first principles calculation will be employed to study the Cu/γ-Fe interfacial properties to make clear the effect of alloying elements on the nucleation of γ-Fe in the early stage of aging of the Cu-Fe alloy. Based on the previous research [12,13], a large database of solute–vacancy binding energies in Cu was calculated, and then the solubility and diffusivities of iron in dilute Cu-Fe-X alloys were studied. These studies provide a basis for the selection of alloying elements; therefore, 14 elements including B, Si, P, Al, Ge, S, Mg, Ag, Cd, Sn, In, Sb, Zr and Bi will be selected for this investigation. Calculation on the substitutional energies, lattice misfits, work of adhesion, interfacial energy and electronic structure will be carried out with reasonably built Cu/γ-Fe interface configuration.

2. Models and Computation Method

Calculations were performed using the pseudopotential plane wave method within the density functional theory (DFT) T = 0 K as implemented in the VASP package [14,15,16,17]. The projector augmented wave (PAW) potentials were used with the generalized gradient approximation (GGA) [18,19,20,21] for exchange–correlation. The electron–ion interaction is described by the 3d, 4s states for Cu, the 3d, 4s states for Fe and the corresponding valence electron for the third elements. The integration over the Brillouin zone integrations was performed with a 6 × 6 × 1 Monkhorst–Pack k-point grid in case of interface and surface structures; a grid of 16 × 16 × 16 was used in case of bulk structures. An energy cutoff of 340 eV was used for the plane wave basis. The convergence of results with respect to energy cutoff and k-points was carefully considered. In order to obtain a precise crystal structure, atomic positions as well as unit cell volume and shape were fully relaxed. Ground-state geometries were determined by minimizing stresses and Hellmann–Feynman forces with the conjugate gradient algorithm, until the forces on all atomic sites were less than 10−3 eV/Å. Convergence of the electronic system was reached when the difference of the total energies between two iterations was below 10−4 eV. For magnetic 3d elements such as Fe, spin-polarized calculations were performed to analyze the influence of magnetic moments on the energy and electronic structure, γ–Fe was spin-polarized (antiferromagnetic ordering).

3. Results and Discussion

3.1. Bulk and Surface Calculation

The calculations on bulk properties of Cu and γ–Fe were carried out by GGA. The calculated lattice constant of Cu and γ–Fe were a = 3.634 Å and a = 3.447 Å, respectively, which were consistent with the experimental value and other calculation results. Table 1 shows the calculated lattice constants and cell volumes.
The free surfaces of the γ-Fe and Cu slabs with the in-plane periodicity were separated by a 15 Å vacuum to prohibit their interactions. Based on this model, the substitutional energies, lattice misfits, work of adhesion, interfacial energy and electronic structure were calculated, respectively. To make sure that both sides of the surface slabs of Cu (001) and γ–Fe (001) are thick enough to show the bulk-like interior, the convergence tests on Cu (001) and γ–Fe (001) surfaces with different thickness were carried out first. The calculated surface energies of Cu(100) and γ-Fe(100) slabs were obtained with the following equation [26,27]:
d E = T d S P d V + σ d A + i = 1 C μ i d n i
where E is the total energy of the system, T is the temperature, P is the pressure, V is volume, σ is the surface energy, A is the surface area, C is the component fraction, μ i is the chemical potential of the component and n i is the stoichiometry of the component. Because of the different substitutional sites of Cu/γ–Fe interface, either the X-occupied Fe site or the X-occupied Cu site, the chemical potentials were chosen carefully. Two possible reference states were considered: (1) a dilute solid solution CunX (n = 107); and (2) a dilute solid solution FenX (n = 107). For the X-occupied Cu site, the chemical potentials were calculated from the dilute solution Cu107X and pure Cu: μ X = E Cu 107 X 107 E Cu . For the X-occupied Fe site, the chemical potentials were calculated from the dilute solution Cu107X and pure Cu: μ X = E Fe 107 X 107 E Fe . Table 2 shows the calculated chemical potentials.
Scaled to 0 K, the σ was given by:
σ = E A S , V , n i
For FCC Cu(100) and γ-Fe(100) slabs, the σ is defined as:
σ = 1 A s u r f a c e E s l a b t o t a l n i μ i
E s l a b t o t a l is the total energy of the surface model. In the calculated surface energies list in Table 3, it was found that the surface energies of the Cu(100) slabs and γ–Fe(100) slabs with more than six layers can converge to 1.48 J/m2 and 3.37 J/m2, respectively, which were in agreement with the experimental value and other calculated values.

3.2. Model Geometry of Interface

The interface properties are directly affected by the position of alloying elements. Therefore, the possible positions of alloy elements in Cu(100)/γ–Fe(100) interface were considered firstly in the calculations. The substitutional energy of alloy elements at different lattice positions in the interface model was calculated to determine the occupying tendency; the corresponding atomic position was shown in Figure 1.
The substitutional energy, E Cu X f c c   Cu in Cu and E Fe X f c c   Fe in Fe are defined as follows [6]:
E Cu X f c c   Cu = E Cu n 1 Fe n X t o t a l + μ Cu E Cu n Fe n t o t a l + μ X
E Fe X f c c   Fe = E Cu n Fe n 1 X t o t a l + μ Fe E Cu n Fe n t o t a l + μ X
where E t o t a l is the calculated total energy before and after substitution, n is the number of substitutional atoms and X is the alloy element. Table 4 shows the substitution energy of alloy elements at different lattice positions in the interface model.
The results showed that all of the alloying elements were concentrated in the Cu-based near the interface during the nucleation. The substitution positions of Ag and Mg atoms were concentrated in the 1 position and far away from the interface. The substitution position of the Cd atom was concentrated in the 3 position, which is the second nearest neighbor interface. The other atoms such as Al, B and Bi were concentrated in the 4 position, which is the nearest neighbor interface. The occupation law could be explained by the bonding strength between Fe and the alloy element (X). The calculation of the X–Fe binding energy studied in previous research [9] showed that the elements with strong X-Fe bonding were more likely to be concentrated near the interface. Therefore, ten atoms there (Al, B, Ge, In, P, S, Si, Zr, Bi and Sb) with a strong binding energy tended to be located in the nearest neighbor interface.

3.3. Work of Adhesion and Interfacial Energy

In the study of Finnis [29], the ideal work of adhesion Wad is defined as the bond energy needed to separate an interface into two free slabs; it is an important and convenient factor to predict the mechanical properties and the chemical bonding strength of an interface [30,31,32]. Wad was calculated by the following equation:
W a d = 1 A i n t e r f a c e E Cu t o t a l + E Fe t o t a l E Cu / γ Fe t a t a l
where E Cu / γ Fe t a t a l is the total energy of interface model, and E Cu t a t a l and E Fe t o t a l are the free surfaces energies of Cu and Fe, respectively.
Table 5 shows the effect of alloying elements on the work of adhesion Wad, the magnetic moment and the interfacial distance of fully relaxed interfaces.
The work of adhesion of binary Cu/γ-Fe is 3.822 J/m2, which is in agreement with the other calculated value, 3.89 J/m2 [20]. Furthermore, the work of adhesion of Cu/γ-Fe decreased by the addition of Ag, Bi, Cd, Ge, In, Mg, S, Sb and Sn, increased by the addition of Al, B, P, Si and Zr. Generally, lattice mismatch is an important reason for the change of work of adhesion, and the lattice mismatch at the interface of Cu/γ-Fe is determined by the volume of the alloy atoms. In order to get the influence of lattice misfits on work of adhesion, the volume of the solute atom V S o l u t e X was given by the volume difference induced by placing a single solute into pure Cu:
V S o l u t e X = V Cu n 1 , X V Cu n
where V Cu n 1 , X and V Cu n is the calculated volume before and after substitution. Then, the lattice misfit is given as:
δ = a Cu a γ Fe / a γ Fe
where a is the lattice constant.
The results list in Figure 2 shows that there was a positive correlation between the work of adhesion and the lattice misfits.
It means that the significant strain and lattice misfits can be reduced by adding relatively small atoms, thus, the work of adhesion of the interface becomes larger and the stable Cu/γ-Fe interface is formed. The lattice misfits at the interface of the Cu/γ-Fe effect by alloy elements vary from 5.3% to 6.1%, except for that of the Bi atom. However, the content of alloying elements is lower in reality, which indicates that the addition of alloying elements will not change the coherent relation of the Cu/γ-Fe interface. In addition, the magnetic moment and interfacial distance effect by alloying elements were given in Table 4. Under the influence of alloying elements, the magnetic moment of Fe varies from 2.367 μ B to 2.528 μ B . The interfacial distance of binary Cu/γ-Fe is 1.83 Å, which is in agreement with another calculated value, 2.55 Å [20]. As the volume of alloying element increases, the interfacial distance increases, which should be due to the lattice distortion caused by the addition of a large atom. According to the thermodynamic theory of nucleation, interfacial energy as nucleation resistance is an important factor that hinders γ-Fe nucleation by influencing nucleation work and a critical nucleation radius. The effect of alloy elements on γ-Fe nucleation was determined by calculating the interfacial energy of Cu/γ–Fe. The interfacial energy can be given by:
γ = 1 A i n t e r f a c e E Cu / γ Fe t a t a l i = 1 C μ i d n i σ Cu σ Fe
where A is the interface area, E Cu / γ Fe t a t a l is the total energy of a fully relaxed interface supercell, and σCu and σFe are the surface energy of the Cu(100) and γ-Fe(100), respectively. Table 6 shows the interfacial energy of Cu/γ-Fe effect by alloy elements.
The interfacial energy of Cu/γ-Fe in the binary Cu-Fe alloy is 0.819 J/m2. Furthermore, the interfacial energy of Cu/γ-Fe was enhanced with the addition of Ag, Mg and Cd. On the contrary, the interfacial energy decreased with the addition of the elements of Al, B, Bi, Ge, In, P, S, Sb, Si, Sn and Zr. It was known that the smaller the interfacial energy is, the more stable the interface structure is, the smaller the energy required to form the Cu/γ-Fe interface is, and then the more favorable to the nucleation of Fe it is.

3.4. Electronic Structure

The interface mechanical strength is closely related to the interfacial atomic bonding. In order to explore the interfacial bonding between Cu(100) and γ-Fe(100), the electron structure studies, including the difference charge density and the density of states, were calculated to obtain the bonding of atoms at the interface. The difference charge density Δ ρ was given by the following equation:
Δ ρ = ρ Cu / γ Fe ρ Cu 100 ρ γ Fe 100
where ρ Cu / γ Fe is the total charge density of Cu/γ–Fe interface, and ρ Cu 100 and ρ γ Fe 100 are the charge densities of isolated Cu(100) and the Fe(100) surface, respectively.
Figure 3a displays the difference charge density for the Cu/γ–Fe interface; three typical elements with different relative positions were also given, as shown in Figure 3 The figures show the charge density difference between the Cu/γ–Fe system and that of the ρ γ Fe 100 system and the ρ γ Fe 100 system. It is obvious that the perturbation created by the Fe-Cu binding was mostly localized between Cu and γ–Fe interface, which indicates that the bonding is characterized by covalent and will not be changed by the addition of the third elements. However, the binding strength of Cu/γ–Fe interface was significantly changed by the addition of the third elements. These variations of the binding strength were mainly caused by charge transfer. It is found that the charge density at the interface decreases obviously and the bonding between atoms decreases with the addition of alloy elements Ag, Cd and Mg, as illustrated in Figure 3b. The charge density at the interface increases slightly and the bonding between atoms increases with the addition of alloy elements B, S, Si, Zr, Al, Bi, Ge, In, P, Sb and Sn, as illustrated in Figure 3c.
To investigate the electron population and bonding strength, local density of states (LDOS) only containing the atoms near the interface was further calculated. Figure 4 shows the LDOS of Cu/γ–Fe interface, the difference ΔN between the LDOS of Cu/γ–Fe system and that of Cu(100) system and Fe(100) system revealed clearly the Fe-Cu bonding and antibonding states, which were just below and above the Fermi level: bonding states were around −5.43~−3.54 eV below and antibonding states were around 1.64~3.25 eV above. For LDOS, the energy range from −3.54 to 1.64 eV and the peaks of the Cu/γ–Fe system coincide with the peaks of LDOS in the inner shell of Cu and γ–Fe, which means that they are atomic ground state electrons. The positive and negative fluctuations in the ΔN indicate electron migration. For the bonding states, which were mainly composed of 3d states of Cu and 3d states of Fe, one peak appeared at −4.21 eV. Obviously, a covalent band formed at this coverage and the calculated charge number of the bonding states is 10.79 (charge numbers are the result of integration of the curves).
For ternary alloy, Cu-Fe-Ag(Si) were taken as an example; the Fe-Cu bonding and antibonding states were determined by the difference ΔN between the LDOS of the Cu-Ag(Si)/γ–Fe system and that of the Cu-Ag(Si) (100) system and Fe(100) system. The Cu-Ag/γ–Fe interface list in Figure 5 shows that the bonding states were around −5.45~−3.51 eV above the Fermi level and the antibonding states were around 1.69~3.29 eV.
Compared with the LDOS analysis of the binary Cu/γ–Fe interface, in the ternary Cu-Ag/γ–Fe interface, the calculated charge number of the bonding states decreased from 10.79 to 10.00, indicating that the less charge transferred from Fe (Cu) to bonding states and the bonding strength decreases.
For the Cu-Si/γ–Fe interface list in Figure 6, the bonding states around −7.67~−3.61 eV was above the Fermi level and antibonding states were around 1.62~3.20 eV.
Compared with the LDOS analysis of the binary Cu/γ–Fe interface, in the ternary Cu-Si/γ–Fe interface, the calculated charge number of the bonding states increased from 10.79 to 12.10, indicating that more charge was transferred from Fe (Cu) to bonding states and that the bonding strength increased. More computational results about the charge number of the bonding states affected by alloying elements were listed in Table 7. The calculation revealed that 5 elements Al, Ag, Mg, Cd and Zr can reduce the charge number of the bonding states, which suggests that the bonding strength will decrease by the addition of these elements. The weakening of interatomic bond at the interface indicates that the higher the energy required to form the interface, the higher the interface energy is.
In summary, according to the analysis of interfacial energy, eleven atoms Al, B, Bi, Ge, In, P, S, Sb, Si, Sn and Zr decrease the interfacial energy, and the addition of large solute atoms increases the interfacial distance. However, by studying the electronic structure, it was found that the interfacial distance should be due to the lattice distortion caused by the addition of a large atom. Although the interfacial distance increases, because the solute atoms participated in the bonding at the interface, the charge numbers of the bonding states were enhanced and the interface energies were reduced. Therefore, it was found that nine alloy elements will promote the precipitation of γ-Fe and form stable γ-Fe phases during nucleation.

4. Conclusions

The influence of the nucleation of Fe was investigated by studying the Cu/γ–Fe interfacial properties. The substitutional energy of alloy elements at different lattice positions of the interface were calculated to determine the occupying tendency. The results showed that the all of the alloying elements were concentrated in the Cu-based, the substitution position of Ag and Mg atoms are concentrated in the most far away from the interface, Cd atom is concentrated in the secondary adjacent interface position, and the rest of the atoms were concentrated in the nearest neighbor interface. Under the condition of coherent interface, the interfacial energy was dominated by the work of adhesion, so the distribution of interfacial energy was proportional to the size of atoms. However, under the influence of the bonding strength at the interface, 3 alloy elements—Ag, Mg and Cd—obviously increase the interface energy. The addition of 11 elements—Al, B, S, P, Si, Ge, Sn, Zr, Bi, Sb and In—can reduce the interfacial energy of Cu/γ –Fe, while among those, Al and Zr can reduce the charge number of the bonding states. Therefore, these nine alloy elements will promote the precipitation of γ-Fe and form stable γ-Fe phases during nucleation.

Author Contributions

Y.W. and M.L. conceived and designed the DFT calculations; Y.W. designed and performed the DFT calculations; H.G. performed the surface and interface analysis, J.W. and B.S. initiated and supervised the project. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by the National Natural Science Foundation of China (Grant No. 51671131), and the calculations were performed on the Dawning 5000A supercomputer at the Shanghai Supercomputer Center.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hong, S.I.; Song, J.S. Strength and conductivity of Cu-9Fe-1.2X (X = Ag or Cr) filamentary microcomposite wires. Metall. Mater. Trans. A 2001, 32, 985–991. [Google Scholar] [CrossRef]
  2. Hong, S.I. Copper-Iron Filamentary Microcomposites. Adv. Eng. Mater. 2001, 3, 475–479. [Google Scholar] [CrossRef]
  3. Gao, H.; Wang, J.; Sun, B. Effect of Ag on the thermal stability of deformation processed Cu-Fe in situ composites. J. Alloys Compd. 2009, 469, 580–586. [Google Scholar] [CrossRef]
  4. Gao, H.; Wang, J.; Da, S.; Sun, B. Microstructure and strength of Cu-Fe-Ag in situ composites. Mater. Sci. Eng. A 2007, 452, 367–373. [Google Scholar] [CrossRef]
  5. Gao, H.; Wang, J.; Da, S.; Sun, B. Microstructure and properties of Cu-11Fe-6Ag in situ composite after thermo-mechanical treatments. J. Alloy. Compd. 2007, 438, 268–273. [Google Scholar] [CrossRef]
  6. Fu, C.L.; Freeman, A.J. Electronic and magnetic properties of the fcc Fe(001) thin films: Fe/Cu(001) and Cu/Fe/Cu(001). Phys. Rev. B 1987, 35, 925–932. [Google Scholar] [CrossRef]
  7. Polzonetti, G.; Di Castro, V.; Furlani, C. The Growth Mode of Copper and Iron for the Cu/Fe and Fe/Cu Interfaces. Surf. Interface Anal. 1994, 22, 211–213. [Google Scholar] [CrossRef]
  8. Mao, Z.; Chen, W.; Seidman, D.N.; Wolverton, C. First-principles study of the nucleation and stability of ordered precipitates in ternary Al-Sc-Li alloys. Acta Mater. 2011, 59, 3012–3023. [Google Scholar] [CrossRef]
  9. Wang, F.; Li, K.; Zhou, N.G. First-principles calculations on Mg/Al2CO interfaces. Appl. Surf. Sci. 2013, 285, 879–884. [Google Scholar] [CrossRef]
  10. Wang, J.; Enomoto, M.; Shang, C. First-principles study on the interfacial segregation at coherent Cu precipitate/Fe matrix interface. Scripta Mater. 2020, 185, 425–446. [Google Scholar] [CrossRef]
  11. Zhao, Y.; Wen, Z.; Hua, H.; Wei, G.; Han, P. Density functional theory study of the interfacial properties of Ni/Ni 3 Si eutectic alloy. Appl. Surface Sci. 2014, 303, 205–209. [Google Scholar] [CrossRef]
  12. Wang, Y.; Gao, H.; Han, Y.; Dai, Y.; Wang, J.; Sun, B. Role of the third element in accelerating Fe diffusivities in Cu from first principles. J. Alloys Compd. 2015, 639, 642–647. [Google Scholar] [CrossRef]
  13. Wang, Y.; Gao, H.; Han, Y.; Dai, Y.; Wang, J.; Sun, B. First-principles study on the solubility of iron in dilute Cu-Fe-X alloys. J. Alloys Compd. 2016, 691, 992–996. [Google Scholar] [CrossRef]
  14. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  15. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar]
  16. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  17. Perdew, J.P.; Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 1992, 45, 13244. [Google Scholar] [CrossRef]
  18. Liu, L.; Wang, R.; Xiaozhi, W.U.; Gan, L.; Wei, Q. Temperature effects on the generalized planar fault energies and twinnabilities of Al, Ni and Cu: First principles calculations. Comput. Mater. Sci. 2014, 88, 124–130. [Google Scholar] [CrossRef]
  19. Abhilash, R.; Gian, B.; Guido, F. Core Level Spectra of Organic Molecules Adsorbed on Graphene. Materials 2018, 11, 518–532. [Google Scholar]
  20. Herrig, F.; Music, D.; Völker, B.; Hans, M. Ab Initio Guided Low Temperature Synthesis Strategy for Smooth Face-Centred Cubic FeMn Thin Films. Materials 2018, 8, 384. [Google Scholar] [CrossRef] [Green Version]
  21. Pang, H.; Li, M.; Gao, C.; Huang, H. Phase Transition of Single-Layer Molybdenum Disulfide Nanosheets under Mechanical Loading Based on Molecular Dynamics Simulations. Materials 2018, 11, 502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Benali, A.; Lacaze-Dufaure, C.; Morillo, J. Density functional study of copper segregation in aluminum. Surf. Sci. 2011, 605, 341–350. [Google Scholar] [CrossRef] [Green Version]
  23. Gale, W.F.; Totemeir, T.C. Smithells Metals Reference Book, 8th ed.; Elsevier: Oxford, UK, 2004. [Google Scholar]
  24. Yu, J.; Xin, L.; Wang, J.; Jing, C.; Huang, W. First-principles study of the relaxation and energy of bcc-Fe, fcc-Fe and AISI-304 stainless steel surfaces. Appl. Surf. Sci. 2009, 255, 9032–9039. [Google Scholar] [CrossRef]
  25. Acet, M.; Auml, H.; Wassermann, H.F.; Pepperhoff, W. High-temperature moment-volume instability and anti-Invar of gamma-Fe. Phys. Rev. B 1994, 49, 6012. [Google Scholar] [CrossRef] [PubMed]
  26. Jiang, Y.; Adams, J.B.; Van Schilfgaarde, M. Density-functional calculation of CeO2 surfaces and prediction of effects of oxygen partial pressure and temperature on stabilities. J. Chem. Phys. 2005, 123, 064701. [Google Scholar] [CrossRef]
  27. Han, H.; Yin, G.; Wang, H.; Wang, C.; Shao, K.; Zhang, W.; Dai, J.; Ping, H. First-principles investigation on the geometries, stabilities and defective properties of fluoride surfaces. Comput. Mater. Sci. 2017, 133, 159–166. [Google Scholar] [CrossRef]
  28. Vitos, L.; Ruban, A.; Skriver, H.; Kollar, J. The surface energy of metals. Surf. Sci. 1998, 411, 186–202. [Google Scholar] [CrossRef]
  29. Finnis, M.W. The theory of metal-ceramic interfaces. J. Phys. Condens. Matter 1996, 8, 5811. [Google Scholar] [CrossRef]
  30. Liu, R.; Yin, X.; Feng, K.; Xu, R. First-principles calculations on Mg/TiB2 interfaces. Comput. Mater. Sci. 2018, 149, 373–378. [Google Scholar] [CrossRef]
  31. Hashibon, A.; Elsässer, C.; Mishin, Y.; Gumbsch, P. First-principles study of thermodynamical and mechanical stabilities of thin copper film on tantalum. Phys. Rev. B 2007, 76, 245434. [Google Scholar] [CrossRef]
  32. Chunmei, L.I.; Chen, Z.Q.; Zeng, S.M.; Cheng, N.P.; Chen, T.X. Intermetallic phase formation and evolution during homogenization and solution in Al-Zn-Mg-Cu alloys. Sci. China 2013, 56, 2827–2838. [Google Scholar]
Figure 1. Structure of Cu(100)/γ–Fe(100) interface. The different layers of the interface area are marked with 1–8.
Figure 1. Structure of Cu(100)/γ–Fe(100) interface. The different layers of the interface area are marked with 1–8.
Materials 13 03112 g001
Figure 2. Work of adhesion and lattice misfit of Cu/γ–Fe interface as a function of solute volume.
Figure 2. Work of adhesion and lattice misfit of Cu/γ–Fe interface as a function of solute volume.
Materials 13 03112 g002
Figure 3. Difference charge density (e/Å3) of Cu/γ–Fe interface. (ac) represent no solute atom, a solute atom Ag and Si, respectively. The dashed lines indicate the location of the interfaces. Areas of electron accumulation (red) and depletion (blue) have positive and negative signs, respectively; contour lines are drawn at 0.002 e/Å3 intervals.
Figure 3. Difference charge density (e/Å3) of Cu/γ–Fe interface. (ac) represent no solute atom, a solute atom Ag and Si, respectively. The dashed lines indicate the location of the interfaces. Areas of electron accumulation (red) and depletion (blue) have positive and negative signs, respectively; contour lines are drawn at 0.002 e/Å3 intervals.
Materials 13 03112 g003
Figure 4. The LDOS of the binary Cu-Fe system, the Fermi energy is set to 0 eV. ΔN is the difference between the LDOS of Cu/γ–Fe system and that of Cu(100) system and Fe(100) system.
Figure 4. The LDOS of the binary Cu-Fe system, the Fermi energy is set to 0 eV. ΔN is the difference between the LDOS of Cu/γ–Fe system and that of Cu(100) system and Fe(100) system.
Materials 13 03112 g004
Figure 5. The LDOS of the binary Cu-Fe-Ag system; the Fermi energy is set to 0 eV. ΔN is the difference between the LDOS of the Cu-Ag/γ–Fe system and that of the Cu-Ag(100) system and the Fe(100) system.
Figure 5. The LDOS of the binary Cu-Fe-Ag system; the Fermi energy is set to 0 eV. ΔN is the difference between the LDOS of the Cu-Ag/γ–Fe system and that of the Cu-Ag(100) system and the Fe(100) system.
Materials 13 03112 g005
Figure 6. The LDOS of the binary Cu-Fe-Si system; the Fermi energy is set to 0 eV. ΔN is the difference between the LDOS of Cu-Si/γ–Fe system and that of Cu-Si(100) system and Fe(100) system.
Figure 6. The LDOS of the binary Cu-Fe-Si system; the Fermi energy is set to 0 eV. ΔN is the difference between the LDOS of Cu-Si/γ–Fe system and that of Cu-Si(100) system and Fe(100) system.
Materials 13 03112 g006
Table 1. The lattice constants and cell volumes.
Table 1. The lattice constants and cell volumes.
ParameterCuγ–Fe
GGA
(this work)
GGA
[18]
GGA
[22]
Exp.
[23]
GGA
(this work)
GGA-PBE [24]GGA-PW91
[24]
Exp.
[25]
a 3.633.633.643.623.453.473.473.65
V3·cell−148.0047.9948.2347.2440.9541.9341.8548.79
Table 2. Chemical potentials of elements (X), a dilute solid solution Cu107X and a dilute solid solution Fe107X.
Table 2. Chemical potentials of elements (X), a dilute solid solution Cu107X and a dilute solid solution Fe107X.
Elements (X) μ X ( eV )
Cu107XFe107X
B4.8027.103
Si5.6597.235
P5.2586.595
Al4.4474.691
Ge4.5145.605
S3.5294.261
Mg1.730.911
Ag2.2941.296
Cd0.320.44
Sn3.563.416
In2.2861.637
Sb3.6183.742
Zr8.1438.071
Bi2.4591.877
Table 3. Surface energies of Cu(100) and γ–Fe(100).
Table 3. Surface energies of Cu(100) and γ–Fe(100).
Layer Number/nσ/J·m−2
Cu(100)γ–Fe(100)
41.453.30
61.483.36
81.483.37
101.473.36
121.463.38
Ref1.485, 1.532 [22] 1.79 J/m2 [28]2.956, 2.973 [14]
Table 4. Substitutional energies (eV) of alloy atoms at different positions. Figure 1 lists the 1–8 representing the different substitutional sites. Underline indicates the smallest substitutional energy.
Table 4. Substitutional energies (eV) of alloy atoms at different positions. Figure 1 lists the 1–8 representing the different substitutional sites. Underline indicates the smallest substitutional energy.
Elements12345678
B1.4311.4391.7510.1880.2760.6630.6160.666
Si−0.240−0.284−0.312−1.165−1.076−0.928−1.099−1.052
P0.1070.0330.010−1.065−0.954−0.784−0.953−0.915
Al−0.763−0.777−0.792−1.235−0.817−0.414−0.621−0.585
Ge0.0840.0390.012−0.673−0.384−0.038−0.251−0.235
S0.6350.5350.769−0.5440.0480.6960.4970.516
Mg−0.121−0.119−0.1110.0000.6851.4381.1711.177
Ag0.5950.6020.6040.7441.5642.4502.2932.258
Cd0.6490.6480.6420.6431.5362.5032.2122.176
Sn0.5180.4760.459−0.0500.5731.2250.8540.862
In0.5030.4840.4840.2551.0421.8771.5381.514
Sb0.8780.8000.7660.0630.5641.0870.7400.727
Zr0.5360.4530.495−0.3360.3280.9660.5980.620
Bi1.9161.8441.8181.3492.0552.8632.4092.395
Table 5. Effects of alloying elements on the work of adhesion, the magnetic moment and the interfacial distance of fully relaxed interfaces.
Table 5. Effects of alloying elements on the work of adhesion, the magnetic moment and the interfacial distance of fully relaxed interfaces.
Elements W a d ( J · m 2 ) Magnetic   Moment   ( μ B ) Interfacial Distance (Å)
Cu3.8222.4411.827
B4.185 2.4091.827
Si4.0452.4181.825
P3.8562.4771.827
Al3.9882.4521.831
Ge3.7932.3671.832
S3.6042.4061.838
Mg3.7862.4791.827
Ag3.7752.5021.827
Cd3.7862.4511.837
Sn3.4882.5071.847
In3.4812.5101.849
Sb3.3942.4941.847
Zr3.9722.4161.844
Bi3.0272.5281.855
Table 6. Effects of alloying elements on interfacial energies of Cu/γ–Fe interface.
Table 6. Effects of alloying elements on interfacial energies of Cu/γ–Fe interface.
ElementCuAgAlBBiCdGeIn
γ /J·m−20.8190.8520.6590.2880.7500.8690.5820.641
ElementMgPSSbSiSnZr
γ /J·m−20.8450.4430.4740.5840.5030.6840.579
Table 7. The charge number of the bonding states affected by alloying elements.
Table 7. The charge number of the bonding states affected by alloying elements.
ElementCuAgAlBBiCdGeIn
10.7910.009.7110.9212.8410.6913.2711.94
ElementMgPSSbSiSnZr
9.7412.4012.1713.1212.1012.628.18

Share and Cite

MDPI and ACS Style

Wang, Y.; Li, M.; Gao, H.; Wang, J.; Sun, B. First-Principles Study on the Cu/Fe Interface Properties of Ternary Cu-Fe-X Alloys. Materials 2020, 13, 3112. https://doi.org/10.3390/ma13143112

AMA Style

Wang Y, Li M, Gao H, Wang J, Sun B. First-Principles Study on the Cu/Fe Interface Properties of Ternary Cu-Fe-X Alloys. Materials. 2020; 13(14):3112. https://doi.org/10.3390/ma13143112

Chicago/Turabian Style

Wang, Yufei, Min Li, Haiyan Gao, Jun Wang, and Baode Sun. 2020. "First-Principles Study on the Cu/Fe Interface Properties of Ternary Cu-Fe-X Alloys" Materials 13, no. 14: 3112. https://doi.org/10.3390/ma13143112

APA Style

Wang, Y., Li, M., Gao, H., Wang, J., & Sun, B. (2020). First-Principles Study on the Cu/Fe Interface Properties of Ternary Cu-Fe-X Alloys. Materials, 13(14), 3112. https://doi.org/10.3390/ma13143112

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop