Magnetic dipole induced guided vortex motion

Magnetic dipole induced guided vortex motion

N. Verellen INPAC-Institute for Nanoscale Physics and Chemistry, Nanoscale Superconductivity and Magnetism Pulsed Fields Group, K. U. Leuven Celestijnenlaan 200 D, B-3001 Leuven, Belgium.    A. V. Silhanek INPAC-Institute for Nanoscale Physics and Chemistry, Nanoscale Superconductivity and Magnetism Pulsed Fields Group, K. U. Leuven Celestijnenlaan 200 D, B-3001 Leuven, Belgium.    V. Metlushko Department of Electrical and Computer Engineering, University of Illinois, Chicago, IL 60607.    W. Gillijns INPAC-Institute for Nanoscale Physics and Chemistry, Nanoscale Superconductivity and Magnetism Pulsed Fields Group, K. U. Leuven Celestijnenlaan 200 D, B-3001 Leuven, Belgium.    F. Gozzini Department of Electrical and Computer Engineering, University of Illinois, Chicago, IL 60607.    B. Ilic Cornell Nanofabrication Facility, School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853    V. V. Moshchalkov INPAC-Institute for Nanoscale Physics and Chemistry, Nanoscale Superconductivity and Magnetism Pulsed Fields Group, K. U. Leuven Celestijnenlaan 200 D, B-3001 Leuven, Belgium.
July 27, 2019

We present evidence of magnetically controlled guided vortex motion in a hybrid superconductor/ferromagnet nanosystem consisting of an Al film on top of a square array of permalloy square rings. When the rings are magnetized with an in-plane external field , an array of point-like dipoles with moments antiparallel to , is formed. The resulting magnetic template generates a strongly anisotropic pinning potential landscape for vortices in the superconducting layer. Transport measurements show that this anisotropy is able to confine the flux motion along the high symmetry axes of the square lattice of dipoles. This guided vortex motion can be either re-routed by 90 degrees by simply changing the dipole orientation or even strongly suppressed by inducing a flux-closure magnetic state with very low stray fields in the rings.

74.78.-w 74.78.Fk 74.25.Dw

During the last years there has been a considerable effort to conceive and realize new superconducting devices that allow to modulate locally the magnetic fields practically at will zhu-2003 (). These systems rely essentially on either static arrays, such as pinning centers baert (); fiory (); martin (); raedts () and vortex-antivortex generators lange (); gillijns (), or components influencing the vortex dynamics like channels martinoli (); pastoriza (); silhanek-03 (); villegas-06 () and ratchets villegas (); clecio (). The ultimate motivation behind the manipulation of the local vortex density is to enhance the performance of superconductor-based devices by reducing the noise in squid-based systems roger-02 (); selders (), gaining control on superconducting THz emitters thz () or even providing a way to predefine the optical transmission through the system cond-mat ().

Unfortunately, for the majority of the components used in fluxonics devices, once they are created there is no margin for further modifications. In some cases this lack of flexibility becomes a limiting factor in the performance of the devices. For instance, a predefined ratchet system designed to effectively remove vortices from a specific location can work properly at low fields but ignores the inevitable reversed ratchet at higher fields thus making its functionality rather impractical villegas (); clecio (); lu (); silhanekapl (). A way to circumvent this shortcoming can be worked out by introducing magnetic pinning centers which have additional internal degrees of freedom not available in conventional nanostructured pinning sites or defects created via irradiation.

In this work we demonstrate that reversible and switchable guidance of the vortex motion can be achieved using a square array of magnetic square rings lying underneath a superconducting film. Transport measurements unambiguously show that the vortex dynamics is fully dominated by the magnetic landscape generated by the ring structures. When the magnetic unit cell is rotated 45 degrees off the Lorentz force the average vortex velocity follows the direction of the principal axis of the magnetic lattice rather than the driving force. This channeling effect can be easily suppressed when inducing a nearly isotropic pinning landscape by setting the square rings in a flux-closure state with very low stray fields.

The two samples used for this investigation consist of a 50 nm thick Al film evaporated on top of a checker-board () patterned array and a close-packed () square array of permalloy square rings with lateral size 1 m, line width 150 nm and thickness 25 nm [Fig. 1(a) and (b)]. The magnetic template is electrically separated from the superconducting film by a 5 nm Si buffer layer added to reduce proximity effects. The ring-shaped patterns and the transport bridge were fabricated with electron-beam lithography and lift-off technique on a silicon substrate. A superconducting coherence length nm for the sample and nm for the sample was estimated from the superconducting/normal phase boundary as determined by 10 normal state resistance criterion of a coevaporated reference film. The superconducting critical temperature at zero field for the CB and CP sample in the onion state (i.e. maximum stray field) is K and K, respectively.

