arxiv: v3 [cond-mat.mtrl-sci] 31 May 2018

Transcription

1 Magnetic order and energy-scale hierarchy in artificial spin ice Henry Stopfel, Erik Östman, Ioan-Augustin Chioar, David Greving, Unnar B. Arnalds, Thomas P. A. Hase, Aaron Stein, 4 Björgvin Hjörvarsson, and Vassilios Kapaklis Department of Physics and Astronomy, Uppsala University, Box 56, SE-750, Uppsala, Sweden Department of Physics, University of Warwick, Coventry, United Kingdom Science Institute, University of Iceland, Reykjavik, Iceland 4 Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 97, USA (Dated: June, 08) arxiv: v [cond-mat.mtrl-sci] May 08 In order to explain and predict the properties of many physical systems, it is essential to understand the interplay of different energy-scales. Here we present investigations of the magnetic order in thermalised artificial spin ice structures, with different activation energies of the interacting Ising-like elements. We image the thermally equilibrated magnetic states of the nanostructures using synchrotron-based magnetic microscopy. By comparing results obtained from structures with one or two different activation energies, we demonstrate a clear impact on the resulting magnetic order. The differences are obtained by the analysis of the magnetic spin structure factors, in which the role of the activation energies is manifested by distinct short-range order. This demonstrates that artificial spin systems can serve as model systems, allowing the definition of energy-scales by geometrical design and providing the backdrop for understanding their interplay. I. INTRODUCTION Initially introduced as a playground for the experimental investigation of magnetic frustration effects, the area of artificial spin ices has evolved into a vibrant field of research. Using designed structures, collective order and dynamics can be controlled and directly visualized by using, for example, nano-characterization techniques. This approach has motivated experimental investigations of celebrated classical spin models 5, some of which display non-conventional and exotic magnetic phases 6 0. Furthermore, given the large freedom in design, new topologies that promote non-conventional emergent order can be manufactured and investigated in thermallyactive artificial spin ice structures, allowing the study of their complex energetic manifolds in a superparamagnetic regime 7. A variety of new artificial spin ice systems have been proposed by Morrison et al. 8, with particular attention given to the so-called Shakti lattice, a variation of the square ice geometry. In this lattice, 5% of the square ice lattice elements have been removed and certain pairs of islands merged as depicted in Fig.. Theoretical and numerical studies have predicted an emergent magnetic order in the Shakti geometry 9 which retrieves the features of the six-vertex model 4 on a larger length-scale. This behavior has been demonstrated experimentally by Gilbert et al. 0 for the Shakti lattice (SH), as well as for a modified Shakti lattice (msh). Given the mesoscopic nature of the structures, different activation energies are expected for the long (T B ) and short islands (T B ), resulting in two blocking temperatures of the elements (T B < T B ). The presence of two activation energies is therefore expected to influence the eventual magnetic correlations and order in the lattice, since the larger elements could act as a source of quenched disorder. Such disorder would hinder the development of any medium- or long-range magnetic order of the smaller islands, as these become arrested at lower temperatures (T B ). Whilst the vertex statistics for both the Shakti lattices have been experimentally investigated for a range of inter-island coupling strengths 0, the thermally driven kinetics in the modified Shakti lattice were only recently reported 7. Consequently, the impact of the hierarchy of the activation energies on the emergent order has not been addressed to date. In this work, we undertake this task for the particular case of the Shakti and the modified Shakti geometries. Using thermallyactive elements, we characterize and compare their final magnetic configurations using identical experimental protocols, highlighting the impact of different activation energies on the obtained magnetic correlations and order. II. EXPERIMENTAL DETAILS A. Sample preparation The samples were produced by a post-growth patterning process, applied to a thin film of δ-doped Palladium (Iron) as described in Pärnaste et al.. A thick bottom layer of Pd (40 nm) is followed by.0 mono layers of Fe, defining the nominal Curie-temperature, T C = 400 K, as well as the thermally active temperature range for the δ-doped Pd(Fe). The post-patterning was carried out at the Center for Functional Nanomaterials (CFN), Brookhaven National Laboratory in Upton, New York. A positive high resolution e-beam resist (ZEP50A) was employed to create a Chromium mask of the nano-structures using e-beam lithography. The magnetic structures were formed by Argon milling. All investigated structures were produced on the same substrate from the same layer, ensuring identical intrinsic material

