# Magnetic Skyrmion Lattice by Fourier Transform Method

###### Abstract

We demonstrate a fast numerical method of theoretical studies of skyrmion lattice or spiral order in magnetic materials with Dzyaloshinsky-Moriya interaction. The method is based on the Fourier expansion of the magnetization combined with a minimization of the free energy functional of the magnetic material in Fourier space, yielding the optimal configuration of the system for any given set of parameters. We employ a Lagrange multiplier technique in order to satisfy micromagnetic constraints. We apply this method to a system that exhibits, depending on the parameter choice, ferromagnetic, skyrmion lattice, or spiral (helical) order. Known critical fields corresponding to the helical-skyrmion as well as the skyrmion-ferromagnet phase transitions are reproduced with high precision. Using this numerical method we predict new types of excited (metastable) states of the skyrmion lattice, which may be stabilized by coupling the skyrmion lattice with a superconducting vortex lattice. The method can be readily adapted to other micromagnetic systems.

## I Introduction

Apart from the more familiar ordered magnetic phases, such as ferro- or antiferromagnetic, in recent years a plethora of more complex magnetically ordered states has been unveiled, both theoretically and experimentally. Among these states, a crucial role is played by magnetic skyrmions,Skyrme (1961) characterized by a whirling pattern of localized Heisenberg spins Bogdanov and Yablonskii (1989); Bogdanov and Hubert (1994); Barrett et al. (1995); Brey et al. (1995); Gervais et al. (2005); Mühlbauer et al. (2009); Neubauer et al. (2009); Münzer et al. (2010); Yu et al. (2010); Nagaosa and Tokura (2013). Magnetic skyrmions can be regarded as localized topological defects of the ferromagnetic state that can only be moved around or deformed in shape and size, but cannot be removed by a smooth deformation of the magnetization pattern due to micromagnetic constraints. They either exist as isolated defects, as, for example, in certain types of topological domain walls found in helical magnets, or they themselves give rise to an extended, energetically favorable phase of helical magnets, in the form of a skyrmion lattice.

Skyrmion lattices were first introduced by Klebanov Klebanov (1985) for neutron crystals, whereas Bogdanov and Yablonskii Bogdanov and Yablonskii (1989) discussed lattices of magnetic vortices (now called skyrmions lattices) for magnetic materials. Specific theoretical predictions for skyrmion lattices in e.g. , , and where made by Bogdanov and Hubert Bogdanov and Hubert (1994), and subsequently skyrmion lattices were reported experimentally in ,Mühlbauer et al. (2009); Neubauer et al. (2009) and in .Münzer et al. (2010); Yu et al. (2010) It is known from both theoretical predictions Bogdanov and Hubert (1994) and experimental evidence Nagaosa and Tokura (2013) that skyrmion lattices are 2D triangular lattices perpendicular to the external field that are translation invariant in the parallel direction.

Magnetic skyrmions arise due to the competition between ferromagnetic exchange and anisotropic exchange in the presence of an external magnetic field. While the exchange interaction tends to align all spins in one direction, Dzyaloshinsky-Moriya interactionDzyaloshinsky (1958); Dzyaloshinskii (1964); Moriya (1960) favors a non-collinear alignment. This effect results from the nature of the super-exchange mechanism in systems lacking inversion symmetry, where spin-orbit interaction rotates the direction of neighboring spins with respect to their perfect alignment governed by the exchange interaction only.

The model Hamiltonian for Heisenberg spins in the presence of a Dzyaloshinsky-Moriya interaction and an external magnetic field is given by,

(1) |

where is the ferromagnetic exchange constant, is the Dzyaloshinsky-Moriya coupling vector, is the external field, is Bohr’s magneton, and is a Heisenberg spin at site . The first term in Eq. (1) represents the Heisenberg exchange interaction, the second term the anisotropic Dzyaloshinsky-Moriya interaction, and the last term the Zeeman energy.

