Abstract:Molecular dynamics simulations of the interface structure in binary AgCu eutectic were performed by using the realistic EAM potential. In simulations, we examined such quantities as the time dependence of the total energy in the process of equilibration, the probability distributions, the composition profiles for the components, and the component diffusivities within the interface zone. It is shown that the equilibrium in the solid state associated with the complete phase separation in binary eutectic is accompanied by formation of the steady disordered diffusion zone at the boundary of the crystalline components. At higher temperatures, closer to the eutectic point, the increase in the width of the steady diffusion zone is observed. The particle diffusivities grow therewith to the numbers typical for the liquid phase. Above the eutectic point, the steady zone does not form, instead, the complete contact melting in the system occurs. The results of simulations indicate that during the temperature increase the phenomenon of contact melting is preceded by the similar process spatially localized in the vicinity of the interface.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is related to the interface structure and its behavior of binary eutectic systems (especially the Ag - Cu eutectic) in the solid state and near the eutectic point. Specifically, the research aims to explore the following aspects:
1. **Interface Region Structure**: Through molecular dynamics simulations, study the structural characteristics of the solid - phase interface region in binary eutectic systems, especially the behavior near the eutectic temperature.
2. **Diffusion Region Formation**: Examine the stable disordered diffusion region formed during the complete phase separation in the solid - state equilibrium process, and study how its width changes with temperature.
3. **Contact Melting Phenomenon**: Explore the contact melting phenomenon that occurs above the eutectic point, as well as the process of local pre - melting near the eutectic point.
4. **Diffusion Coefficient**: Analyze the diffusion coefficients of components in the interface region, especially the diffusion characteristics of copper and silver atoms.
### Research Background
Binary eutectic systems are widely present in nature and industrial applications, such as sintering in powder metallurgy, synthesis of electronic and photonic materials, etc. Understanding the interface properties of these systems is crucial for grasping the structure formation process, especially phenomena near the eutectic temperature such as contact melting (CM). Traditionally, it is believed that a binary eutectic will completely decompose into ideal crystal phases under solid - state thermal equilibrium, and the interface width is zero. However, the observed diffusion regions and contact melting phenomena in experiments indicate that this is not the case in reality, and the actual interface width may be non - zero.
### Research Methods
To verify these phenomena and provide microscopic explanations, the author used the molecular dynamics (MD) simulation method, combined with the embedded - atom - method (EAM) potential function to model the Ag - Cu eutectic system. The main simulation parameters include:
- Temperature range: 600 K to 1050 K
- Lattice type: face - centered cubic (FCC)
- Lattice constants: \(a_{\text{Cu}} = 3.61 \, \text{\AA}\), \(a_{\text{Ag}} = 4.085 \, \text{\AA}\)
- Atomic masses: \(m_{\text{Cu}} = 63.5 \, u\), \(m_{\text{Ag}} = 107.9 \, u\)
### Main Findings
1. **Below the Eutectic Temperature**: When the temperature is below the theoretical eutectic point (\(T_{eut} = 935 \, \text{K}\)), the system reaches a stable state and forms a disordered diffusion region with a width of 5 - 20 Å. At this time, the diffusion coefficient increases significantly, approaching the liquid - phase level.
2. **Near the Eutectic Temperature**: As the temperature rises to near the eutectic point, the width of the diffusion region increases, the degree of atomic disorder intensifies, and the diffusion coefficient further increases.
3. **Above the Eutectic Temperature**: After exceeding the eutectic point, the width of the diffusion region continues to grow, eventually leading to complete contact melting, that is, the system transforms into the liquid phase.
### Conclusion
Through molecular dynamics simulations, the research shows that the contact between the ideal crystal phases of the binary Ag - Cu eutectic system is thermodynamically unstable in the solid state, and a stable disordered diffusion region will be formed at the interface. This phenomenon is similar to the pre - melting phenomenon in single - component systems, but in binary eutectics, it is manifested as contact melting. In addition, the diffusion coefficient in the diffusion region is significantly higher than that inside the crystal, approaching the liquid - phase level. These results provide a new perspective for understanding the interface behavior of eutectic systems and are helpful for explaining the simulation results based on the phase - field theory model.
### Formula Example
When calculating the three - dimensional diffusion coefficient of copper atoms, the Einstein relation was used:
\[ D = \frac{\langle \Delta r^2 \rangle}{6 \Delta t} \]
where \(\langle \Delta r^2 \rangle\) is the mean - square displacement (MSD), and \(\Delta t\) is the time interval.