arxiv: v1 [cond-mat.mes-hall] 7 Oct 2011

Transcription

1 Disorder regimes and equivalence of disorder types in artificial spin ice Zoe Budrikis, 1,2,3,a) Paolo Politi, 2,4 and R.L. Stamps 3 1) School of Physics, The University of Western Australia, 35 Stirling Hwy, arxiv: v1 [cond-mat.mes-hall] 7 Oct 2011 Crawley 6009, Australia 2) Istituto dei Sistemi Complessi CNR, Via Madonna del Piano 10, Sesto Fiorentino, Italy 3) SUPA School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, United Kingdom 4) INFN Sezione di Firenze, via G. Sansone 1, Sesto Fiorentino, Italy a) Electronic mail: zoe.budrikis@gmail.com 1

2 ABSTRACT The field-induced dynamics of artificial spin ice are determined in part by interactions between magnetic islands, and the switching characteristics of each island. Disorder in either of these affects the response to applied fields. Numerical simulations are used to show that disorder effects are determined primarily by the strength of disorder relative to inter-island interactions, rather than by the type of disorder. Weak and strong disorder regimes exist and can be defined in a quantitative way. 2

3 Artificial spin ice 1,2 is a two-dimensional frustrated system, consisting of arrays of singledomain sub-micron islands of magnetic material, which are coupled by stray magnetic fields. It was conceived as an experimental model of three dimensional spin ices 3 such as Dy 2 Ti 2 O 7 and Ho 2 Ti 2 O 7, with the advantage that its microstates can be imaged directly due to the relatively large island volume. A consequence of this large volume is that artificial spin ice is typically athermal, and dynamics must be induced by an external magnetic field, which acts to tilt the complex energy landscape and allow transitions between configurational states. Disorder in the system affects field-driven dynamics via two mechanisms. Disorder in island interactions (energetic disorder) alters the energies of configurational states. On the other hand, disorder in island magnetization switching characteristics (switching disorder) changes the barriers between states, altering the pathways the system takes through its space of configurational states 5. There are many possible sources of quenched disorder. In this paper we study four disorder types and show that their effects are similar. We also distinguish two regimes of disorder strength, namely the weak disorder regime and the strong disorder regime. Previous experimental studies of disorder have shown that the hysteresis of artificial spin ice reveals the importance of disorder 6. The switching characteristics of connected 7,8 and disconnected 9 kagome ice arrays have been used to estimate the strength of disorder in those systems. In simulation studies of an analog to artificial spin ice based on vortices in patterned superconductors, disorder is shown to nucleate grain boundaries between different ground state orderings 4. The four disorder types studied here are disorder in island positions, disorder in island orientations, a random perturbation to pairwise interactions and a random perturbation to island switching fields. Pairwise energy disorder is a simplified representation of energetic disorder in particular, it does not lead to correlations in island pair energies, unlike positional and orientational disorder. In all cases, disorder is assumed to have a Gaussian distribution. For example, positional disorder is implemented by perturbing the x and y coordinates of each island i by δ (i) x,y, where the δ come from a Gaussian distribution with mean 0 and standard deviation σ. We limit our considerations to square artificial spin ice 1. In square geometry, four islands meet at each vertex of the array, as shown in Fig. 1. The interactions of the four islands are 3