2 FIG.. A schematic representation of the Shakti (SH) and modified-shakti (msh) lattices along with magnetic states representative of each temperature regime, generally defined by the Curie (T c) and the blocking temperatures (T B/B ). Elements in their paramagnetic state are gray, superparamagnetic elements are light-blue and the frozen elements are represented with a magnetic north (blue) and south (red) pole, reflecting their final magnetic orientation. The gray-scale images to the left are representative PEEM-XMCD snapshots obtained at 65 K, the lowest temperature reached during the cooling. properties as well as the same thermal history for all the investigated structures. Each array had a spatial dimension of µm. The lengths of the large and small islands, were chosen to be 050 nm and 450 nm, respectively, while the width was kept at 50 nm. The lattice spacing between two parallel neighboring short islands was set to 600 nm. B. Thermal protocol Similar to previous studies 5,6, the thermal protocol involves a gradual cooling from a superparamagnetic state of the elements, towards a completely arrested state at the lowest temperatures. This process is schematically represented in Fig. for both the Shakti and the modified Shakti lattices. Notice the two distinctive stages of the superparamagnetic regime for the Shakti lattice because to the presence of elements with different intrinsic activation energies. When the temperature is lowered below the Curie temperature (T c ), the elements become magnetic and can be considered as mesospins. The stray field of these islands is resulting in a temperature-dependent magnetostatic coupling. The coupling between the elements has a direct impact on the involved activation energies, biasing the otherwise equally-favoured two magnetic states of the islands. On the lattice scale, this effectively results in a dynamic distribution of blocking temperatures around the intrinsic energy barriers, a distribution whose width depends on the strength of the inter-island interactions. C. Interaction strength and their impact There are three special cases to be considered: the strong, intermediate and weak coupling limits. In the strong coupling limit, the intrinsic activation energy of the mesospins is much smaller than the inter-island interactions. Consequently, the emergent order will not depend on the differences in the activation energies, but significant effects on the dynamic response cannot be ruled out. In the intermediate coupling limit, the strength of the inter-island couplings is comparable to the intrinsic activation energy of the mesospins. This will give rise to a complex interplay, in which the relation between the activation and interaction energies will dictate both the order and dynamics of the magnetic phases. The individual activation energies depend on the dynamically developing magnetic configurations, for both long- and short-mesospins. In other words, these couplings can yield a strong overlap in the two blocking temperature distributions of the short and long islands. Therefore, the configurations of the long- and the short-mesospins will exhibit different degrees of order, all depending on the relation of the activation and interaction energies. Conversely in the weak coupling limit, the intrinsic activation energies of the long-mesospins are much larger than the inter-island interactions. Around the intrinsic blocking temperature of the long islands (T B ), the available thermal energy is higher than the energy difference given by any magnetization reversal of the short islands. Therefore the lattice as a whole will be in a disordered state when the long-mesospins freeze (T B )

3 and consequently the magnetic orientation of the long elements will be randomly distributed, with no short- or long-range correlations. However, at the blocking temperature for the short-mesospins (T B ) the inter-island coupling strength can be comparable to their intrinsic activation energies, potentially resulting in a short-range order below T B. Thus, having two distinct activation energies is expected to have a significant impact on the emergent order, for both the intermediate and the weak coupling limits. By the same token, the presence of two or more activation energies is expected to have a marginal impact in the strong coupling limit. Here, we focus on the weak coupling limit and, as previously stated, we compare the obtained magnetic order in SH and msh lattices, ascertaining the impact of the hierarchical relation of inter-element interactions and the activation energies involved. D. Determination of the magnetization direction The magnetic states were determined as a function of temperature using Photo Emission Electron Microscopy (PEEM) capturing the X-ray Magnetic Circular Dichroism (XMCD) contrast. The experimental studies were performed at the.0. PEEM beamline at the Advanced Light Source, in Berkeley, California, USA. For the SH and msh lattices the islands were oriented 45 with respect to the incoming X-ray beam, to enable identification of the magnetization direction of all elements in the lattices. For these lattices multiple, partially overlapping, PEEM-XMCD images were acquired and stitched together to form one large image containing over 4000 mesospins for each of the investigated lattices, see Supplementary Materials Fig.. The images used for determining the magnetic states were taken at 65 K, far below the Curie-temperature and the blocking temperatures of the short- and long-mesospins. The sampling time (t s = 00 s) required for obtaining a PEEM-XMCD image defines a time window linked to the thermal stability of the mesospins: At the lowest temperatures, all magnetic states are stable (frozen) and the lattice can be imaged for practically arbitrary long times, yielding the same result. However, if the reversal times are much smaller than t s, no magnetic contrast is obtained. This allowed us to ask at which temperature the spins have reversed their magnetization direction within a given time window, thereby providing an estimation on the involved energy scales. III. RESULTS A. Activation energies of the mesospins Magnetic interactions can affect the inferred activation energies in a profound way. We therefore fabricated two additional lattices consisting of long- and short-mesospins, 0 nm and 450 nm respectively, in which the distance was large enough to allow us to ignore any effects arising from the interactions between the elements (see Supplementary Materials). The islands were oriented parallel with respect to the incoming X-ray beam, to determine the magnetization direction of the mesospins. Following the line of argumentation above, we determined that 50% of the short islands reverse their magnetization direction within 00 s at K, while a temperature of 49 K is required for obtaining the same condition for the long islands. The determination of the activation energies is model dependent. Using an expression for the probability of a magnetic element not having altered its magnetization direction after a time period t 6, s, as described in the Supplementary Materials, yields a ratio for the activation energies of 0.49(4). This approach provides support for well separated energy-scales, which will be assumed to be the case in the remainder of this communication. B. Correlations of the long-mesospins in the Shakti lattice The correlations between the long-mesospins were determined from the PEEM-XMCD images recorded at 65 K, after cooling the samples from room temperature with a rate of K/minute. The results from analysis along different directions in the SH lattice are illustrated in Fig.. As seen in the figure, no correlations between long-mesospins are observed. The cooling of the sample below the blocking temperature of the shortmesospins does not result in any observed order for the long-mesospins. Because of this, the effective interaction mediated by the short-mesospins must be much smaller than the activation energy of the long-mesospins. Correlation number 0 String Diagonal 4 Neighbor String Diagonal 4 Neighbor Fully ferromagnetic Random Fully antiferromagnetic FIG.. Long-mesospin correlations in the Shakti lattice. The network-averaged magnetic correlations for longmesospins spatially oriented vertically/horizontally are shown in the left/right graph. In blue/red the correlations of long-mesospins which are oriented in a string and in lightblue/orange the correlations of the long-mesospins diagonally oriented are represented, as indicated in the inset. All correlations show that already the first neighbor is randomly oriented and therefore the magnetic orientations of the long-mesospins are not correlated. Error bars represent one standard deviation.

