Michael R. Fellinger,Hyoungki Park,John W. Wilkins
Abstract:Large-scale simulations of plastic deformation and phase transformations in alloys require reliable classical interatomic potentials. We construct an embedded-atom method potential for niobium as the first step in alloy potential development. Optimization of the potential parameters to a well-converged set of density-functional theory (DFT) forces, energies, and stresses produces a reliable and transferable potential for molecular dynamics simulations. The potential accurately describes properties related to the fitting data, and also produces excellent results for quantities outside the fitting range. Structural and elastic properties, defect energetics, and thermal behavior compare well with DFT results and experimental data, e.g., DFT surface energies are reproduced with less than 4% error, generalized stacking-fault energies differ from DFT values by less than 15%, and the melting temperature is within 2% of the experimental value.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a reliable and transferable embedded - atom - method (EAM) potential function for large - scale simulations of niobium (Nb). Specifically, by optimizing the parameters of the EAM potential function to match the forces, energies, and stresses obtained from density - functional - theory (DFT) calculations, the researchers constructed a potential function that can accurately describe the structure, elastic properties, defect energies, and thermal behavior of niobium under different conditions. This potential function is intended for use in molecular - dynamics simulations, especially for simulations of processes involving long - range strain fields or long - wavelength fluctuations, which require a large number of atoms, and first - principles methods are difficult to handle such large - scale systems.
### Main objectives of the paper
1. **Develop EAM potential function**: Construct a niobium potential function based on the embedded - atom method that can accurately describe various physical properties of niobium, including structure, elasticity, defect energy, and thermal behavior.
2. **Optimize potential - function parameters**: By comparing with DFT calculation results, optimize the parameters of the EAM potential function to ensure its reliability and transferability under different conditions.
3. **Verify the accuracy of the potential function**: Verify the accuracy of the constructed EAM potential function by comparing it with DFT calculation results and experimental data.
### Specific problems
- **Structure and elastic properties**: Accurately describe the lattice parameters, elastic constants, and other structural and elastic properties of niobium.
- **Defect energy**: Accurately predict the defect energy in niobium, such as the vacancy - formation energy.
- **Thermal behavior**: Accurately describe the thermal behavior of niobium, such as the melting point.
- **Long - range strain fields and long - wavelength fluctuations**: Handle processes involving long - range strain fields or long - wavelength fluctuations, which require large - scale atomic simulations.
### Solutions
- **Force - matching method**: Use the force - matching method proposed by Ercolessi and Adams, and optimize the parameters of the EAM potential function to match the forces, energies, and stresses obtained from DFT calculations.
- **Database generation**: Generate a database containing niobium - atom configurations under different temperatures and strain conditions through DFT calculations.
- **Optimization algorithm**: Use the `potfit` program to optimize the parameters of the EAM potential function to ensure its accuracy and transferability under different conditions.
### Results
- **Structure and elastic properties**: The EAM potential function can accurately describe the lattice parameters, elastic constants, and other structural and elastic properties of niobium with an error within 4%.
- **Defect energy**: The EAM potential function can accurately predict the defect energy in niobium, such as the vacancy - formation energy.
- **Thermal behavior**: The EAM potential function can accurately describe the thermal behavior of niobium, such as the melting point, with an error within 2%.
- **Long - range strain fields and long - wavelength fluctuations**: The EAM potential function exhibits good performance when handling processes involving long - range strain fields or long - wavelength fluctuations.
In conclusion, this paper solves the problem of accurately describing the structure, elastic properties, defect energy, and thermal behavior of niobium under different conditions by developing and optimizing a reliable EAM potential function, providing a powerful tool for large - scale molecular - dynamics simulations.