# Melting artificial spin ice

## Abstract

Artificial spin ice arrays of micromagnetic islands are a means of engineering additional energy scales and frustration into magnetic materials. Here we demonstrate a magnetic phase transition in an artificial square spin ice and use the symmetry of the lattice to verify the presence of excitations far below the ordering temperature. We do this by measuring the temperature dependent magnetization in different principal directions and comparing with simulations of idealized statistical mechanical models. Our results confirm a dynamical pre-melting of the artificial spin ice structure at a temperature well below the intrinsic ordering temperature of the island material. We thus create a spin ice array that has real thermal dynamics of the artificial spins over an extended temperature range.

###### pacs:

75.10.Hk, 75.78.-nGeometric frustration is observed in many physical systems. The textbook example is the frustration of proton interactions in water ice, giving rise to proton disorder, as revealed in the pioneering experimental work of Giauque and Stout (1936) and the theoretical interpretation of Pauling (1935). Frustration in antiferromagnets analogous to the ice model was predicted theoretically by Anderson in 1956 Anderson (1956). Only much later was ice-type disorder observed in magnets, and most surprisingly, this was in ferromagnetic materials that were thus named ‘spin ice’ Harris *et al.* (1997); Ramirez *et al.* (1999); Bramwell and Gingras (2001). The spin ice phenomenon relies on dipole-dipole interactions and the concept was later generalized to include arrays of magnetic islands with a spin ice type geometry, the artificial spin ices Wang *et al.* (2006); Lammert *et al.* (2010); Ladak *et al.* (2010); Mengotti *et al.* (2010); Cumings (2011); Schumann *et al.* (2010); Morgan *et al.* (2011). Beyond the ice-type systems, the property of frustration is responsible for the occurrence of thermodynamically metastable phases in a variety of systems, including structural (amorphous materials) Senkov and Miracle (2001), magnetic (spin glasses) Kirkpatrick (1977) and polymeric systems Menon and Pandit (1995). In general, the absence of unique ground states Lammert *et al.* (2010) and the presence of kinetic constraints Ritort and Solich (2003) makes the study of excitations extremely important. These, in combination with the degree of frustration, are the driving forces for disorder in such systems Fecht (1992). An order-disorder transition has been proposed for colloid spin ice model systems with adjustable interactions Libál *et al.* (2006) and dynamical studies have already proven the importance of distortions Han *et al.* (2008).

The aim of the present work is to realize a thermally induced transition in an artificial spin ice array. In previous work on artificial magnetic spin ice, dynamics have been induced by changes in an applied magnetic field, rather than by thermal excitations Wang *et al.* (2007); Nisoli *et al.* (2007); Ke *et al.* (2008); Nisoli *et al.* (2010). The use of field-cycling protocols has enabled the realization of quasi-degenerate spin states reminiscent of real spin ice Wang *et al.* (2006); Cumings (2011), as well as ordered states that are in static equilibrium Schumann *et al.* (2010). All these results are obtained at temperatures far below the ordering temperature of the patterned structures. Recent work by Morgan et al. Morgan *et al.* (2011) explored for the first time ground states in artificial spin ice, obtained during the growth of the pattern.

Here we used -doped Pd(Fe) to enable the exploration of the thermally driven order-disorder transition of the artificial spin ice structure. The Curie temperature of -doped Pd(Fe) films can be adjusted by variation of the iron (Fe) layer thickness in the palladium (Pd) Pärnaste *et al.* (2007); Papaioannou *et al.* (2010). A 1.2 monolayer thick Fe layer was embedded between 10 monolayers of Pd, which resulted in a Curie temperature of K. The film was patterned by an electron beam lithography and Ar ion milling into an artificial square spin ice array, of 750250 nm islands with a period of 1000 nm, extending over a 1.51.5 mm area (Fig. 1). The total magnetic moment for one island was determined to be Am at 5 K.

