Quantitative decoding of interactions in tunable nanomagnet arrays using first order reversal curves.

To develop a full understanding of interactions in nanomagnet arrays is a persistent challenge, critically impacting their technological acceptance. This paper reports the experimental, numerical and analytical investigation of interactions in arrays of Co nanoellipses using the first-order reversal curve (FORC) technique. A mean-field analysis has revealed the physical mechanisms giving rise to all of the observed features: a shift of the non-interacting FORC-ridge at the low-HC end off the local coercivity HC axis; a stretch of the FORC-ridge at the high-HC end without shifting it off the HC axis; and a formation of a tilted edge connected to the ridge at the low-HC end. Changing from flat to Gaussian coercivity distribution produces a negative feature, bends the ridge, and broadens the edge. Finally, nearest neighbor interactions segment the FORC-ridge. These results demonstrate that the FORC approach provides a comprehensive framework to qualitatively and quantitatively decode interactions in nanomagnet arrays.

Nanomagnet arrays are basic building blocks1 for key technologies such as ultrahigh density magnetic recording media234, magnetic random access memory (MRAM)56, and logic devices789. Interactions within the arrays critically affect the functionalities of the nanomagnets as well as enable new device concepts. For instance, dipolar interactions may trigger an analog memory effect in nanowire arrays10, enable digital computation in magnetic quantum-dot cellular automata systems7811, lead to frustrations in artificial “spin ice”121314, or adversely affect thermal stability and switching field distribution in magnetic recording media or MRAM elements5151617. Probing and managing these interactions is often difficult because they are long-ranged, anisotropic, and configuration-dependent17.

The first-order reversal curve (FORC) method1819 has provided detailed characterization for a variety of magnetic2021222324252627 and other hysteretic systems2829. However, a coherent framework to interpret the features of the FORC diagrams and extract quantitative information is still lacking despite decades of effort by numerous groups. In this work, using the FORC method we have quantitatively investigated tunable interactions in model systems of single domain nanomagnet arrays. With mean-field level simulations, supplemented with a cluster extension, we have reproduced all features and trends of the experimental FORC diagrams quantitatively and identified their physical origins. Our approach decodes interactions in nanomagnet arrays, even disordered arrays, and also presents a pathway to evaluate (de)stabilizing interactions in other hysteretic systems.

Results

Rectangular arrays of polycrystalline Co ellipses were fabricated with varying center-to-center spacing by magnetron sputtering, in conjunction with e-beam lithography and lift-off. Details are presented in Methods. The ellipses have major/minor axes of 220/110 nm, with a structure of Ta(1 nm)/Co(9 nm)/Ta(1 nm), forming 50 × 50 μm2 arrays. In arrays A1/2/3 the minor-axis spacings of 150/200/250 nm are less than the major-axis spacing of 500 nm. Thus the mean dipolar interactions are demagnetizing, favoring anti-parallel alignment. In arrays B1/2/3 the minor-axis spacing of 500 nm exceeds the major-axis spacings of 250/300/350 nm. Therefore, the mean dipolar interactions are magnetizing, favoring parallel alignment. Scanning electron microscopy (SEM) and magnetic force microscopy (MFM) images, at remanence after DC demagnetization, of arrays A1 and B1 are shown in Figs. 1 and 2, respectively. MFM image contrast indicates the out-of-plane stray fields and confirms the ellipses' single domain state.

Figure 1

Demagnetizing Arrays: (a) SEM and (b) MFM image of the DC-demagnetized A1 array. Dashed ovals outline single ellipses, while the dashed box outlines an example of the checkerboard pattern. (c) Experimental and (d) simulated FORC distributions for the A2 array.

Figure 2

Magnetizing Arrays: (a) SEM and (b) MFM image of the DC-demagnetized B1 array. Dashed ovals outline single ellipses. (c) Experimental and (d) simulated FORC distributions for the B2 array.