In the absence of an external magnetic field (or if the field is sufficiently weak so that the Zeeman energy can be neglected), the interplay of Dzyaloshinsky-Moriya and exchange interactions leads to a helical (spiral) spin order.Bak and Jensen (1980) For higher fields the Zeeman interaction becomes relevant and the spiral order gives way to a more complex magnetic state (the skyrmion lattice) until eventually at highest fields a ferromagnetic state is stabilized. Skyrmion magnetic order (see Fig. 1) occurs at intermediate external fields as a compromise between Dzyaloshinsky-Moriya interaction and Zeeman energy. In the center of a skyrmion, the spin points in the direction opposite to the magnetic field. When moving radially from the center towards the periphery, the direction of the spin rotates in a similar fashion as in the helical state, until it points in the direction of the magnetic field at the skyrmion boundary. In the regions between the skyrmions the spins stay aligned with the external magnetic field. By varying the inter-skrymion distance it is possible to arrive at a net energy gain resulting from both the Dzyaloshinsky-Moriya and Zeeman interactions, stabilizing a skyrmion lattice as opposed to a spiral or ferromagetic order.

The goal of this paper is to develop a fast and accurate numerical procedure for the calculation of skyrmion lattice or helical order by utilising a method similar to the one used by Brandt for the effective numerical implementation of superconducting vortex lattices.Brandt (1995) We additionally increase the speed of the numerical convergence by invoking a virial theorem in order to find the equilibrium lattice constants for the skyrmion lattice. This general procedure also allows for the study of metastable solutions and can be adapted to combine skyrmion lattices with superconducting vortex lattices. Such a combined system may stabilize new types of skyrmion lattices that otherwise would be metastable.

In Section II we present our analytical and numerical approach to the skyrmion lattice problem, including the results obtained for the phase diagram of a 2D non-centrosymmetric ferromagnet and the field dependence of the inter-skyrmion distance; in Section III we discuss alternative metastable solutions for the skyrmion lattice and the possible ways to stabilize them, focusing in particular on the honeycomb skyrmion lattice.

## Ii Skyrmion Lattice Phase

### ii.1 Euler-Lagrange Equations for a Magnetic System in Fourier Space

In this section we write the free energy functional for a two-dimensional non-centrosymmetric magnetic system in an external field and derive the corresponding Euler-Lagrange equations for the magnetization. Afterwards, we Fourier transform the obtained equations and solve them in Fourier space.

We start with the free energy functional written in the continuum limit. This approximation is legitimate as long as the magnetization varies on a length scale much larger than the crystal lattice spacing , i.e., when . In Cartesian coordinates, the free energy per unit area reads

(2) | |||||

where is the ferromagnetic exchange constant, is the Dzyaloshinsky-Moriya interaction, and is the external magnetic field applied. Here, , obeying , and , with , corresponding to a uniform external field in the -direction. We use the notation for the ratio between the Dzyaloshinsky-Moriya interaction and the ferromagnetic exchange coupling. As has dimensions of an inverse length, dimensionless coordinates result from . This leads to the dimensionless form of the free energy functional , Han et al. (2010)

where is the dimensionless magnetic field. For simplicity, we shall drop all tildes from now on, remembering that lengths are measured in units of and energies in units of .

Since the magnetization obeys the micromagnetic constraint at every point in space, we introduce a Lagrange multiplier field, , in order to enforce the aforementioned condition. The total functional density then reads

(4) | |||||

We now write the Euler-Lagrange equations for each component of the magnetization, ,

(5) |

with summation convention applied for , with . Explicitly, this leads to

(6) |

(7) |

(8) |

Being inspired by Brandt’s approach for Abrikosov vortices in type-II superconductors, Brandt (1995) we Fourier transform the magnetization components (and the Lagrange multiplier) and solve the problem in Fourier space. As the types of magnetic order promoted by the Dzyaloshinsky-Moriya interaction, namely the skyrmion lattice and helical order, are periodic structures, the Fourier approach is preferable. The discrete Fourier transform of is written as

(9) |