Figure 1: (color on-line) Micromagnetic simulations of the out of plane component of the stray field at remanence for the checkerboard (a) and close-packed (b) arrays of square rings after applying a saturation field along their diagonal (onion state). The thick broken line indicates the magnetic unit cell for each pattern. The orientation of the external current (), field (), driving force (), and vortex drift () is clearly indicated. The lower panels schematically show a polar plot of the pinning force only considering the lattice symmetry (broken line), or the local magnetic pinning force (off centered circle), and including both the lattice and the local symmetry (solid line) for the checker-board (c) and the close-packed (d) arrays.

It has been shown recently that multiply connected magnetic structures make it possible to readily switch between different magnetic states imre (); Xiaobin-03 (). For the particular square geometry chosen in this work, an in-plane external field along the diagonal of the squares can induce six different domain distributions Vavassori-2003 (), namely four dipolar states (onion states) obtained at remanence after saturation and two flux-closure states with opposite chirality when the field is reduced from saturation to 36 mT and then set to zero. Interestingly, when these rings are placed in close proximity to a superconducting layer, either a strong vortex pinning or weak pinning can be obtained by simply switching from onion to flux-closure state, respectively. Although the influence of these magnetic templates on the vortex pinning is relatively well understood silhanek-06 (), little is known about their influence on the vortex dynamics.

In order to address this issue we carried out transport measurements recording simultaneously the electric field parallel () and perpendicular () to the external current. This allows us to estimate the direction (see Fig.1) of the average vortex motion with respect to the Lorentz force . Notice that in contrast to previous works studying the in-plane angular dependence of the vortex motion pastoriza (); silhanek-03 (); villegas-06 () where the orientation of the driving force was changed by combining two independent current sources, our ability to change the potential landscape in-situ allows us to use a simpler setup where the current orientation is kept constant. In all cases unwanted contributions to the measured signal like thermoelectric coupling, Hall resistance and misalignment of the transverse voltage contacts are small and have been properly taken into account.

Figure 2: (color on-line) Parallel () and transverse () electric field as a function of current density (upper row) and field (lower row) at . The first column [panels (a) and (b)] corresponds to the square rings magnetized at 45 as indicated in the inset of panel (a). The second column [panels (c) and (d)] presents the data for the square rings magnetized at 135, as indicated in the inset of panel (c). The third column [panels (e) and (f)] corresponds to the squares in the vortex state, as indicated in the inset of panel (e).

A representative set of experimental data for the sample [Fig.1(a)] and different magnetic states of the square rings is shown in Fig.2 at . Fig.2(a) shows the electric field-current density () characteristic for the case when the squares have been set in the onion state magnetized 45 away from the current direction [see inset in panel (a)] and setting an out-of-plane field where , with the flux quantum and the period of the lattice. From Fig.2(a) it can be seen that for low enough currents ( MA/cm) the vortex lattice remains pinned as no dissipation is detected in either direction (). Surprisingly, for currents higher than the critical current the direction of vortex motion does not coincide with the Lorentz force. More specifically since vortices move at an angle 45 as indicated in the inset of Fig.2(a). This is a clear evidence that the net displacement of the vortices is along the high symmetry axis of the magnetic pinning landscape indicated with a dashed line in Fig.1(a) and (b). Essentially this effect is the result of an anisotropic depinning force which reaches its maximum (minimum ) value perpendicular (parallel) to the channel direction defined by the principal axes of the magnetic landscape. In Fig.2(a), since the applied force is in this case 45 away from the channel direction, the critical current at which vortices start moving is determined by the condition . This situation persists up to MA/cm where and gradually separate from each other which means that now the component of the Lorentz force perpendicular to the channels is large enough to overcome , i.e. , and vortices flow along the direction of the Lorentz force. This allows us to estimate the ratio of depinning forces as .

By keeping a constant current density and progressively increasing the applied field it is possible to realize the same dynamic behaviors, as shown in Fig.2(b). For the low current density MA/cm vortices remain pinned up to a field where the vortex-vortex interaction becomes stronger than the vortex-pinning interaction. At higher fields clear commensurability effects can be seen for and . In this field region () the fact that indicates that vortices are guided by the pinning potential as described above.

