First Principles Force Field for Metallic Tantalum

Alejandro Strachan,Tahir Cagin,Oguz Gulseren,Sonali Mukherjee,Ronald E. Cohen,William A. Goddard III
DOI: https://doi.org/10.48550/arXiv.cond-mat/0208027
2002-08-02
Abstract:We propose a general strategy to develop accurate Force Fields (FF) for metallic systems derived from ab initio quantum mechanical (QM) calculations; we illustrate this approach for tantalum. As input data to the FF we use the linearized augmented plane wave method (LAPW) with the generalized gradient approximation (GGA) to calculate: (i) the zero temperature equation of state (EOS) of Ta for bcc, fcc, and hcp crystal structures for pressures up to ~500 GPa. (ii) Elastic constants. (iii) We use a mixed-basis pseudopotential code to calculate volume relaxed vacancy formation energy also as a function of pressure. In developing the Ta FF we also use previous QM calculations of: (iv) the equation of state for the A15 structure. (v) the surface energy bcc (100). (vi) energetics for shear twinning of the bcc crystal. We find that withappropriate parameters an embedded atom model force field (denoted as qEAM FF) is able to reproduce all this QM data. Thus, the same FF describes with good accuracy the bcc, fcc, hcp and A15 phases of Ta for pressures from ~ -10 GPa to \~ 500 GPa, while also describing the vacancy, surface energy, and shear transformations. The ability of this single FF to describe such a range of systems with a variety of coordinations suggests that it would be accurate for describing defects such as dislocations, grain boundaries, etc. We illustrate the use of the qEAM FF with molecular dynamics to calculate such finite temperature properties as the melting curve up to 300 GPa; we obtain a zero pressure melting temperature of T_{melt}=3150 +/- 50 K in good agreement with experiment (3213-3287 K). We also report on the thermal expansion of Ta in a wide temperature range; our calculated thermal expansivity agrees well with experimental data.
Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an accurate force field (Force Field, FF) to describe the inter - atomic interactions of tantalum (Ta) under different temperature and pressure conditions. Specifically, the author proposes a general strategy based on quantum mechanics (QM) calculations to develop this force field and applies it specifically to tantalum. The main objectives include: 1. **Accurate description of multiple crystal structures**: The developed force field needs to be able to accurately describe the properties of tantalum in body - centered cubic (bcc), face - centered cubic (fcc), hexagonal close - packed (hcp) and A15 phases under different temperature and pressure conditions (from about - 10 GPa to about 500 GPa). 2. **Calculation of elastic constants**: Calculate the elastic constants by quantum mechanical methods and ensure that the force field can reproduce these data. 3. **Calculation of vacancy formation energy**: Calculate the vacancy formation energy after volume relaxation using the mixed - basis pseudopotential method and ensure that the force field can accurately describe this process. 4. **Calculation of surface energy**: Calculate the unrelaxed (100) surface energy and incorporate it into the training data of the force field. 5. **Thermodynamics of shear twinning**: Calculate the thermodynamic properties of bcc crystals under uniform shear, especially the energy change in the shear twinning mode. 6. **Calculation of melting curve**: Use molecular dynamics (MD) simulations to calculate the melting curve of tantalum under pressures up to 300 GPa. 7. **Calculation of coefficient of thermal expansion**: Calculate the coefficient of thermal expansion of tantalum in a wide temperature range by MD simulations and compare it with experimental data. Through these objectives, the author hopes to develop a force field that can be efficiently used in classical molecular dynamics simulations, so that the dynamic behavior of large systems (containing millions of atoms) on long - time scales can be studied, especially in terms of mechanical properties such as plasticity and failure.