where is the Fourier coefficient for , and the discrete wave vectors are given by

(10) |

with and integer indices and, for a triangular lattice with lattice spacing , , , . In order to shorten the notation, we write

(11) |

and by analogy for the remaining components of the magnetization. For the Lagrange multiplier we write

(12) |

The Fourier transformed Euler-Lagrange equations, (6)–(8), then read

(13) |

(14) |

(15) |

The above equations can be rewritten in matrix notation as

(16) |

where the term on the right hand side is proportional to the external magnetic field,

(17) |

is the vector containing the Fourier coefficients of the magnetization components for a given ,

(18) |

and the matrix is given by

(19) |

with being an identity matrix, and

(20) |

Equation (16) can be solved analytically for via matrix inversion,

(21) |

yielding a magnetization functional that depends on the Lagrange multiplier field. This solution, however, does not necessarily obey the micromagnetic constraint . In order to enforce it, we must then solve the constraint equation numerically to find the appropriate Lagrange multiplier. In Fourier space, this amounts to solve the following equations for ,

(22) | ||||

(23) |

where we use the -dependence of the magnetization, , from Eq. (21). We find that the Lagrange multiplier takes the largest values in regions where the spins are aligned parallel or anti-parallel to the external field and the lowest values in regions where the spins are perpendicular to it (parallel to the -plane).

### ii.2 Skyrmion Lattice

Typical solutions obtained by the procedure described in the previous section are shown in Fig. 2. Keeping in mind that the applied field is directed in positive -direction, a typical skyrmion lattice solution appears with spins aligning along the field direction, , in the region between skyrmions, and in the opposite direction in the centers of skyrmions. Hence, the region between skyrmions is mostly ferromagnetic, while in moving away from the center of a skyrmion along any radial direction, the spin orientation rotates continuously within the plane perpendicular to the given radial direction, acquiring a component along the - plane, in a similar fashion as found in helical magnetic patterns. In this way, the system is able to gain energy from the Dzyaloshinsky-Moriya interaction, which favors the twisting of spins, while retaining a large region of ferromagnetic alignment, promoted by the Zeeman energy.

Shown in Fig. 2 are two self-consistently determined solutions for a triangular skyrmion lattice, one close to the phase transition to a spiral state (top panels), and one close to the phase transition to the homogeneous ferromagnetic state (bottom panels). In the first case the spatial modulation of the magnetization vector along a path going from one skyrmion center to a neighboring one approaches that of a spiral modulation, whereas in the latter case the skyrmions become isolated from each other and the inter-skyrmion regions expand. It is also interesting to observe, in these two extreme cases, the different spatial distribution of the -component of the magnetization within a single skyrmion, whose radius can be conventionally identified in the red dashed contour (). In the first case the area of negative magnetization (blue region, enclosed by the thick black contour) is approximately equal to the area of positive magnetization (red region, enclosed between the red dashed contour and the thick black contour). Hence, skyrmions in this case carry, as a whole, essentially no net magnetic moment, as in the case of a helix. On the other hand, in the vicinity of the ferromagnetic transition the area of positive magnetization is nearly three times that of negative magnetization, so that skyrmions do carry a finite magnetic moment along the external field direction. As we shall see later on, the shrinking of the negative magnetization region plays an important role in determining the magnetic field dependence of the inter-skyrmion distance.

Solutions to the Fourier transformed Euler-Lagrange equations, (13)–(15), satisfy the condition of stationarity of the functional within a given lattice geometry. However, there is no guarantee that such a solution is truly stationary with respect to arbitrary variations of the magnetization pattern and, even if this is the case, that the solution represents a global minimum of the functional and not just a local one. Hence, it is necessary to devise an additional procedure to distinguish stable solutions from those metastable or not stable at all.