4 not all equivalent: four of the six pairwise interactions are stronger than the other two, due to the difference in distances and orientations between islands. As a result, there are four energetically-degenerate types of vertex configuration, and the ground state of the system is a tiling of Type 1 vertices. The population fraction of Type 1 vertices, n 1, is a good measure of the energy of the system. In the absence of disorder, the four vertex types have distinct, single-valued energies E 1 < E 2 < E 3 < E 4. However, when energetic disorder is present the energy of a vertex depends on the particular disorder realization. Considering many disorder realizations leads to a distribution of energies for each vertex type. The width of the distribution enables us to characterize the strength of disorder. We make this characterization in a model system in which the island magnetizations are Ising point dipoles, interacting so that the dipolar energy of island i is E (i) dip = h (i) dip M i = 1 4πµ 0 j i ( Mi Mj r 3 ij 3 ( M i r ij )( M ) i r ij ) rij 5 (1) with the island magnetization M and the nearest-neighbor distance both set to unity, so that the nearest-neighbor coupling has strength 1.5, in units of 1/(4πµ 0 ). In a perfect system, this gives E 1 = 6+ 2, E 2 = 2, E 3 = 0, E 4 = Shown in Fig. 2(a c) are energy bands for each vertex type, defined by the mean energy plus/minus one standard deviation, plotted against disorder strength σ for the three energetic disorder types. The bands become wider as the disorder strength increases, and eventually overlap occurs. We take the weak disorder regime to be that where there is no overlap, and the threshold of the strong disorder regime being the point where two bands first cross (in practice, the first crossing is always of the Type 2 and Type 3 bands, since they are closest initially). The transition occurs at σ(e pair ) for pairwise energy disorder, σ(φ) 0.3 for island orientational disorder and σ(r i ) 0.05 for island positional disorder. (Note that r i = r x,y refers to each component of an island s position.) The fourth type of disorder we discuss has a different nature, because it is disorder in island magnetization switching characteristics. We model the switching of the Ising moment of island i as occuring when h (i) tot ˆm i < h c, (2) 4

5 where h (i) tot is the sum of the external field h and the dipolar field h (i) dip and ˆm i is the direction of the island s magnetic moment. Unlike energetic disorder, switching disorder leaves the energies of configurational states unchanged and cannot be characterized in terms of vertex energies. However, an analogous measure can be constructed by considering the applied field required for two single spin flip processes. When all switching fields are equal, the applied field required to convert a Type 2 vertex into a Type 3 vertex is h c +1/ 2, while the field required for the reverse process is h c 1/ 2. These processes are shown as insets to Fig. 2(d). We take the boundary between weak and strong disorder to be the meeting point of the plus/minus one standard deviation bands for the fields required for the two processes. The standard deviation of the switching fields, σ(h c ), is equal to the standard deviation in the applied field required, so the crossing point can be easily determined to be σ = 1/ 2 0.7, as seen in Fig. 2(d). Note that this measure of disorder strength is relative to the inter-island coupling, not the mean switching field, since changing h c does not affect the disorder strength at which the bands meet. Thus, both energetic and switching disorder are measured relative to the same quantity and can be compared meaningfully. To see that the strength of disorder is more important than its origin, we compare in Figures 3 and 4 the final Type 1 population of a island square ice, attained under a rotating applied field for both a perfect system and systems subject to different types of weak disorder (Fig. 3) and strong disorder (Fig. 4). In the simulations, the field rotates with constant amplitude and angular step dθ = 0.01 radians for 10 cycles, long enough to obtain a steady state. At each field application, the system relaxes by flipping single spins according to criterion (2) until no further flips are possible. We have previously shown 10 that, in the absence of disorder, dynamics that create Type 1 vertices can only occur for a narrow range of field amplitudes h 2 and dynamics always start at array edges. Smaller fields are unable to flip spins, and larger fields force the magnetization to track the field. Between these limits, there are two field regimes. In the low field regime, h 11.25, dynamics proceed in a regular process of spin flips that invade from the array edges. In the high field regime, h 12, more dynamical processes are possible, but the final Type 1 vertex population is limited because dynamics must start at the array edges, and Type 1 vertices that are created near the array edges block this from happening. 5

