We address the question of whether results obtained for small indenters scale to indenter sizes in the experimental range. The quasicontinuum method is used in order to extend the computational cell size to 2 2 1 m 3 , nominally containing of order 2:5 10 11 atoms, and in order to permit consideration of indenter radii in the range 70 ÿ 700 A. The dislocation structures for the large indenter are found to be less sharp and to extend over a larger region than for the small indenter. In addition, the large-indenter force-displacement curve differs from that corresponding to the small indenter in one important respect, namely, the absence of force drops during indentation, despite profuse dislocation activity. Based on these observations, we conclude that the indenter force is not a reliable indicator of the onset of dislocation activity and plastic deformation for indenter sizes in the experimental range. DOI: 10.1103/PhysRevLett.90.226102 PACS numbers: 02.70.-c, 31.15.-p, 46.15.-x The objective of this work is to ascertain the effect of indenter radius on the mechanics of nanoindentation of ductile fcc crystals at zero temperature (see, e.g., [1][2][3][4][5][6][7][8][9] for experimental background). For definiteness, we specifically focus on the indentation of Au (001) by spherical indenters of tip radii 70 and 700 A. One reason why the tip radius size is of concern is that straight molecular dynamics calculations have often relied on artificially small indenters in order to reduce the size of the computational cell [10 -20], and the scaling of the results to larger indenters is not straightforward. In particular, it is not immediately clear how dislocation activity depends onand how effective properties such as the force vs depthof-indentation relation scales with -indenter size.The technique we use in order to sidestep the size restrictions of straight molecular dynamics is the quasicontinuum method of Tadmor and co-workers [21][22][23]. Thus, by coarsening the level of spatial resolution of the computational mesh away from the indenter we are able to span realistic material samples of the order of 2 2 1 m 3 , which nominally contain 2:5 10 11 atoms. This helps to capture the elastic field of the indenter without introducing spurious or parasitic effects associated with periodicity or small sample sizes. By adapting the mesh size to the deformation field, the calculations provide full atomistic resolution over an appropriate region under the indenter. In this manner, the calculation is reduced to a coarse-grained system of size several orders of magnitude smaller than the original one, without appreciable loss of accuracy. In addition, this model reduction enables us to consider a wide range of indenters (70 to 700 A in the work presented here).The particular implementation of the quasicontinuum method used in the calculations has been described in [23]. In particular, following [21-23] we adopt as adaption indicatorwhere K denotes a simplex in the triangulation, II E K denotes the second invariant of the L...