We recall that the method described above is based on a Fourier expansion of the functional, where quantities such as the magnetization and the Lagrange multiplier are assumed to be periodic on a given lattice structure (triangular, for the time being) with lattice spacing . When seeking a solution to the Euler-Lagrange equations we must therefore remember that we are only exploring the subspace of magnetization profiles with a fixed given periodicity. The classes of solutions that we can obtain in this way include (but are not limited to) the helical magnetic order, the triangular skyrmion lattice, and the ferromagnetic order. Hence, in order to find the globally stable solution, i.e., the ground state of the system, we must first minimize, with respect to the lattice spacing , the free energy corresponding to a given class of solutions, and later on compare the optimal energies of each class of solutions, identifying the lowest one.

A typical profile of the free energy as a function of the lattice spacing is shown, for the skyrmion lattice class of solutions, in Fig. 3 (a). The free energy potential has a distinct minimum, corresponding to the optimal spacing of a triangular skyrmion lattice. The solutions shown in Fig. 2 are calculated for such free energy minima.

Although the value of the optimal spacing can be, in principle, extracted from this curve alone, one must be aware that its precise determination requires the evaluation of the potential at many points in around the minimum, due to the quadratic dependence of . On the other hand, there is a far more efficient procedure to extract the optimal spacing, and it is based on the virial theorem. Even more importantly, the virial theorem method allows one to discern whether or not the optimal spacing solution for a given lattice is truly stationary with respect to variations in the lattice geometry.

### ii.3 Virial theorem

In order to formulate the virial theorem for our problem, we rescale the coordinates inside the spatial integrals via a dimensionless parameter , writing . This leads to the following rescaling of the dimensionless free energy, Eq. (II.1),

(24) | ||||

where is the ferromagnetic exchange term, is the Dzyaloshinsky-Moriya term, is the Zeeman term, and . Since the above rescaling is just the result of a change of variables in the integrals, Eq. (24) holds for any value of . Hence, the right-hand side of Eq. (24) must be independent of and, accordingly, its derivative with respect to must vanish,

(25) | ||||

We can now use the arbitrariness in the value of and set in the above equation, in order to obtain the following relation between free energy terms,

(26) |

where we have used . The above equation holds for any magnetization profile , which need not be a minimum of the free energy functional. However, when represents a stationary solution of the Euler-Lagrange equations, the additional condition

(27) |

holds. Therefore, when evaluated on a genuinely stationary solution , Eq. (26) reads

(28) |

This is the virial theorem Doria et al. (1989) formulated for our problem, which can be rewritten in terms of the following “virial ratio”,

(29) |

satisfied by stationary magnetization profiles.

In Fig. 3 (b) we plot the virial ratio (29) against the lattice spacing. As expected, the minimum of the free energy potential coincides precisely with the value of spacing for which we have , indicating that the solution is a genuinely stable configuration. Furthermore, we clearly see that the virial ratio has a linear dependence on the lattice spacing, and this allows for a much more efficient and precise numerical determination of its optimal value , which we can now define through the condition , rather than through the minimum of the potential.

### ii.4 Phase Diagram and Critical Points

The free energy can be expressed in terms of the Fourier components via

This expression allows for a fast and precise evaluation of the free energy without the need to evaluate numerical derivatives.

In the following, we discuss solutions minimizing the free energy for different external magnetic fields, starting from . In order to get a better understanding of the transitions, it is also instructive to include the free energy functional of metastable solutions next to the stable ones. This is shown in Fig. 4, where energies of metastable helices, skyrmions, and ferromagnetic state are plotted alongside with the stable solutions.

Figure 4 also shows the critical values of for phase transitions. In our units the critical field for the helical-skyrmion phase transition is and the critical field for the skyrmion-ferromagnet phase transition is . This is in excellent agreement with the experimentally observed phase diagram of the 2D chiral magnet,Yu et al. (2010) and with previous theoretical investigations such as Monte-Carlo simulations Yi et al. (2009) and analytical approximate methods.Han et al. (2010) We find that the helical-skyrmion transition at is clearly first-order, as underlined by the possibility to extend, as metastable configurations, the skyrmion or spiral order into the neighboring phase (blue and red dotted lines in Fig. 4). The skyrmion-ferromagnet transition at , instead, appears to be second-order in the continuum limit. This is suggested by the impossibility to find metastable skyrmion lattices above and by the diverging behavior of the inter-skyrmion distance in approaching the transition, as discussed hereafter. However, we have to keep in mind that the continuum limit is only an approximation and that in any actual material the spins originate from a crystal lattice. Hence, for real materials the second-order transition could easily turn into a weakly first-order one.