4 4 The msh lattice has an additional degree of freedom as compared to the SH lattice, arising from the two-fold coordinated vertex having two possible states: Type-I and Type-II. The vertex statistics obtained for the three- and fourfold coordinated vertices are presented in Fig. (d) and (e), respectively. Comparing the abundance of the different states with published data from Gilbert et al.0, we can again confirm that our results are representative for the weak coupling regime, as e.g. the previous works on ASI structures,6. Notice that these vertex statistics do not display a significant difference between the common vertex states of the SH and msh lattices, similar to previous reports0. This implies that the final vertex populations are independent of the presence of the long islands in the Shakti geometry and the related energy hierarchy. However, the vertex statistics yield only a partial description of the final magnetic state and are insensitive to magnetic order established over spatial extensions larger than the vertex size. D. FIG.. Vertex type definition and statistics for SH and msh lattices. SEM image of the SH (a) and msh (b) lattice. The different vertex types are indicated by their respective colors. The real space base vectors for the Shakti lattice are defined along the x- and y-direction and each have a length of four times the lattice spacing (d = 600 nm). (c), Definition of the vertex states for all different vertex types. While both lattices share the vertex type with a coordination number of four (green), the msh lattice contains vertices with coordination number three (red) and two (orange). In the SH lattice the three-fold coordinated vertices (blue) include a long-mesospin. The four- (d) and three-fold (e) coordinated vertex statistics for SH (dark-blue) and msh (dark-red) lattices show no significant differences. Error bars represent one standard deviation. C. Beyond the vertex-level We use the calculated spin structure factor (SSF) to determine correlations beyond nearest neighbours7,8,0,. This invokes the use of a Fourier transformation of pairwise correlations of the magnetic states in the following form: I(q) = N X iq (ri rj ) Si Sj e N () (i,j=) where I(q) is the intensity for a given scattering vector q in reciprocal space, N denotes the number of spins over which the SSF is calculated, S i the spin component perpendicular to q and ri is the position of the corresponding spin in the real-space lattice. The SH and the msh lattices have different lattice points, ri, for the Magnetic order on the vertex-level Having established the randomness of the magnetic states of the long-mesospins in the Shakti lattice, we now turn our attention to the overall magnetic order established within the SH and msh lattices. For square ice related geometries, a commonly-employed tool for characterizing the obtained magnetic state is the determination of populations of the various types of vertex configurations. The different vertex types for both the SH and the msh lattices are illustrated in Fig. (a)-(c). In the nomenclature used here, a roman numeral denotes the vertex type along with a subscript representing the number of mesospins forming the vertex. As an example, I4 represents four elements forming a Type-I vertex. FIG. 4. Spin structure factor maps computed for the fourfold coordinated vertex sub-lattices of the msh (left) and the SH (right) lattices. The real space lattice vectors as well as the sub-set of four-fold coordinated vertices (red islands) are indicated in the insets.

