The shape evolution of two-dimensional islands through periphery diffusion biased by an electromigration force is studied numerically using a continuum approach. We show that the introduction of crystal anisotropy in the mobility of edge atoms induces a rich variety of migration modes, which include oscillatory and irregular behavior. A phase diagram in the plane of anisotropy strength and island size is constructed. The oscillatory motion can be understood in terms of stable facets which develop on one side of the island and which the island then slides past. The facet orientations are determined analytically.PACS numbers: 66.30.Qa, The manipulation of nanostructures by macroscopic forces is likely to become a key ingredient in many nanotechnology applications. Understanding the influence of external fields on the shape evolution of nanoscale surface features is therefore of considerable importance. As a first step in this direction we analyze here the effects of an electric current on single-layer islands on a crystalline surface. The islands evolve under surface electromigration, the directed motion of adsorbed atoms due to the slight force transmitted by collisions with the conduction electrons in the sample [1].Electromigration along interfaces and grain boundaries is the most persistent and menacing reliability problem in integrated circuit technology [2,3]. Correspondingly, much work has been devoted to electromigration-induced void formation and breakdown in metallic conductor lines [3], and the capacity for quantitative numerical modeling has been demonstrated at least for simple void geometries [4,5]. A major obstacle to achieving predictive power in such studies, however, is the insufficient control over the complex internal structure of the polycrystalline samples. Hence an important motivation for investigating electromigration-induced effects on simple, well-controlled nanoscale morphologies, such as step patterns on vicinal surfaces [6] and single layer islands [7], is to bridge the gap between the microscopic mechanisms of electromigration and their consequences on technologically relevant length and time scales.Electromigration of islands has been modeled previously using Monte Carlo simulations [8] and continuum theory [9]. The continuum approach to island shape evolution, which treats the island edge as a smooth curve, has been successfully applied to a range of problems including the diffusion [10] and sintering [11] of islands, and the pinch-off of vacancy clusters [12]. Here we focus on the regime of periphery diffusion (PD), where the dominant kinetic process is the migration of atoms along the island boundary. The shape then follows a local evolution law, without coupling to the adatom concentration on the surrounding terrace.We extend the model of [9] by including crystal anisotropy in the adatom mobility. It was observed recently in the context of step flow growth [13] that crystalline anisotropy can change the behavior of step patterns in a qualitative way. In the present case, it lea...