The lattice spacing that provides an energetically optimal configuration varies with the external magnetic field, a behavior that recalls the analogy with Abrikosov vortex lattices in type-II superconductors. However, while the lattice spacing of a superconducting vortex lattice is a monotonic function of the external field, scaling as the inverse square root of the magnetic induction and being determined by the magnetic flux quantization,Abrikosov (1957) the optimal spacing of a skyrmion lattice, shown in Fig. 5, exhibits a non-monotonic dependence upon the external field. Starting from the helical-skyrmion transition, the top panels of Fig. 2) indicate that in the vicinity of skyrmions are close packed, so that the lattice spacing is roughly twice the skyrmion radius. Moreover, as already discussed in Section II.2, in this region a given skyrmion is characterized by approximately an equal amount of negative and positive , resulting in a zero net magnetic moment. If we now increase the magnetic field, the area of negative magnetization around the core of the skyrmion starts to shrink with respect to the positive one, giving the skyrmions a finite net magnetic moment along the field direction, and at the same time reducing the skyrmion radius. Since skyrmions at this stage are still close packed, the slightly smaller skyrmion radius results in a smaller lattice spacing, which reaches a minimum at . Upon further increasing the magnetic field, the skyrmion radius slowly continues to decrease. However, skyrmions cease to be close packed and become more and more isolated, separated from each other by increasingly large ferromagnetic domains, as shown in the bottom panels of Fig. 2). Accordingly, the skyrmion lattice spacing increases rapidly when moving towards the ferromagnetic phase. Exactly at the transition, the behavior shown in Fig. 5 seems to indicate a divergent lattice spacing, signaling a continuous vanishing of the skyrmion density. This scenario would imply a second-order phase transition, as mentioned earlier. Unfortunately, with our calculations we are not able to assert this with complete certainty due to the finite truncation of the Fourier series required by the numerical implementation. In fact, Eq. (10) shows that, for an increasingly large lattice spacing , the discrete vectors defining the Fourier series become closer and closer to each other, so that to achieve the same spatial resolution one has to include more and more vectors. In our numerical implementation, instead, we limited the number of vectors to a grid.

In addition to the skyrmion lattice spacing, in Figure 5 we present the field dependence of , the period of a helix. At zero magnetic field one can prove analytically that (in dimensionful units, ), so that the optimal value for a helix is Dzyaloshinsky (1958); Dzyaloshinskii (1964) . From the low- branch of the curve in Fig. 5 we can see that, at , obtained from our calculation is exactly equal to , showing that our calculations match the analytical predictions. As long as the field is turned on, the helix gets deformed, rotating faster in the regions where the magnetization is opposite to the field direction, and slower where the magnetization points along the field, so that the system can gain some Zeeman energy from the latter regions. Overall, this leads to a larger period of the helix, which thereby increases with the external field.

## Iii Metastable Skyrmion Lattices

In addition to the stable triangular skyrmion lattice phase illustrated above, we shall now discuss a few examples of alternative metastable skyrmion lattice geometries, and a possible way to stabilize them. If for a stable triangular lattice the virial ratio, and specifically its linear dependence on the lattice spacing, provided a numerically efficient and accurate tool to find the optimal spacing, in the case of metastable solutions it plays an even more important role. In fact, for lattice geometries different from the triangular one, there is no guarantee that the solution obtained by minimising the free energy with respect to the lattice spacing represents a true metastable solution. Changing the lattice spacing is indeed only one of the possible ways to vary the underlying lattice structure, the others being a change in the ratio and in the angle between the unit vectors of the lattice—namely, a change in the lattice geometry. As discussed in Section II.3, any truly stationary solution of the Euler-Lagrange equations, including metastable skyrmion lattices, should fulfill the virial theorem. Hence, we can use this criterion to distinguish true metastable solutions from those that are local minima of the free energy only within the subset of solutions with a given lattice symmetry.