5 5 (a) a a (b) q y -[r.l.u.] (d) q y -[r.l.u.] q x -[r.l.u.] q x -[r.l.u.] (c) q y -[r.l.u.] (e) q y -[r.l.u.] (f) q x -[r.l.u.] q x -[r.l.u.] q,-q -[r.l.u.] x y ---SH msh I I 4 I 9 I q,-q -[r.l.u.] x y FIG. 5. Spatially-limited spin structure factor as indicator for short-range order. (a), Illustration of the mesospin sub-set used in the SSF calculations. The real space lattice vectors (a, a ) which are used for the reciprocal lattice units (r.l.u.) are also shown. The real space input window size for the averaged SSF I (red SSF (b/d) ) and I 9 (blue SSF (c/e) ) are highlighted. b-e, The averaged SSF maps for SH (b/c) and msh (d/e) lattices, with reduced real space input size according to equation (). (f), Diagonal line-profiles for the SSF of different input sizes for the SH (dark-blue) and msh (dark-red) lattices. The modulation of the intensity in the SSF maps at (q x, q y) = (, ) [r.l.u.] is clearly distinguishable for the case of the SH lattice and an indication of a different order for these length-scales. long- and two short-mesospins, the differences in the SSF of the SH and msh data is dominated by this structural difference. To facilitate a proper comparison between the two cases, only the mesospin sub-set found in both geometries, namely the four-fold coordinated vertices, has been taken into account for the SSF computation. Therefore we calculate the SSF according to equation (), in the same manner as is explained in more detail by Östman et al.0, but only using the mesospins of the vertices with a coordination number of four. Since this sub-lattice of mesospins is the same in both the SH and msh lattices, differences in the SSF represent different ordering among these sub-sets of spins. Characteristic SSF maps for the SH and msh networks are presented in Fig. 4. Although neither of the two maps presents specific features, such as pinch points or Bragg peaks, their diffuse signals systematically exhibit clear differences. This is indicative of differences in the global spin order, which have not been captured by the statistics determined on the vertex-level (presented in Fig. (d) and (e) ). The striking feature is the presence of a hollow region (e.g. at the (±,±) (r.l.u.) points in reciprocal space) in the SH lattice SSF map, with the same region having a more uniform intensity distribution in the structure factor map of the msh. To understand the origin of this effect and identify the length-scales that are at play, we have computed magnetic structure factors for networks of gradually-increasing size and network-averaged these individual signals for each chosen network-size. Since the sample is large enough to assume ergodic conditions, this size reduction of the sampling window and the corresponding network-averaging, provides statistical information about the magnetic state on the length-scale of the sampling window rather than the high spatial frequencies of the speckle pattern (Fig. 4) obtained from the specific global microstate. Hence, we define I L = (4L) M M 4L m= (i,j=) S i S j e iq (ri rj) () as an incoherent 4 spatial average over sampling windows with L four-fold coordinated vertices. Each sampling window therefore consists of 4L spins (taking into account only the short-mesospins). The SSF I L is the average of all M sampling windows present in the measured global microstate. Fig. 5(b) and 5(d) show a comparison of I, i.e. the averaged SSF for a single four-fold coordinated vertex (four mesospins) for the SH and msh lattices. Both are virtually indistinguishable, which is not surprising since I is merely a different representation of the vertex statistics shown in Fig. (d). In other words, as the vertex statistics showed no significant differences, the SSF maps in Fig. 5(b) and Fig. 5(d) should also be similar, traceable also in a comparison of the top-most diagonal line-profiles in Fig. 5(f). Expanding the SSF input window to sixteen mesospins (I 4, four four-fold coordinated vertices), a clear difference between the SH and msh appears, with a more pronounced dip in the intensity at the (,) reciprocal space position for the SH lattice case (Fig. 5(f) ). Expanding the real space input window to thirty-six mesospins (I 9, nine four-fold coordinated vertices), the differences in the SSF maps between the SH (Fig. 5(c) ) and msh (Fig. 5(e) ) become even more pronounced (see also the line-profiles in Fig. 5(f) ). The line-profiles in Fig. 5(f) are diagonal cuts through the respective SSF maps I, I 4, I 9 and I 6

6 6 from reciprocal space position (q x, q y ) = (0, 0) [r.l.u.] to (q x, q y ) = (4, 4) [r.l.u.]. The results from these spatially-limited SSF maps in Fig. 5 are qualitatively similar to the features observed for the SSF of the full microstate in Fig. 4 and capture differences in the magnetic short-range order, exceeding the extension of just one vertex. (a) A A A E. Short-range magnetic order 5 4 We attribute the origin of the features observed in the SSF maps to the long elements in the SH lattice and the difference in activation energies between the short and long islands. To elaborate more on this premise, we have extracted the abundance of various spin arrangements that connect four-fold coordinated vertices across a Type-I vertex in the msh lattice or a long-mesospin in the SH lattice. Here, we focus on three of the different configurations of interest for such spin arrangements, as presented in Fig. 6(a). In the three vertex coupling classes, A, A and A, a total of nine different ferromagnetic spin arrangements are investigated. A detailed analysis of all nine spin arrangements can be found in the Supplementary Materials Fig.. In Fig 6(b) we report their abundances for both lattices, along with the estimated values given by a random distribution, in the absence of any interaction between the considered mesospins. A and A are averages of four different spin arrangements (C -C 4 and C 5 -C 8, see Supplementary Materials Fig. ), which have the same characteristic. In the case of A, one mesospin (mesospin or in schematic 6(a) ) is ferromagnetically aligned to the long-mesospin (or the two ferromagnetically aligned mesospins of the two-fold coordinated vertex) and one mesospin on the other side (mesospin or 4 in Fig. 6(a) ). For A one mesospin on one side is ferromagnetically aligned with two mesospins on the other side of the long-mesospin (or the two ferromagnetically aligned mesospins of the twofold coordinated vertex). A represents the lowest energy configuration for these five (six in the msh) mesospins and arranges all mesospins (mesospin,, and 4 in schematic Fig. 6(a) ) ferromagnetically with respect to the long-mesospin (or the two ferromagnetically aligned mesospins of the two-fold coordinated vertex). IV. DISCUSSION The spin arrangements in Fig. 6 highlight the development of a stronger correlation between four-fold coordinated vertices across a long-mesospin in the SH lattice than in the msh, in which case the coupling is mediated via a two-fold coordinated vertex. In the Shakti ground state manifold as described by Chern et al. 9, the fourfold coordinated vertices are always in their lowest energy state, Type-I 4 (see definition Fig. c). This state (b) Abundance A A A Ferromagnetic spin arrangement SH msh random FIG. 6. Spin arrangements across the long-mesospins / twofold coordinated vertices. (a), Schematic representation of all mesospins contained within two neighboring three-fold coordinated vertices (gray numbered islands). We classify the arrangements of all these mesospins in three classes, denoted as A, A and A. The abundance values of A and A are averages of the four possible spin arrangements in these classes. One representative spin arrangement for each of the three classes is depicted in (a). The coupled four-fold coordinated vertices are indicated by the light-green islands. The number of coupled four-fold coordinated vertices across the Type-I vertex or a long-mesospin is two for A, three for A and all four for A. (b), The abundance of these three classes is plotted for the SH (dark-blue) and the msh (dark-red) lattices, the abundance for each spin arrangement in a random mesospin state is represented with a gray inside bar. The abundance revealing the distinct differences in the arrangement of the mesospins around the Type-I vertex and is strongly indicating a defined tiling for the four-fold coordinated vertices around the long-mesospins and the build-up of a short-range order. Error bars represent one standard deviation. has a two-fold degeneracy, which we denote by Type- I 4 (A) and Type-I 4 (B) (see Supplementary Materials Fig. 4(b) ). A transition from Type-I 4 (A) to a Type-I 4 (B) is given by a magnetization reversal of all four mesospin. Two Type-I 4 vertices which are ferromagnetically cou-