These findings clearly demonstrate that under certain conditions the vortex motion can be effectively guided along the principal axes of the underlying pinning potential, but tell us little about the origin of this pinning potential. In other words, it is necessary to find out whether the guided vortex motion (GVM) is resulting from the magnetic landscape or from the corrugation of the superconducting film deposited on top of the squares. In order to answer this question we magnetized the square rings with an in-plane field 135 away from the current direction (0) as shown in the inset of Fig.2(c). In this way the topographic pinning remains the same whereas the magnetic pinning is rotated 90 with respect to the previous one. If the magnetic landscape is responsible for the GVM then a channeling along the -45 direction with respect to should be detected. This is in agreement with the results shown in Fig.2(c) and (d), where remains unchanged while reverses sign. The most convincing evidence however comes from panels (e) and (f) corresponding to the square rings in the flux-closure state. In this case the stray field generated by the domain walls in each square ring is reduced to a minimum level and as a consequence it is expected a minor influence from the magnetic landscape on the vortex dynamics. This is indeed corroborated by the lack of guided motion as is evidenced by the condition comment ().

Figure 3: (color on-line) Parallel () and transverse () electric field as a function of current density (upper panel) and field (lower panel) at for square rings in a close-packed array and magnetized at 45 as indicated in the inset.

It is worth to note that the fact that for the same geometry of the pinning lattice the GVM reverses sign when switching the dipole orientation from to , clearly points out the relevance of the local symmetry of the pinning centra. However, the question arises whether the geometry of the pinning lattice plays any role in the vortex guidance. In order to tackle this issue we have changed the lattice geometry from the checker-board pattern to a close-packed array [see Fig.1(b)]. In this way we maintain a square lattice of nearly point-like dipoles oriented at away from the current direction but change the lattice period from 1.5 m to 1.07 m. In addition, the unit cell of the dipolar lattice is rotated by as illustrated in Fig.1(a) and (b) with dashed lines.

The characteristics and the field dependence of the electric field for this sample with the square rings in the onion state are shown in Fig. 3 for similar experimental conditions shown in Fig. 2. Strikingly, the resulting measurements indicate the lack of vortex guidance (i.e. ) at all temperatures and fields explored experimentally. This clearly demonstrate the relevance of the particular orientation of the magnetic unit cell for the GVM.

In order to understand the origin of the guided vortex motion lets first hypothetically separate the different contributions arising from the lattice symmetry and from the local symmetry of the dipoles. This is schematically illustrated in panel (c) and (d) of Fig.1 for the checker-board and close-packed arrays, respectively. In order to unveil the lattice symmetry contribution we can assume that each dipole behaves as an effective isotropic pinning center, very much like a lattice of circular antidots or magnetic dots. In this case it has been previously shown silhanek-03 () that the minimum pinning force is along the principal axes of the pinning lattice [broken line in Fig.1(c) and (d)]. On the other hand, ignoring the presence of the lattice and just considering the pinning force produced by one single dipole gives rise to a strongly anisotropic pinning force carneiro () [schematically displayed as an off-centered circle in Fig.1(c) and (d)]. It is precisely the lack of inversion symmetry of the dipolar stray field which makes these structures efficient ratchet systems as well carneiro (); clecio-2007 (). The solid line in Fig.1(c) and (d) shows the resulting pinning force combining both effects: the four-fold symmetric lattice contribution and the non-symmetric local magnetic field. It is clear from this analysis that the presence of the non symmetric magnetic landscape lifts the degeneracy expected in the channeling of the sample [see Fig.1(c)] as experimentally observed. Moreover, the lack of guided motion in the sample for the used current direction is also expected since the minimum of the pinning force is along the direction of the Lorentz force [see Fig.1(d)], irrespective of the asymmetries introduced by the magnetic dipoles. It should be emphasized that although the above description helps to identify the basic mechanisms behind the guided motion of vortices, it is derived from the unphysical assumption that magnetic and non-magnetic contributions are separable.

In summary, we have reported experimental evidence indicating that the anisotropic potential landscape of a square array of magnetic dipoles is capable of directing the vortex motion 45 degrees away from the driving force direction. The multiple states of the used magnetic rings allow us to switch this direction ( for two different dipolar states and for the flux-closure state) resulting in a control of the transverse voltage signal. The contrasting behavior of the two different lattices explored (checker-board and close-packed arrays) demonstrate that the guidance of vortices is not only a consequence of the dipole orientation relative to the driving force direction but that also the orientation of the square lattice unit cell is important.