The first type of metastable skyrmion lattice that we are going to discuss is a square lattice, shown in the top panels of Fig. 6. This pattern can be realized, as an alternative to the triangular lattice, in the whole skyrmions phase between and , and it represents a genuine metastable solution in the sense that it is a local minimum of the free energy functional with respect to arbitrary variations of the magnetization pattern. This is confirmed by the virial ratio, which takes the value of at the minimum of the free energy as a function of the lattice spacing. In the phase diagram of Fig. 4, the energy of the square lattice (not shown to avoid an overburdening of the plot) is only slightly higher than the ground state energy, lying between the triangular (blue full curve) and honeycomb (black double-dot-dashed curve) skyrmion lattice energies. In approaching the skyrmion-ferromagnet transition, all skyrmion energies merge together at . Such a behavior is precisely what one would expect from a second-order phase transition, characterized by a vanishing skyrmion density (or, equivalently, a diverging lattice spacing). Indeed, it is clear that for nearly isolated skyrmions the details of the lattice symmetry become unimportant and the energy of the system is controlled only by the skyrmion density.

The second type of metastable phase that we would like to present is, instead, a peculiar metastable solution that was found while investigating the system, namely a honeycomb skyrmion lattice, depicted in the bottom panels of Fig. 6. As it is well known for honeycomb lattices, the underlaying Bravais lattice for this solution is still a triangular lattice, however the unit cell contains now two skyrmions and the lattice spacing is no longer the distance between the centers of nearest skyrmions, but the distance between the centers of honeycombs. In Fig. 7 we show the free energy of a honeycomb skyrmion lattice as a function of the lattice spacing for a given . As in the case of a triangular lattice (shown in the same figure for comparison), the free energy has a clear minimum at a certain optimal spacing. However, in stark contrast with the triangular and square lattice solutions presented earlier, the honeycomb lattice does not fulfill the virial theorem at the minimum of the corresponding free energy. This behavior entails a strong message, telling us that this solution cannot be a true local minimum of the free energy functional. Indeed, the free energy of such a configuration does have a local minimum as a function of the lattice spacing, but only as long as one considers lattice configurations with trigonal symmetry. Instead, if the geometry of the lattice is allowed to change, this solution can relax into a lower energy configuration. At first, one might think that configurations with high symmetries, such as the square or trigonal one, should be protected against deformations of the lattice geometry, so that the corresponding magnetization patterns should always be at least stationary points of the free energy functional. However, this is true only for simple Bravais lattices such as the square or triangular lattices, but no longer holds in the case of Bravais lattices with a basis such as the honeycomb one. In the latter case, for example, one realizes that the honeycomb lattice can be smoothly turned into a square lattice, which has a lower energy, by modifying the angle and ratio between the unit vectors in a certain direction. A modification of the angle and ratio in the opposite direction, instead, leads to a very different configuration with a higher energy. The honeycomb lattice would then relax into the square lattice if the lattice symmetry is left unconstrained. On the other hand, it still represents a valid metastable solution if the space of solutions is constrained, by some external mechanism, to the subset of configurations with trigonal symmetry.

As in the case of a square skyrmion lattice, the free energy of this peculiar metastable solution lies just above the free energy of the stable ground state configuration and, in approaching , the difference between them becomes almost negligible. Hence, while this type of solution has not been experimentally observed in pure skyrmion materials, there is hope that the honeycomb skyrmion lattice might be stabilized by some sort of magnetic perturbation that establishes an underlying periodic structure with trigonal symmetry. In particular, by looking at the spin texture of this solution, we realize that skyrmions are arranged so as to surround and enclose large ferromagnetic domains. These domains are themselves arranged in a triangular lattice with vertices located at the centers of honeycombs. Since we already know that ferromagnetic domains are stabilized by a large magnetic field, a promising way of stabilising this honeycomb skyrmion lattice is to make the spins experience an “external” magnetic field that is no longer constant, but modulated periodically so that it is largest at the vertices of a triangular lattice.