7 7 pled via a Type-I or a long-mesospin will always be of opposite states, Type-I 4 (A) and Type-I 4 (B). Returning to the SSF maps in Fig. 4 and Fig. 5, the distinct shape in the SSF can be explained by the enhanced effective coupling of Type-I 4 (A) to Type-I 4 (B) vertices via the longmesospins in the SH lattice. In other words, although the long-mesospins freeze in a random configuration, they are more effective in transmitting magnetic correlations from one end of a Shakti plaquette to the other, than the twofold vertex of the msh lattice. While this is a rather straightforward result for the case of one plaquette, it is a priori expected that, on the lattice scale, the prior establishment of a random distribution of long-mesospin states would yield an even more randomized configuration of short islands, a fact that the msh lattice might circumvent given the absence of these constraints. To further discuss our results we would like to compare them with the ones previously reported by Gilbert et al. 0. For the strong coupling regime, they observed an almost perfect ground-state ordering of the four-fold coordinated vertices for the msh lattice, while the SH lattice has only about 80% of these vertices in the lowest energy state, Type-I 4. These differences, in the light of our results, can be attributed to the distinct energyscales. As the long-mesospins freeze into a pre-defined disordered sub-lattice, the Shakti ground state manifold is restricted and shaped by the long-mesospin sub-lattice. When coupling four Type-I 4 vertices via a long-mesospin, 6 different Type-I 4 -tilings are possible. The Shakti ground state, as defined by Chern et al. 9, can be obtained with only out of the total 6 tilings. The msh lattice, with its mono-sized short-mesospins, can access the entire ensemble of tilings. However, the pre-defined disordered long-mesospin arrangement effectively reduces the possible states in the Shakti ground state manifold (see Supplementary Materials Fig. 4). The frozen longmesospin restricts the Type-I 4 -tiling amongst four fourfold coordinated vertices to 7 out of 6 possible arrangements. As the system grows in size, the Shakti ground state manifold is further reduced by each added frozen long-mesospin (for further information see Supplementary Materials). This is a plausible cause for the higher order found in the msh lattice for the strong-coupling regime reported by Gilbert et al. 0 compared to the SH lattice. However, as the inter-island coupling strength is reduced, the long-mesospins, although a source of randomness, provide a better locally-correlated state than the two-fold coordinated vertices of the msh lattice. V. CONCLUSIONS Using thermally-active nano-patterned arrays of weakly-coupled magnetic mesospins, we have highlighted the impact of distinct length- and energy-scales on the development of magnetic order in artificial spin ice. While our study is restricted to the Shakti geometry and the weak-coupling regime, it captures the importance of the energy hierarchy in mesoscopic architectures and its role for the obtained magnetic order. Given the vast opportunities in geometrical design and the available choices of magnetic materials, the energy hierarchy can be engineered, and, in combination with suitable thermal or even demagnetization protocols, used as a tool for tailoring the magnetic order in mesoscopic magnetic structures. The hierarchy of the energy-scales, as well as the interplay between them, lie behind the complex emergent behavior and properties seen frequently across a range of natural systems 5. Artificial spin systems provide, therefore, a unique opportunity to study emergence, by offering a designer fabrication route and parametrization of the relevant variables. This opportunity potentially allows for a better understanding of problems in physics requiring multiple scale-analysis 6. VI. ACKNOWLEDGMENTS The authors would like to thank C. Nisoli and G.W. Chern for valuable discussions. The authors acknowledge support from the Knut and Alice Wallenberg Foundation, the Swedish Research Council and the Swedish Foundation for International Cooperation in Research and Higher Education. The patterning was performed at the Center for Functional Nanomaterials, Brookhaven National Laboratory, supported by the U.S. Department of Energy, Office of Basic Energy Sciences, under Contract No. DE-SC This research used resources of the Advanced Light Source, which is a DOE Office of Science User Facility under contract No. DE- AC0-05CH. U.B.A. acknowledges funding from the Icelandic Research Fund grants Nr. 458 and 548. R. F. Wang, C. Nisoli, R. S. Freitas, J. Li, W. McConville, B. J. Cooley, M. S. Lund, N. Samarth, C. Leighton, V. H. Crespi, and P. Schiffer, Artificial spin ice in a geometrically frustrated lattice of nanoscale ferromagnetic islands, Nature 49, 0 06 (006). C. Nisoli, R. Moessner, and P. Schiffer, Colloquium: Artificial spin ice: Designing and imaging magnetic frustration, Rev. Mod. Phys. 85, (0). I. Syôzi, Statistics of Kagomé Lattice, Progress of Theoretical Physics 6, (95). 4 E. H. Lieb, Exact Solution of the F Model of An Antiferroelectric, Phys. Rev. Lett. 8, (967). 5 E. H. Lieb, Residual Entropy of Square Ice, Phys. Rev. 6, 6 7 (967). 6 S. Zhang, I. Gilbert, C. Nisoli, G.-W. Chern, M. J. Erickson, L. O Brien, C. Leighton, P. E. Lammert, V. H. Crespi,