Magnetic hysteresis loops were recorded using the magneto-optical Kerr effect (MOKE) at temperatures from K to K Papaioannou *et al.* (2010) (see inset in Fig. 2). The field was applied in two directions: (i) parallel to the [10] and (ii) parallel to the [11] directions of the patterned array, see Fig. 1. When the field is removed, different populations of magnetic states for the different directions are obtained. For example, in the absence of thermal excitations, such a field cycling effectively removes vertices of type I, due to symmetry reasons. The normalized remanent magnetization, , of the magnetic island array is compared to the continuous sheet material in Fig. 2. The remanent magnetizations follow a similar trend up to K, but at higher temperatures the magnetic island array crosses over to a regime in which the magnetization is strongly reduced, approaching zero at around K, while that of the continuous film falls to zero at K.

In order to discuss the disordering of the magnetic island array, we distinguish between the microscopic magnetic moments within the islands (micro-spins, ) and the island macro-spins, composed of all moments within an island (). The coupling between the micro-spins is due to the exchange interaction, with a coupling constant related to the of the magnetic material, while the coupling between the islands is magnetostatic. These interactions are long ranged but the characteristic energy scale can be defined through the interaction between macro-spins, , on neighboring islands. The most relevant interaction is between the nearest neighbours, coupling spins along rows or columns of the array. We define this coupling as , where is of order and .

Introducing the thermally averaged quantity implies that thermal fluctuations of the macro and micro degrees of freedom are partially decoupled, allowing the two sets to order at different temperatures. This assumption could break down for smaller units, in which case the two ordering processes could be renormalized into a single transition, but it should be valid for the islands studied here, consisting of micro-spins. Ordering of the macro-spins due to magnetostatic interactions should therefore occur at a temperature , as even in the limit , the term ensures that at . Thus we expect disordering of the macro-spins at and a range of temperatures , in which the islands are decoupled.

The macro-spin is confined to lie along the long axis of the island through the magnetic shape anisotropy , which minimizes the collective dipolar energy of an island for the moment in one of the two configurations with Ising symmetry. The energy barrier for reversing the direction of can be approximated by , with , and is temperature dependent through . Here is the permeability of free space, is the demagnetization factor of the island, and its volume. The active magnetic material is assumed to have a thickness of 1 nm Osborn (1945) when estimating . The macro-spin reversal time scale can be assumed to follow a Néel-Arrhenius law where is the Boltzmann constant and for inter-island magnetic fluctuations to occur, this barrier must be thermally accessible.

The temperature dependence of the effective moment of the islands is captured by using , with K, and , which matches the temperature variation of the remanent magnetization of the continuous film shown in Fig. 2 to a good enough approximation for our purposes here. This leads to , with K and at K, which ensures accessible thermal fluctuations for macro-spins below .

The possible vertex configurations for artificial square spin ice are summarized in Fig. 1(d). Two dimensional square ice does not have the same narrow band of low energy states as three dimensional spin ice, as the long range part of the dipolar interaction is less well screened than in the three dimensional counterpart Isakov *et al.* (2005). As a result, the band of Pauling states cannot be considered quasi-degenerate at any level of approximation. The lowest energy configurations are made up of type I vertices, leading to an ordered ground state, with staggered magnetic moments. The type I vertices, which carry no moment, are effectively removed during field cycling and are therefore irrelevant in the absence of thermal excitations. The lowest energy states which both satisfy the ice rules and carry a magnetic moment are the subset of type II vertices. Hence the simplest model, that captures the essential physics observed in this multi-scale problem is a sixteen-vertex model in which the type I states are raised in energy, leaving a four-vertex (type II) manifold of ground states.

Within this scenario the application of a saturation field along the [10] direction would result in a partially lifted degeneracy at low temperature, in favor of the ensemble of states shown in Fig. 1(b). Applying the field along the [11] direction, followed by a return to remanence would lead to the selection of the fully ordered state illustrated in Fig. 1(c). To capture the temperature dependence of the excitations we can therefore compare the ratio of the remanent and saturation magnetization, , for these two directions. The order parameter of the islands, , cancels out, leaving the ordering behavior of the macro-spins only. The resulting ratio should therefore be for the [10] direction and for the [11] direction. As seen in Fig. 3), the measured ratio falls on the prediction up to 125 K, illustrating collective behavior of the macro spins, with the single magnetic domain structure illustrated in Fig. 1(b, c). Remarkably, the four vertex behavior, in which vertex defects can be neglected, remains up to temperatures where thermal fluctuations have reduced the squared order parameter of the islands, , to less than one half of its zero temperature value. In the case of the [11] field, the data falls below the value at 125 K, while that for the [10] field stays on the line up to a slightly higher temperature. The magnetization, for the field applied in both directions, finally collapses as the ordering temperature of the material, K, is approached.