A noteworthy example of such a spatial pattern is given by the microscopic magnetic field a type-II superconductor in the Abrikosov vortex phase, Abrikosov (1957) where the magnetic field is largest at the vortex cores and weaker in the interstitial superconducting region between vortices. Hence, provided that the length and magnetic field scales of the skyrmion material and type-II superconductor are compatible, we argue that a promising way to stabilize the honeycomb skyrmion lattice is to realize a bi- or tri-layer consisting of a non-centrosymmetric ferromagnet and a strong type-II superconductor, electromagnetically coupled, where the Abrikosov vortices would be pinned at the centers of the skyrmion honeycombs. Schematic image of such a system can be found in Fig. 8.

## Iv Summary

In this paper we have introduced a numerically fast and stable approach for modeling magnetic systems by combining a Fourier-transform method with a virial theorem. We have demonstrated the accuracy of the method by applying it to a magnetic skyrmion lattice, reproducing the optimal skyrmion spacing and critical fields.

We then have discussed a novel metastable honeycomb lattice solution, and suggested to combine a skyrmion lattice with a superconducting vortex lattice in order to stabilize this solution. The new phase arranges superconducting vortices in the inter-skyrmion regions, whereas skyrmions are placed in the regions where the local magnetic field is suppressed due to superconductivity. Both of these requirements can be fulfilled in a state hosting a skyrmion honeycomb lattice in combination with a triangular Abrikosov vortex lattice. This system would also lead to interesting potential ways of manipulating skyrmions, which could be moved around by drifting the Abrikosov vortices.Rabinovich et al. (2018); Reichhardt and Reichhardt (2018); Becerra and Milošević (2017) Moreover, the inhomogeneous magnetic pattern of skyrmions could be able to induce odd-frequency spin-triplet correlations in the superconductor already in a bilayer system, without the need of two separate ferromagnetic layers as in the case of superconducting spin valves.Kalenkov et al. (2011); Yokoyama and Linder (2015)

Finally, the proximity between skyrmion textures and either s-wave or p-wave superconductors is a promising tool to generate and manipulate topologically non-trivial excitations such as Majorana bound states and tuneable Weyl points.Hals et al. (2016); Yang et al. (2016); Zhang et al. (2016); Pershoguba et al. (2016); Pöyhönen et al. (2016); Takashima and Fujimoto (2016)

###### Acknowledgements.

ME acknowledges support by EPSRC grant EP/N017242/1.## References