8 8 and P. Schiffer, Crystallites of magnetic charges in artificial spin ice, Nature 500, (0). 7 B. Canals, I.-A. Chioar, V.-D. Nguyen, M. Hehn, D. Lacour, F. Montaigne, A. Locatelli, T. O. Menteş, B. S. Burgos, and N. Roug le, Fragmentation of magnetism in artificial kagome dipolar spin ice, Nature Communications 7, 446 (06). 8 Y. Perrin, B. Canals, and N. Roug le, Extensive degeneracy, Coulomb phase and magnetic monopoles in artificial square ice, Nature 540, 40 4 (06). 9 C. Nisoli, V. Kapaklis, and P. Schiffer, Deliberate exotic magnetism via frustration and topology, Nature Physics, 00 0 (07). 0 E. Östman, H. Stopfel, I.-A. Chioar, U. B. Arnalds, A. Stein, V. Kapaklis, and B. Hjörvarsson, Interaction modifiers in artificial spin ices, Nature Physics 4, (08). J. P. Morgan, A. Stein, S. Langridge, and C. H. Marrows, Thermal ground-state ordering and elementary excitations in artificial magnetic square ice, Nature Physics 7, (0). U. B. Arnalds, A. Farhan, R. V. Chopdekar, V. Kapaklis, A. Balan, E. T. Papaioannou, M. Ahlberg, F. Nolting, L. J. Heyderman, and B. Hjörvarsson, Thermalized ground state of artificial kagome spin ice building blocks, Applied Physics Letters 0, 404 (0). A. Farhan, P. M. Derlet, A. Kleibert, A. Balan, R. V. Chopdekar, M. Wyss, L. Anghinolfi, F. Nolting, and L. J. Heyderman, Exploring hyper-cubic energy landscapes in thermally active finite artificial spin-ice systems, Nature Physics 9, 75 8 (0). 4 A. Farhan, P. M. Derlet, A. Kleibert, A. Balan, R. V. Chopdekar, M. Wyss, J. Perron, A. Scholl, F. Nolting, and L. J. Heyderman, Direct Observation of Thermal Relaxation in Artificial Spin Ice, Physical Review Letters, (0). 5 V. Kapaklis, U. B. Arnalds, A. Harman-Clarke, E. T. Papaioannou, M. Karimipour, P. Korelis, A. Taroni, P. C. W. Holdsworth, S. T. Bramwell, and B. Hjörvarsson, Melting artificial spin ice, New Journal of Physics 4, (0). 6 V. Kapaklis, U. B. Arnalds, A. Farhan, R. V. Chopdekar, A. Balan, A. Scholl, L. J. Heyderman, and B. Hjörvarsson, Thermal fluctuations in artificial spin ice, Nature Nanotechnology 9, (04). 7 Y. Lao, F. Caravelli, M. Sheikh, J. Sklenar, D. Gardeazabal, J. D. Watts, A. M. Albrecht, A. Scholl, K. Dahmen, C. Nisoli, and P. Schiffer, Classical topological order in the kinetics of artificial spin ice, Nature Physics (08), in press. 8 M. J. Morrison, T. R. Nelson, and C. Nisoli, Unhappy vertices in artificial spin ice: new degeneracies from vertex frustration, New Journal of Physics 5, (0). 9 G.-W. Chern, M. J. Morrison, and C. Nisoli, Degeneracy and Criticality from Emergent Frustration in Artificial Spin Ice, Physical Review Letters, 770 (0). 0 I. Gilbert, G.-W. Chern, S. Zhang, L. O Brien, B. Fore, C. Nisoli, and P. Schiffer, Emergent ice rule and magnetic charge screening from vertex frustration in artificial spin ice, Nature Physics 0, (04). M. Pärnaste, M. Marcellini, E. Holmström, N. Bock, J. Fransson, O. Eriksson, and B. Hjörvarsson, Dimensionality crossover in the induced magnetization of Pd layers, Journal of Physics: Condensed Matter 9, 46 (007). W. Wernsdorfer, E. Orozco, K. Hasselbach, A. Benoit, B. Barbara, N. Demoncy, A. Loiseau, H. Pascard, and D. Mailly, Experimental Evidence of the Néel-Brown Model of Magnetization Reversal, Physical Review Letters 78, (997). A. Farhan, C. F. Petersen, S. Dhuey, L. Anghinolfi, Q. H. Qin, M. Saccone, S. Velten, C. Wuth, S. Gliga, P. Mellado, M. J. Alava, A. Scholl, and S. van Dijken, Nanoscale control of competing interactions and geometrical frustration in a dipolar trident lattice, Nature Communications 8, 995 (07). 4 With the term incoherent, we mean that we do not account for the phase relationship between the averaged sampling windows. The correlations represented in the calculated SSF maps, are only up to lateral extensions of the respective sampling window size and are insensitive to correlations between or across multiple sampling windows. 5 P. W. Anderson, More Is Different, Science 77, 9 96 (97). 6 K. G. Wilson, Problems in Physics with many Scales of Length, Scientific American 4, (979). VII. AUTHOR CONTRIBUTIONS H.S. D.G. and A.S. fabricated the samples. H.S., E.Ö., U.B.A., and V.K. performed the PEEM-XMCD experiments. H.S., E.Ö., I.A.C., D.G, T.P.A.H., B.H. and V.K. analyzed the data and contributed to theory development. H.S., I.A.C., B.H. and V.K. wrote the manuscript. All authors discussed the results and commented on the manuscript.