We performed Monte Carlo simulations on a system with macro spins, to model the experimental results. A sixteen-vertex ice model of Ising spins was used, arranged on a square lattice. To mimic the experimental conditions the energy of type I vertices was raised above those of type II and for simplicity we set them equal in energy to the type IV vertices, giving a single energy scale, , . This model is equivalent to a grill of interpenetrating Ising spin chains coupled uniquely through the first neighbours along the chains [see Fig. 1(d)]. The ground state is thus an degenerate manifold of type II states with ferromagnetic ordering along chains (rows) of length but with disorder among the chains. The parameter defines the effective coupling of macro-spins within a chain and its temperature dependence is obtained by scaling with . As expected, there is no phase transition to an ordered state at finite temperature and in zero field the system remains disordered to . However the model does have remarkable finite size effects related to the zero temperature critical point of the one dimensional Ising model Ising (1925); fsI (). This effect is illustrated in Fig. 4, which shows the temperature evolution of the magnetic moment for macro spins on a single row (or column) of spins, for different values of . In each case the data are compared with , calculated for a one dimensional Ising chain of macro-spins fsI (), with a coupling constant , which is exact to leading order in . The Monte Carlo results are excellently reproduced by the analytic expression for all values of , with small differences coming from corrections to scaling Hof (). The chain moment approaches saturation below a crossover temperature, which we can interpret as when modeling the array. It scales logarithmically to zero with system size, thus for any reasonable system size, there is a finite region of temperatures for which each chain of the vertex model is ordered.

This finite size chain ordering is at the origin of the observed pre-melting as it ensures that the field hysteresis protocol, leading to , will order the chains giving a non-zero total moment for the vertex model below this crossover temperature scale. For a field along the [10] direction, the chains parallel to the field will be ordered but those perpendicular will be free, while for the field along [11] all chains should be ordered, as in the experiment. This field induced ordering is confirmed in Fig. 3, where we compare calculated and experimental remanent to saturation magnetization ratio values. The numerical data was taken from simulations using two fitted parameters: K, corresponding to (see Fig. 4), and a remanent applied field T. The simulation and experimental data are in good agreement up to around K, above which the experimental data falls more rapidly towards zero. This could arise due to the development of magnetic fluctuations perpendicular to the Ising spin axis or thermal macro-spin flips becoming more favorable. These are not taken into account numerically.

The best fit value of compares favorably with calculated values: the direct interaction of the island magnetic poles (treated as point charges) yields K. In addition, this value of yields a crossover temperature, which is compatible with the experimental data. The general agreement between experiment and simulation and the closeness of the fitted parameters to expected values therefore suggests that the sixteen vertex model proposed here, is indeed an accurate description of the experimental data.

In both experiment and simulation, the crossover from the ordered spin ice regime to the pre-melting regime begins at a lower temperature for the [11] field direction than for the [10] direction. Pre-melting begins as vertex excitations disobeying the ice rules occur, breaking up the correlated chains of macro-spins. In such an excitation, a vertex of type II becomes one of type III through the flipping of the direction of a macro-spin. In the case of remanent magnetization along the [10] direction, there are vertex excitations that preserve the magnetic moment, for example type II type III , as illustrated in Fig. 1(d). However, for along [11] all excitations are moment reducing. This suggests that the magnetization ratio will remain near the upper bound in the first stages of pre-melting for ordering along the [10] direction while the ratio will be reduced from this bound as soon as defects appear, for ordering along the [11] direction. These defects carry accumulations of effective magnetic charge and in the case of three dimensional spin ice, have been successfully described using the language of emergent magnetic monopoles Castelnovo *et al.* (2008); Jaubert and Holdsworth (2009); Giblin *et al.* (2010). In this case, the term ‘monopole’ should be used with great caution, as long ranged interactions are less well screened for four-fold symmetry in a plane, than in spin ice. Here we have demonstrated the possibility to measure the order to disorder transition in artificial magnetic spin ice. The possibility of observing emergent monopoles is therefore conceivable, following the general approach which we describe in the design of spin ice arrays.