FORC measurements were performed to obtain magnetization M(H, H) under different reversal field H and applied field H. The FORC-distribution is then extracted18, , where M is the saturation magnetization. In loose analogy to the Preisach model, the FORC-distribution in certain simple cases can be interpreted as the 2-dimensional distribution of elemental hysteresis loops with unit magnetization called hysterons on the (H, H) plane, or on the corresponding (H) plane, defined by local coercivity H = (H-H)/2 and bias/interaction field H = (H+H)/2.

Good agreement between measured and simulated FORC-distributions was obtained for all studied arrays [e.g., Figs. 1(c) and (d) for array A2]. All demagnetizing FORCs A1/2/3 exhibited a ridge with the high-H end on the H axis and the low-H end shifted in the +H direction30. Increasing interactions in A3→A2→A1 increased the low-H shift and the length of the ridge. In addition, an edge emerged from the low-H end towards negative H, highlighted by the arrows, forming a “wishbone” or boomerang structure24. A negative feature is observed at negative H values near the high-H end. Similarly, FORC distributions for the magnetizing arrays also exhibit a ridge with the high-H end on the H axis, but with the low-H end shifted towards −H [e.g., Figs. 2(c) and (d) for array B2]. Increasing interactions again increase the low-H shift, but reduce the length of the ridge30. A negative feature below the ridge is more prominent than that in the demagnetizing case.

Non-interacting case

A FORC gives a non-zero contribution to ρ(H, H) only if dM/dH along the FORC depends on H. We first show that in the non-interacting case ρ coincides with the coercivity distribution D(H), spread along the H axis. Indeed, particle P with coercivity H down-flips at H and up-flips at H. Therefore, on a FORC starting at H > −H, P starts out up-flipped and remains so, not contributing to dM/dH nor ρ.

In contrast, on FORC(H, H), P is the last to down-flip along the major loop, and has the highest coercivity among the down-flipped particles. Therefore, P is the last to up-flip as H increases past H on the same FORC, causing a dM/dH > 0 jump. Since this dM/dH jump is unmatched by the neighboring FORC(H, H), dM/dH exhibits a dependence on H, making d(dM/dH)/dH non-zero. dM/dH increases as Hdecreases, making a positive contribution to ρ at (H = H).

For all subsequent FORCs at H < −H, P starts out down-flipped but still up-flips at H = H. The dM/dH jumps on these FORCs are matched since they occur at the same field on each FORC(H, H < −H). Thus dM/dH is independent of H, and doesn't contribute to ρ. This reasoning highlights that only dM/dH jumps on individual FORCs that are unmatched by neighboring FORCs contribute to ρ. Each particle P contributes to ρ only once, at (H = H) or equivalently at (H). The contributions of all particles gives rise to a ridge along the H (H) axis, which reflects D(H).

Interacting case