9 Magnetic order and energy-scale hierarchy in artificial spin ice Henry Stopfel,, Erik Östman, Ioan-Augustin Chioar, David Greving, Unnar B. Arnalds, Thomas P. A. Hase, Aaron Stein, 4 Björgvin Hjörvarsson, and Vassilios Kapaklis arxiv: v [cond-mat.mtrl-sci] May 08 Department of Physics and Astronomy, Uppsala University, Box 56, SE-750, Uppsala, Sweden Department of Physics, University of Warwick, Coventry, United Kingdom Science Institute, University of Iceland, Reykjavik, Iceland 4 Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 97, USA (Dated: June, 08)

10 I. ACTIVATION ENERGIES To determine the activation energies or the ratio between the activation energies for shortand long-mesospins, we fabricated arrays of short- and long-mesospins with an edge to edge distance of 000 nm. For these islands this distance is large enough to consider the islands to be non-interacting. The fractions of identifiable islands for short (red) and long islands (blue) can be fitted with a double exponential of the following form: P (T ) = exp A B exp T A= 00 s ts = τ0 τ0 () B= E kb () () Equation () represents the probability of a magnetic element not having altered its magnetization direction after a time period ts, T is the temperature and A, B are fitting parameters. The measurement time for the PEEM-XMCD imaging was ts = 00 s. In equation (), τ0 is the inverse attempt frequency,, while in equation (), E is the activation energy and kb is Supplementary Figure. Temperature dependent measurements of the thermal dynamics of uncoupled mesospins. (a), Fractions of identifiable magnetic islands during a PEEM-XMCD measurement of 00 s and their corresponding fits. Both island types have a width of 50 nm, while the short islands have a length of 450 nm and the long of 0 nm. Examples of the PEEM-XMCD images for the short (b) and long islands (c) are presented for the measurement at 5 K.

11 Boltzmann s constant. The ratio of the activation energies for short (Eshort /kb = 560(40) K) and long islands (Elong /kb = 40(40) K) can be determined from the fits shown in Suppl. Fig. (a), and is 0.49(4). II. REAL SPACE IMAGES Real space images of the SH and msh lattices that are created by stitching together individual images (each recorded for 60 s) recorded in the limited field of view of the microscope. The composite images shown in Fig. (a) for the SH and Fig. (b) for the msh lattices represent the raw data from which the vertex statistics and spin structure factors were calculated. The individual islands can be seen as either black/white/gray depending on whether the islands are aligned parallel/anti-parallel/reversing with respect to the x-ray propagation direction. The particular composite images shown in Fig. were recorded at 65 K. Supplementary Figure. PEEM-XMCD measurements of the Shakti lattices. Composite PEEMXMCD images of the SH (a) and msh (b) lattices.