The authors acknowledge the support of the Swedish Research Council (VR), the Knut and Alice Wallenberg Foundation (KAW), the Swedish Foundation for International Cooperation in Research and Higher Education (STINT), the Icelandic Research Fund for Graduate Students and the Institut Universitaire de France (IUF). S.T.B. thanks the UK EPSRC for financial support. The patterning was performed at the Center for Functional Nanomaterials (CFN), Brookhaven National Laboratory, which is supported by the U.S. Department of Energy, Office of Basic Energy Sciences, under Contract No. DE-AC02-98CH10886. V.K., U.B.A. and B.H. thank Aaron Stein for support received during patterning at CFN. We thank Simon Banks and Michel Gingras for valuable discussions.

### References

- W. F. Giauque and J. W. Stout, Journal of the American Chemical Society, 58, 1144 (1936).
- L. Pauling, Journal of the American Chemical Society, 57, 2680 (1935).
- P. W. Anderson, Phys. Rev., 102, 1008 (1956).
- M. J. Harris et al., Phys. Rev. Lett., 79, 2554 (1997).
- A. P. Ramirez et al., Nature, 399, 333 (1999).
- S. T. Bramwell and M. J. P. Gingras, Science, 294, 1495 (2001).
- R. F. Wang et al., Nature, 439, 303 (2006).
- P. E. Lammert et al., Nature Physics, 6, 786 (2010).
- S. Ladak et al., Nature Physics, 6, 359 (2010).
- E. Mengotti et al., Nature Physics, 7, 68 (2010).
- J. Cumings, Nature Physics, 7, 7 (2011).
- A. Schumann et al., Appl. Phys. Lett., 97, 022509 (2010).
- J. P. Morgan et al., Nature Physics, 7, 75 (2011).
- O. N. Senkov and D. B. Miracle, Materials Research Bulletin, 36, 2183 (2001).
- S. Kirkpatrick, Phys. Rev. B, 16, 4630 (1977).
- G. I. Menon and R. Pandit, Phys. Rev. Lett., 75, 4638 (1995).
- F. Ritort and P. Solich, Adv. Phys., 52, 219 (2003).
- H. J. Fecht, Nature, 356, 133 (1992).
- A. Libál, C. Reichhardt, and C. J. O. Reichhardt, Phys. Rev. Lett., 97, 228302 (2006).
- Y. Han et al., Nature, 456, 898 (2008).
- R. F. Wang et al., J. Appl. Phys., 101, 09J104 (2007).
- C. Nisoli et al., Phys. Rev. Lett., 98, 217203 (2007).
- X. Ke et al., Phys. Rev. Lett., 101, 037205 (2008).
- C. Nisoli et al., Phys. Rev. Lett., 105, 047205 (2010).
- M. Pärnaste et al., Journal of Physics: Condensed Matter, 19, 246213 (2007).
- E. Th. Papaioannou et al., Journal of Physics: Condensed Matter, 22, 236004 (2010).
- J. A. Osborn, Phys. Rev., 67, 351 (1945).
- S. V. Isakov, R. Moessner, and S. L. Sondhi, Phys. Rev. Lett., 95, 217201 (2005).
- E. Ising, Z. Phys., 31, 253 (1925).
- The finite size Ising chain is often discussed in passing as a prelude to results in the thermodynamic limit. Two recent papers that focus on the model in its own right are Y. Li and B-G. Liu, Phys. Rev. B 73, 174418 (2006) and V. I. Belokon’ et al., Bulletin of the Russian Academy of Sciences: Physics 74, 1413 (2010). The analytical expression is: .
- J. D. Hofgartner, Physics-437B freport, University of Waterloo (2011).
- C. Castelnovo, R. Moessner, and S. L. Sondhi, Nature, 451, 42 (2008).
- L. D. C. Jaubert and P. C. W. Holdsworth, Nature Physics, 5, 258 (2009).
- S. R. Giblin et al., Nature Physics, 7, 252 (2010).