6 The four different types of disorder studied all have similar effects on the n 1 vs h curves. In particular, switching disorder has the same effect as energetic disorder. This points to the fact that dynamics are determined through the spin flip criterion(2), which is affected equally by changes to the energy of configurations (through h dip ) and changes to the switching fields h (i) c. We expect similar results to hold for other switching criteria, such as Stoner-Wohlfarth switching, since the island magnetization reversal always depends on both coupling to other spins and the intrinsic properties of each island. On the other hand, even though different disorder types have similar effects, the difference between weak and strong disorder is striking. As expected, in the weak disorder regime disorder acts as a small perturbation on the perfect system. The two non-trivial field regimes arestill observed, and h ischangedlittlefromtheperfect system. Incontrast, inthestrong disorder regime the n 1 vs h curve is strongly altered from the perfect system. h increases to 4: at very low fields, spins with unfavourable h dip or small h (i) c can flip, whereas at very high fields the magnetization cannot completely track the rotating field due to pinning of spins with favourable h dip or large h (i) c. The other significant change is that n 1 is reduced in the low field regime and increased in the high field regime. This is a result of dynamics being allowed to start in the bulk. In the low field regime, this blocks the regular invasion process that would otherwise lead to large n 1. In the high field regime, this allows Type 1-creating dynamics to continue longer, increasing n 1 over the perfect case. In conclusion, different types of disorder in artificial spin ice have very similar effects on simulated dynamics and can all be compared to an energy scale set by the nearest-neighbor coupling. We show elswhere 11 that this allows disorder in an experimental system to be measured in terms of an effective switching field disorder. Z.B. and R.L.S. acknowledge the Australian Research Council for funding. Z.B. acknowledges funding from INFN and the Hackett Foundation. REFERENCES 1 R. F. Wang, et al., Nature 439, (2006). 2 Y. Qi, T. Brintlinger, and J. Cumings, Phys. Rev. B 77, (2008). 3 M. J. Harris, S. T. Bramwell, D. F. McMorrow, T. Zeiske, and K. W. Godfrey, Phys. Rev. Lett. 79, 2554 (1997). 6

7 4 A. Libál, C. J. Olson Reichhardt, and C. Reichhardt, Phys. Rev. Lett. 102, (2009). 5 Z. Budrikis, P. Politi and R. L. Stamps, arxiv: (unpublished). 6 K. K. Kohli, et al., arxiv: (unpublished). 7 P. Mellado, O. Petrova, Y. Shen, and O. Tchernyshyov, Phys. Rev. Lett. 105, (2010). 8 S. Ladak, et al., Nature Phys. 6, (2010). 9 E. Mengotti, et al., Nature Phys. 5, (2010). 10 Z. Budrikis, P. Politi, and R. L. Stamps, Phys. Rev. Lett. 105, (2010). 11 Z. Budrikis, J. P. Morgan, J. Akerman, A. Stein, R. L. Stamps, P. Politi, S. Langridge, C. H. Marrows, manuscript in preparation. 7

8 Type 1 Type 2 Type 3 Type 4 FIG. 1. The sixteen vertex configurations of square artificial spin ice can be classified into four types, based on energy. a E 4 b E 3 E 1 E 2 c d FIG. 2. The mean energy plus/minus one standard deviation of a (top-to-bottom) Type 4, 3, 2, 1 vertex, subject to (a) pairwise energy, (b) orientational and (c) positional disorder. The point wherethe Type 2 and Type 3 bandscross is taken to bethe transition between the weak and strong disorder regimes. The distributions are calculated over 200 independent disorder realisations. (d) The mean plus/minus one standard deviation of the external field required to convert a Type 2 to a Type 3 vertex (upper band) and the reverse (lower band). The insets show the two processes. 8

9 n1, final h Orientation Position Pairwise Switching No disorder FIG. 3. Final Type 1 population vs applied field h for a perfect system (black open circles) and four different types of weak disorder: island position disorder with standard deviation in each component of the positions of 0.01 (green filled triangles); island orientation disorder with standard deviation 0.06 (red filled circles); island pair energies perturbed with standard deviation (blue filled squares); and switching field disorder with standard deviation 0.14 (orange empty triangles). In all cases, the disorder strength σ = 0.2σ, where σ is the transition from weak to strong disorder. n1, final Orientation Position Pairwise Switching No disorder h FIG. 4. Final Type 1 population vs applied field h for a perfect system with open edges (black open circles) and four different types of strong disorder: island position disorder with standard deviation in each component of the positions of 0.06 (green filled triangles); island orientation disorder with standard deviation 0.36 (red filled circles); island pair energies perturbed with standard deviation 0.27 (blue filled squares); and switching field disorder with standard deviation 0.84 (orange empty triangles). In all cases, the disorder strength σ = 1.2σ, where σ is the transition from weak to strong disorder. 9