12 III. SPIN ARRANGEMENTS All nine possible spin arrangements C C 9, for the short- and long-mesospins (depending on the lattice) of two neighboring three-fold coordinated vertices are depicted in Fig. (a). The abundance of the spin arrangements is presented in Fig. (b). These spin arrangements can be categorized into three groups marked with orange, yellow and green. Group one (orange A ), are arrangements C C 4 connecting one mesospin, and the respective vertex, with another mesospin on the other side of the long-mesospin (two short-mesospins). The second group (yellow A ), are arrangements C 5 C 8 in which one mesospin on one side is in ferromagnetic alignment with the long-mesospin (two short-mesospins) and two mesospins on the other side. Finally, the last arrangement (green), C 9, is the lowest energy state for this sub-set of mesospins. Roughly 5% of all investigated long-mesospin connections of the SH have at least two vertices being correlated to each other (C C 4 ), but only % of them are correlated in the msh via the two short-islands. Similar differences in the correlations can also be found (a) C C C 4 C 5 C 7 C 8 (b) SH msh C C 6 C 9 normalizedbcorrelation C C C C 4 C 5 C 6 C 7 C 8 C 9 CorrelationbType Supplementary Figure. Spin arrangement arround the long-mesospin (Type-I ). (a), Illustration of the five (six in msh) islands which are considered in the spin arrangements. Adjacent islands are indicated by light red and in the different representation of the spin arragements C -C 9 with light green if correlated over the long-mesospin (Type-I ). (b), Abundance of the different spin arrangements, for the SH (dark-blue squares) and msh (dark-red circles) lattices. 4

13 for the other spin arrangements (C 5 C 9 ). This can be attributed to the high abundance (4%) of Type-II excitations in vertices with coordination number two, which are therefore unable to mimic the SH order. IV. TYPE-I 4 -TILING (a) (b) (c) (d) Type-I 4 (A) Type-I 4 (B) Supplementary Figure 4. Shakti ground state manifold in plaquettes. (a), Illustration of the Shakti ground state manifold in a plaquette as introduced by Chern et al. 4. In the Shakti ground state all four-fold coordinated vertices are in their lowest energy state, Type-I 4, and 50% of the three-fold coordinated vertices are in their lowest energy state, Type-I, and 50% are in the excited state, Type-II. These excitations are represented with a circle in the plaquette illustration by Chern et al.. (b), Equivalent representation of the Shakti ground state manifold for a plaquette, focusing on the Type-I 4 tiling at the four-fold coordinated vertices. All Type-I 4 tilings for a pre-defined long-mesospin direction are illustrated for one (c) and two plaquettes (d). The ground state manifold of the Shakti lattice was described by Chern et al. 4 and is best explained on the size of a plaquette. In Fig. 4(a) the extent of a plaquette is illustrated, consisting of 9 islands in the SH and 0 islands in the msh case. A Shakti plaquette has 5

14 one long island (or a two-fold coordinated vertex in the msh) in its center and 8 short islands at its borders, sharing 4 three-fold coordinated vertices and 4 four-fold coordinated vertices with the neighbouring plaquettes. Fig. 4(a) also shows the representation of the Shakti ground state manifold of such a plaquette as introduced by Chern et al. 4. Solid lines illustrate the lowest energy states, but the circles at the three-fold coordinated vertices mark the position of Type-II excitations. In the Shakti ground state a plaquette inhabits two of such Type-II excitations, at any of their 4 three-fold coordinated vertices. At the same time all the four-fold coordinated vertices are in their lowest energy state, Type-I 4. Therefore we represent the tiling of the two degenerated Type-I 4 (A) and Type-I 4 (B) vertices on the bases of a plaquette in Fig 4(b). The excitations in the Shakti lattice are located at the boundaries between Type-I 4 (A) and Type-I 4 (B) vertices. We discussed in the main text, that the pre-defined long-mesospins in the SH lattice will restrict the possible Type-I 4 -tilings. This effect is illustrated in Fig. 4(c)-(d) for one and two plaquettes. All 7 (out of 6) Type-I 4 -tilings accessible, by simultaneously remaining in the Shakti ground state manifold, regarding the given longmesospin direction are illustrated in Fig. 4(c). For two pre-defined long-mesospins (two neighbouring plaquettes Fig. 4(d) ) just (out of 64) Type-I 4 -tilings remain in the Shakti ground state manifold. henry.stopfel@gmail.com W. Wernsdorfer, E. Orozco, K. Hasselbach, A. Benoit, B. Barbara, N. Demoncy, A. Loiseau, H. Pascard, and D. Mailly, Experimental Evidence of the Néel-Brown Model of Magnetization Reversal, Physical Review Letters 78, (997). V. Kapaklis, U. B. Arnalds, A. Farhan, R. V. Chopdekar, A. Balan, A. Scholl, L. J. Heyderman, and B. Hjörvarsson, Thermal fluctuations in artificial spin ice, Nature Nanotechnology 9, (04). S. Bedanta and W. Kleemann, Supermagnetism, Journal of Physics D - Applied Physics 4 (009). 4 G.-W. Chern, M. J. Morrison, and C. Nisoli, Degeneracy and Criticality from Emergent Frustration in Artificial Spin Ice, Physical Review Letters, 770 (0). 6