This work was supported by the Fund for Scientific Research-Flanders FWO-Vlaanderen, the Belgian Inter-University Attraction Poles IAP, the Research Fund K.U. Leuven GOA/2004/02, the European ESF NES programs, the U.S. NSF, grant ECS-0202780 and CNM ANL grants Nr.468 and Nr.470. A.V.S. is grateful for the support from the FWO-Vlaanderen. We acknowledge useful discussions with C.C. Souza Silva and G. Carneiro.


  • (1) B. Y. Zhu, F. Marchesoni, F. Nori, Physica E 18, 318 (2003); ibid. Phys. Rev. Lett. 92, 180602 (2004).
  • (2) A.T. Fiory, A.F. Hebard, S. Somekh, Appl.Phys.Lett. 32, 73 (1978).
  • (3) M. Baert, V.V. Metlushko, R. Jonckheere, V.V. Moshchalkov, and Y. Bruynseraede, Phys. Rev. Lett. 74, 3269 (1995).
  • (4) J. I. Martín, M. Velez, J. Nogués, and Ivan K. Schuller, Phys. Rev. Lett. 79, 1929 (1997).
  • (5) S. Raedts, A.V. Silhanek, M.J. Van Bael, R. Jonckheere, V.V. Moshchalkov, Physica C 404, 298 (2004)
  • (6) M. Lange, M.J. Van Bael, Y. Bruynseraede and V.V. Moshchalkov, Phys. Rev. Lett., 90, 197006 (2003).
  • (7) W. Gillijns, A.V. Silhanek, V.V. Moshchalkov, Phys. Rev. B74, 220509(R) (2006).
  • (8) P. Martinoli, Phys. Rev. B 17, 1175 (1978).
  • (9) H. Pastoriza, S. Candia, and G. Nieva, Phys. Rev. Lett. 83, 1026 (1999).
  • (10) A. V. Silhanek, L. Van Look, S. Raedts, R. Jonckheere, and V. V. Moshchalkov, Phys. Rev. B 68, 214504 (2003); ibid, Physica C 404, 340 (2004).
  • (11) J.E. Villegas, E.M. Gonzalez EM, M.I. Montero, I.K. Schuller, J.L. Vicent, J. Phys. Chem. Solids 67, 482 (2006).
  • (12) J. E. Villegas, S. Savel’ev, F. Nori, E.M. Gonzalez, J.V. Anguita, R. Garcia, J.L. Vicent, Science 302, 1188 (2003).
  • (13) C. C. de Souza Silva, J.V. de Vondel, M. Morelle, V.V. Moshchalkov, Nature 440, 651 (2006).
  • (14) R. Wordenweber, P. Lahl, P. Dymashevski, Physica C 369, 141 (2002).
  • (15) P. Selders, R. Wordenweber, Appl. Phys. Lett. 76, 3277 (2000).
  • (16) S. Savel’ev, V. Yampol’skii, A. Rakhmanov, F. Nori, Phys. Rev. B 72, 144515 (2005).
  • (17) Oleg L. Berman, Yurii E. Lozovik, Maria V. Bogdanova, Rob D. Coalson, cond-mat 0709.2198.
  • (18) Q.M. Lu, C. J. Olsen Reichhardt, C. Reichhardt, Phys. Rev. B 75, 054502 (2007).
  • (19) A. V. Silhanek, W. Gillijns, V. V. Moshchalkov, V. Metlushko, F. Gozzini, B. Ilic, W.C. Uhlig, J. Unguris, Appl. Phys. Lett. 90, 182501 (2007); W. Gillijns, A.V. Silhanek, V.V. Moshchalkov, C. J. Olson Reichhardt and C. Reichhardt, Phys. Rev. Lett. 99, 247002 (2007).
  • (20) A. Imre, E. Varga, L. L. Ji, B. Ilic, V. Metlushko, G. Csaba, A. Orlov, G. H. Bernstein, W. Porod, IEEE Trans. Mag. 42, 3641 (2006).
  • (21) X.B. Zhu, P. Grutter, V. Metlushko, B. Ilic, J. Appl. Phys. 93, 7059 (2003).
  • (22) P. Vavassori, M. Grimsditch, V. Novosad, V. Metlushko, B. Ilic, Phys. Rev. B 67, 134429 (2003).
  • (23) A. V. Silhanek, W. Gillijns, V. V. Moshchalkov, V. Metlushko, and B. Ilic, Appl. Phys. Lett. 89, 182505 (2006)
  • (24) The small but finite transverse field observed when the rings are set in the vortex state results from the fact that not all the rings switch to that state but some (vast minority) remain in the polarized state and therefore are still able to guide the vortices but with much less efficiency.
  • (25) G. Carneiro, Phys. Rev. B 72, 144514 (2005); ibid. Physica C 432, 206 (2005).
  • (26) C. C. de Souza Silva, A. V. Silhanek, J. van de Vondel, W. Gillijns, V. Metlushko, B. Ilic, and V. V. Moshchalkov, Phys. Rev. Lett. 98, 117005 (2007).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description