Atomistic simulation of tantalum nanoindentation: Effects of indenter diameter, penetration velocity, and interatomic potentials on defect mechanisms and evolution

Nanoindentation simulations are a helpful complement to experiments. There is a dearth of nanoindentation simulations for bcc metals, partly due to the lack of computationally efficient and reliable interatomic potentials at large strains. We carry out indentation simulations for bcc tantalum using three different interatomic potentials and present the defect mechanisms responsible for the creation and expansion of the plastic deformation zone: twins are initially formed, giving rise to shear loop expansion and the formation of sequential prismatic loops. The calculated elastic constants as function of pressure as well as stacking fault energy surfaces explain the significant differences found in the defect structures generated for the three potentials investigated in this study. The simulations enable the quantification of total dislocation length and twinning fraction. The indenter velocity is varied and, as expected, the penetration depth for the first pop-in (defect emission) event shows a strain rate sensitivity m in the range of 0.037-0.055. The effect of indenter diameter on the first pop-in is discussed. A new intrinsic length-scale model is presented based on the profile of the residual indentation and geometrically necessary dislocation theory. (C) 2014 Elsevier B.V. All rights reserved.