Next, we introduce interactions between nanomagnets on the mean-field level by including an interaction field H, where α<0 for demagnetizing systems and α>0 for magnetizing ones21. Fig. 3(a) shows a sequence of FORCs for a demagnetizing system with a rectangular D(H), and a zoom-in view of the boxed region (right panel). The three FORC segment-pairs (1)/(2)/(3) show that the last dM/dH jump on each FORC(H, H = −H-αM(H) - caused by the last up-flipping particle P(H) - is unmatched by the neighboring FORC(H, H > −H-αM(H).

Figure 3

(a) Schematic illustration of family of FORCs for arrays with a flat coercivity distribution and mean-field demagnetizing interactions, with bold lines and numbers indicating unmatched dM/dH jumps. Calculated FORC distributions are shown in (b) illustrating the construction of the ridge and (c) the edge. (d) FORC distribution with the same interactions, but a Gaussian coercivity distribution; emergent negative feature is indicated by the dashed boundary. Similar panels are shown in (e–h) for the magnetizing case.

Fig. 3(b) shows that with interactions the unmatched dM/dH jumps still generate the ridge, but at shifted H values. Importantly, on the mean-field level all particles experience the same interaction field and thus the order of flips continues to be governed by the order of the coercivities: along the major loop the particles down-flip in the order of their coercivities, lowest (highest) coercivity particle P(H) first [P(H) last]. Starting at the low-H end, the lowest coercivity particle P(H) down-flips at H, as no other particles have flipped yet: M(H. Increasing H along FORC(H, H = H), P(H) up-flips at H, causing a positive jump dM/dH > 0 as shown by the lower FORC(H, H = H)-segment of pair (1). This jump, caused by P(H), is absent on the upper FORC(H, H)-segment and is thus unmatched, contributing to ρ at (H = H), or similarly at (H), defining the low-H end. These flipping fields are shifted from their non-interacting values, as shown by the arrow set (1) in Fig. 3(b). Since P(H) defines the low-H end of the FORC-ridge, one concludes that interactions shift the low-H end of the FORC-ridge to the +H direction by −αM (recall α < 0), but leave its H coordinate un-shifted at H.

FORC-segment pairs (2)/(3) in Fig. 3(a) illustrate the up-flip of higher coercivity particles P(H). P(H)s create unmatched dM/dH jumps on FORCs where they were the last to down-flip at H = −H [vertical arrows in Fig. 3(b)] and also are the last to up-flip at H = H [horizontal arrows in Fig. 3(b)], since M(H when P(H) up-flips.

The high-H end of the FORC-ridge is defined by P(H) that is the last to down-flip when the rest of the system is already negatively saturated (M = −M), thus H. P(H) up-flips along the FORC(H, H = H) only after the rest of the system is positively saturated: H. Accordingly, the unmatched dM/dH jumps caused by P(H) contributes to ρ only at (H = H), or similarly at (H, H), defining the high-H end. As observed before, the high-H end of the FORC-ridge remains on the H axis24, but stretched along the H axis by −αM (α<0).

Note that interactions shift the FORC ridge feature in H uniformly [Fig. 3(b)], i.e., the resultant projection of FORC distribution onto the H-axis is only displaced from its intrinsic values, but not distorted. This H-projection therefore mirrors the non-interacting case, where H, reflecting the intrinsic coercivity distribution, simply displaced by −αM. Thus the intrinsic coercivity distribution can be - without distortion from interactions – directly identified from the FORC distribution.

Figs. 1(c,d) and 3(b,c) show that besides the ridge, ρ exhibit an edge as well with interactions27. As discussed earlier, in the absence of interactions, the dM/dH jumps along a FORC(H, H = −H) are matched by the jumps on the subsequent FORC(H, H < −H)s, not contributing to ρ. The arrows of Fig. 3(a) show that the interactions destroy this matching specifically at the low-H by shifting the first up-flip field H of each FORC, caused by P(H), by −αM(H. These shifts make the dM/dH jumps misaligned, see FORC-segment-pairs (4) and (5) (above and below), thus contributing to ρ at (H = H, H). These unmatched jumps give rise to the edge in Fig. 3(c). The end-points of the edge are (H = H) and (H = H), or alternatively (H) and (H). Accordingly, the tilt and asymmetry of the edge provide a direct measure of the width of coercivity distribution D(H) = H. In the extreme case of nearly identical nanomagnets with D(H) = H≈0, the edge is vertical in the H plane, at H and H within ± αM, as observed experimentally in Ni nanowire arrays10.

In short, on the mean-field level the FORC-distribution of a system with demagnetizing interactions exhibits an edge and a ridge, shifted by the unmatched first and last dM/dH jumps along each FORC. The FORC-distribution vanishes between them for the considered flat coercivity distribution, because jumps between the first and the last jumps along each FORC are matched by jumps on the neighboring FORCs. Here, the matched jumps are not caused by the same particles, as in the non-interacting system, but rather by different particles whose up-flipping fields were shifted into alignment by the interactions. Still, the jumps are matched because the values of the aligned jumps are the same on neighboring FORCs for a flat distribution D(H). Visibly, the flat coercivity distribution on the mean-field already reproduces most features of the measured FORC-distribution.

To improve our model we introduce a more realistic Gaussian D(H) to elucidate the origin of the negative features, which represent a clear distinction between FORC and a literal Preisach interpretation. A Gaussian breaks the matching of jumps as now shown by examining the set of particles {P(H)} around the center of the coercivity distribution. A P(H) particle is the last to down-flip at H, where M(H, and the last to up-flip on the FORC(H, H), contributing to ρ at (H = H. On FORC(H, H)s P(H) is no longer the last to up-flip. On subsequent FORCs, up-flip jumps from different particles get shifted into alignment with this jump. For the flat distribution, the number of particles shifted into alignment is steady, making the match complete and thus zero contribution to ρ. However, for the Gaussian distribution, the particles shifted into alignment come from the decreasing slope of the Gaussian, leading to dM/dH jumps with a decreasing magnitude, providing only a partial-match and generating a negative contribution to ρ. Fig. 3(d) shows the FORC-distribution for a Gaussian model that indeed develops a negative region specifically tracking the decreasing slope of the Gaussian, highlighted by the dashed line. Analogous arguments show that a Gaussian D(H) also bends the ridge and broadens the edge.

Nearest-neighbor correlations

The last unaccounted feature of the measured ρ is the segmenting of the FORC-ridge into separate low-H and high-H ends, with different amplitudes. To explain this we include the nearest-neighbor interaction fields H as the first terms of a systematic cluster-expansion.

Decreasing the field from positive saturation, the weakest coercivity particles down-flip first. For demagnetizing interactions, H of these down-flipped particles stabilize their nearest-neighbors in their up state. Therefore, for a sufficiently narrow D(H, the magnetization decreases towards zero by developing a checkerboard pattern [Fig. 1(b)]. The checkerboard naturally forms defects where the sequence of increasing coercivities selects third-nearest neighbor particles to down-flip. Still, the dominant reversal mechanism for nearly half of the particles is the checkerboard formation: down-flipping with all-up neighbors. Accordingly, the (nearly) half of the FORC ridge with H gets shifted along the +H axis by H, where H is the interaction field for the all-neighbors-up configuration. Once the checkerboard pattern is formed, the rest of the particles flip with neighbors in various intermediate configurations. Therefore, the H half of the ridge is broken into several pieces, shifted by varying H fields. Consequently, the nearest neighbor interactions manifest as segmenting of the ridge30.

To reiterate, the demagnetizing interactions (α < 0) (1) shift the non-interacting FORC-ridge at the low-H end to the +H direction by −αM; (2) stretch the non-interacting FORC-ridge at the high-H end along H by −αM without shifting it off the H axis; and (3) form a tilted edge connected to the ridge at the low-H end. Changing from flat to Gaussian D(H distribution (4) produces a negative feature, bends the ridge, and broadens the edge. Finally, (5) nearest neighbor interactions segment the FORC-ridge.

Magnetizing interactions

Adapting the above arguments for magnetizing interactions (α > 0): (1) the low-H end is shifted in the −H direction, and (2) the high-H end is compressed without shifting it off the H axis [Fig. 3(f)]. (3) Regarding the edge, the first up-flip along each branch is shifted by interactions in the opposite direction as the demagnetizing case [Fig. 3(e) right panel]. Therefore, the first dM/dH jumps are unmatched, decreasing in magnitude with more negative H, thus negatively incrementing the FORC, forming a negative edge [Fig. 3(g)]. Changing from flat to Gaussian D(H distribution (4) the negative edge gets pressed towards the positive ridge, and the FORC-ridge becomes curved [Fig. 3(h)]. The inclusion of nearest neighbor terms leads to (5) an avalanche reversal, collapsing the FORC-ridge to a single-value31.

Quantifying Interaction Fields

Finally, we demonstrate the quantitative predicting power of the above considerations. The lowest H, which is shifted in H by αM, is extracted from the FORC ridge by selecting an H such that 10% of the particles have H, and averaging ρ over the H→H range. The averaged (dM/dH)' are shown in insets of Figs. 4(a) and 4(b). The H shift is determined by linearly extrapolating (dM/dH)' at the high |H| end to zero. The interaction field is calculated by a finite element method (using the NIST OOMMF code) for the nearest and next nearest neighbors, and treating the remainder of the array as point dipoles. The experimental H shifts and the calculated interaction fields agree remarkably well (Fig. 4), confirming the validity of the mean-field description of the FORC-distribution and its quantitative predictive power, making the FORC technique a tool to extract numerical values of interaction fields. This is particularly important for disordered arrays where calculations of interactions are not easily achievable.

Figure 4

Calculated (open symbols) and experimentally determined (solid symbols) interaction field for (a) demagnetizing arrays A1/2/3 and (b) magnetizing arrays B1/2/3.

Averaged FORC-distribution utilizing the H are shown in insets for the (a) A2 and (b) B2 array, where the linear extrapolation is illustrated by the dashed line and open circle.

Discussion

In this work, systems of interacting nanomagnets were examined experimentally, numerically, and analytically, using the FORC technique. A mean-field analysis based on the concept of unmatched jumps accounted for all experimentally observed features of the FORC diagram, including its shifted ridge-and-edge structure and negative features. The tilting, shifting, and stretching of these structures were identified as tools to extract quantitative information about the system, demonstrating the predictive power of the FORC technique. Construction of the FORC distribution through unmatched jumps, and recognizing the (de)magnetizing interactions as a particular case of (de)stabilizing interactions, presents an approach which can be used to evaluate any hysteretic system with the FORC technique.

Methods

Arrays of Co ellipses were fabricated by DC magnetron sputtering, in a vacuum chamber with a base pressure of 1 × 10−8 Torr and Ar sputtering pressure of 2 × 10−3 Torr, on Si (100) substrates, in conjunction with electron beam lithography and lift-off techniques. Magnetic hysteresis loops were measured at room temperature using the magneto-optical Kerr effect (MOKE) magnetometer with a 632 nm HeNe laser having a 30 μm spot-size, capturing the reversal behavior of ~5,000 ellipses32. The magnetic field was applied parallel to the major axis of the ellipses. Each measurement was averaged over ~103 cycles at a rate of 11 Hz. The arrays were coated with a 60 nm ZnS layer to improve the signal-to-noise ratio33. FORC measurements were performed as follows22: from positive saturation the magnetic field is swept to a reversal field H, where the magnetization M(H, H) is measured under increasing applied field H back to saturation, tracing out a FORC. The process is repeated for decreasing reversal field H30.

Ellipses were modeled as dipoles oriented parallel to their major axes. The inter-dipole spacing and magnetic moment per dipole in the 100 × 100 array were representative of the experimental system. Each dipole i was assigned an intrinsic coercivity H with a distribution experimentally determined from the sample having the weakest interactions, A3. The H dipolar interaction fields at dipole i were calculated on the mean-field level as αM(H), where α was calibrated such that αM equals the analytically calculated H at saturation. This mean-field formulation was extended by the first term of a cluster expansion, representing the nearest neighbor dipole interaction H explicitly: H = αM(H)+H. At each field step (ΔH = 1Oe) the total field H was compared to H, down-flipping occurred when H+Hi < −H and up-flipping occurred when H+Hi > H, until all dipoles became stable.

Author Contributions

D.A.G. obtained the experimental and simulation results, and wrote the first draft of the paper. G.T.Z. and M.W. participated in the simulation design. R.K.D., A.G., N.E., J.L.V. and K.L. participated in the experimental design, fabrication and characterization. K.L. and G.T.Z. designed and coordinated the whole project. All authors contributed to analysis, discussion and revision of the paper.