- Skyrme (1961) T. H. R. Skyrme, Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 260, 127 (1961).
- Bogdanov and Yablonskii (1989) A. N. Bogdanov and D. A. Yablonskii, Journal of Experimental and Theoretical Physics 68, 101 (1989).
- Bogdanov and Hubert (1994) A. Bogdanov and A. Hubert, Journal of Magnetism and Magnetic Materials 138, 255 (1994).
- Barrett et al. (1995) S. E. Barrett, G. Dabbagh, L. N. Pfeiffer, K. W. West, and R. Tycko, Phys. Rev. Lett. 74, 5112 (1995).
- Brey et al. (1995) L. Brey, H. A. Fertig, R. Côté, and A. H. MacDonald, Phys. Rev. Lett. 75, 2562 (1995).
- Gervais et al. (2005) G. Gervais, H. L. Stormer, D. C. Tsui, P. L. Kuhns, W. G. Moulton, A. P. Reyes, L. N. Pfeiffer, K. W. Baldwin, and K. W. West, Phys. Rev. Lett. 94, 196803 (2005).
- Mühlbauer et al. (2009) S. Mühlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, and P. Böni, Science 323, 915 (2009).
- Neubauer et al. (2009) A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, P. G. Niklowitz, and P. Böni, Phys. Rev. Lett. 102, 186602 (2009).
- Münzer et al. (2010) W. Münzer, A. Neubauer, T. Adams, S. Mühlbauer, C. Franz, F. Jonietz, R. Georgii, P. Böni, B. Pedersen, M. Schmidt, A. Rosch, and C. Pfleiderer, Phys. Rev. B 81, 041203 (2010).
- Yu et al. (2010) X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Y. Matsui, N. Nagaosa, and Y. Tokura, Nature 465, 901 EP (2010).
- Nagaosa and Tokura (2013) N. Nagaosa and Y. Tokura, Nature Nanotechnology 8, 899 EP (2013).
- Klebanov (1985) I. Klebanov, Nuclear Physics B 262, 133 (1985).
- Dzyaloshinsky (1958) I. Dzyaloshinsky, Journal of Physics and Chemistry of Solids 4, 241 (1958).
- Dzyaloshinskii (1964) I. E. Dzyaloshinskii, Journal of Experimental and Theoretical Physics 19, 960 (1964).
- Moriya (1960) T. Moriya, Phys. Rev. 120, 91 (1960).
- Bak and Jensen (1980) P. Bak and M. H. Jensen, Journal of Physics C: Solid State Physics 13, L881 (1980).
- Brandt (1995) E. H. Brandt, Reports on Progress in Physics 58, 1465 (1995).
- Han et al. (2010) J. H. Han, J. Zang, Z. Yang, J.-H. Park, and N. Nagaosa, Phys. Rev. B 82, 094429 (2010).
- Doria et al. (1989) M. M. Doria, J. E. Gubernatis, and D. Rainer, Phys. Rev. B 39, 9573 (1989).
- Yi et al. (2009) S. D. Yi, S. Onoda, N. Nagaosa, and J. H. Han, Phys. Rev. B 80, 054416 (2009).
- Abrikosov (1957) A. A. Abrikosov, Journal of Experimental and Theoretical Physics 5, 1174 (1957).
- Rabinovich et al. (2018) D. S. Rabinovich, I. V. Bobkova, A. M. Bobkov, and M. A. Silaev, arXiv:1808.06102 (2018).
- Reichhardt and Reichhardt (2018) C. Reichhardt and C. J. O. Reichhardt, arXiv:1808.08214 (2018).
- Becerra and Milošević (2017) V. F. Becerra and M. V. Milošević, Physica C: Superconductivity and its Applications 533, 91 (2017), ninth international conference on Vortex Matter in nanostructured Superdonductors.
- Kalenkov et al. (2011) M. S. Kalenkov, A. D. Zaikin, and V. T. Petrashov, Phys. Rev. Lett. 107, 087003 (2011).
- Yokoyama and Linder (2015) T. Yokoyama and J. Linder, Phys. Rev. B 92, 060503 (2015).
- Hals et al. (2016) K. M. D. Hals, M. Schecter, and M. S. Rudner, Phys. Rev. Lett. 117, 017001 (2016).
- Yang et al. (2016) G. Yang, P. Stano, J. Klinovaja, and D. Loss, Phys. Rev. B 93, 224505 (2016).
- Zhang et al. (2016) L.-F. Zhang, V. F. Becerra, L. Covaci, and M. V. Milošević, Phys. Rev. B 94, 024520 (2016).
- Pershoguba et al. (2016) S. S. Pershoguba, S. Nakosai, and A. V. Balatsky, Phys. Rev. B 94, 064513 (2016).
- Pöyhönen et al. (2016) K. Pöyhönen, A. Westström, S. S. Pershoguba, T. Ojanen, and A. V. Balatsky, Phys. Rev. B 94, 214509 (2016).
- Takashima and Fujimoto (2016) R. Takashima and S. Fujimoto, Phys. Rev. B 94, 235117 (2016).