Stable and manipulable Bloch point

The prediction of magnetic skyrmions being used to change the way we store and process data has led to materials with Dzyaloshinskii-Moriya interaction coming into the focus of intensive research. So far, studies have looked mostly at magnetic systems composed of materials with single chirality. In a search for potential future spintronic devices, combination of materials with different chirality into a single system may represent an important new avenue for research. Using finite element micromagnetic simulations, we study an FeGe disk with two layers of different chirality. We show that for particular thicknesses of layers, a stable Bloch point emerges at the interface between two layers. In addition, we demonstrate that the system undergoes hysteretic behaviour and that two different types of Bloch point exist. These ‘head-to-head’ and ‘tail-to-tail’ Bloch point configurations can, with the application of an external magnetic field, be switched between. Finally, by investigating the time evolution of the magnetisation field, we reveal the creation mechanism of the Bloch point. Our results introduce a stable and manipulable Bloch point to the collection of particle-like state candidates for the development of future spintronic devices.

The prediction of magnetic skyrmions being used to change the way we store and process data has led to materials with Dzyaloshinskii-Moriya interaction coming into the focus of intensive research. So far, studies have looked mostly at magnetic systems composed of materials with single chirality. In a search for potential future spintronic devices, combination of materials with different chirality into a single system may represent an important new avenue for research. Using finite element micromagnetic simulations, we study an FeGe disk with two layers of different chirality. We show that for particular thicknesses of layers, a stable Bloch point emerges at the interface between two layers. In addition, we demonstrate that the system undergoes hysteretic behaviour and that two different types of Bloch point exist. These 'head-to-head' and 'tail-to-tail' Bloch point configurations can, with the application of an external magnetic field, be switched between. Finally, by investigating the time evolution of the magnetisation field, we reveal the creation mechanism of the Bloch point. Our results introduce a stable and manipulable Bloch point to the collection of particle-like state candidates for the development of future spintronic devices.
In recent years, materials with the Dzyaloshinskii-Moriya interaction (DMI) 1,2 came back into the focus of intensive research, mostly due to the discovery of magnetic skyrmions 3,4 that promise to revolutionise the way we store and process data 5 . Contrary to the ferromagnetic exchange which tends to align all spins parallel to each other, the DMI energy is minimised when neighbouring spins in the crystal lattice are perpendicular. Their mutual competition results in a twist between neighbouring spins which eventually leads to the helical magnetic order 6 over the entire crystal lattice. Although the angle between neighbouring spins, and consequently the helical period, depend on the relative strength of DM compared to the ferromagnetic exchange energy, the handedness of the emerged helical order in a helimagnetic material depends on the material's chirality -sign of the DM energy constant. So far, studies have looked mostly at magnetic systems composed of materials with single chirality. In a search for potential future spintronic devices and new particle-like magnetisation states, combination of materials with different chirality may represent an important new avenue for research.
To guide our investigation we start with a thought experiment. In the absence of chiral interactions, in thin film disks, a vortex configuration can emerge due to the competition between ferromagnetic exchange and demagnetisation energies. A vortex can have two different kinds of handedness as well as two different core polarisations. This results in four possible states of a magnetic vortex in a disk as shown in Fig. 1(a). In order to avoid possible cacophony with the terminology established in previous studies, throughout this work, we refer to the +z or −z core orientation of the vortex state as its polarisation, and to the clockwise or counterclockwise chirality/winding as its handedness. Let us now imagine we have two disks where vortices with the same handedness were formed, but with different polarisation. We stack these on top of each other as though they were parts of a much thicker disk as shown in Fig. 1(b). This would result in a magnetisation configuration that is continuous at all points of the sample, except at the centre, where two vortex cores point in the opposite direction 7 . This discontinuity in the magnetisation field is called a Bloch point 8,9 . There are two possible discontinuous polarisations of the cores, and consequently two different types of a Bloch point, as we show in Fig. 1(b). A Bloch point has a high magnetic energy, especially because of the high ferromagnetic exchange energy at the point where two spins point in the opposite direction. Consequently, the entire system tends to reduce its energy by expelling the high energy Bloch point out of the sample. In our system, this would occur by changing the polarisation of either of the two vortices so that both vortex cores have the same polarisation but keep the same handedness.
In the system of two stacked disks we analysed, the reorientation (change of polarisation) of one of the vortices to expel the Bloch point is possible because in a magnetic vortex both polarisations are allowed for a particular handedness. However, if the Dzyaloshinskii-Moriya energy is present in the system, the polarisation and handedness of a vortex (more precisely, vortex-like state) are always related to each other and only two possible states of a magnetic vortex are allowed 10 . In our thought experiment, this would imply that changing the vortex core polarisation while keeping the handedness fixed in order to expel the Bloch point from the system is not possible. Therefore, we can assume that a Bloch point will remain present in the sample. In this work, we explore this idea.
Using finite element micromagnetic simulations, we demonstrate the existence and stability of a Bloch point at the interface of two materials with different chirality. In addition, by exploring the hysteretic behaviour, we find two different Bloch point configurations and show how by using an external magnetic field, one can switch between them. Finally, by analysing the time evolution of a uniformly magnetised system, we reveal the creation mechanism of a Bloch point. We believe that the existence and stability of a Bloch point in the studied system, as well as the possibility of its manipulation, introduces a new particle-like state to the collection of candidates for the development of future spintronic devices.

Results
Stability. We begin our study by exploring whether a Bloch point emerges in the system shown in Fig. 1(c), obtained as a result of a thought experiment we conducted in the introduction. In order to base our study on realistic material parameters and therefore encourage the experimental verification of our predictions, we use cubic B20 FeGe material parameters (see Methods section). We set the DMI constant in the bottom layer to be negative (D = −|D|) and fix its thickness h bottom to be 20 nm. In the top layer we impose a positive DMI constant (D = |D|) and vary its thickness h top between 2 nm and 18 nm in steps of 1 nm. The diameter of the disk sample is 150 nm and no external magnetic field is applied. For each thickness of the top layer, we begin by initialising the system with a uniform magnetisation configuration in the positive out-of-plane (+z) direction. This initial configuration corresponds to the fully saturated magnetisation configuration feasible in an experimental setup. After that, we relax the system by integrating a set of dissipative time-dependent equations until the condition of vanishing torque (m × H eff ) is satisfied. After the relaxed (equilibrium) state is obtained, we compute a three-dimensional topological charge in the top and the bottom layer individually using 3 as suggested by Lee et al. 11 . Because Q is proportional to the volume (and accordingly the thickness) of individual layers, we normalise it by the thickness =  Q Q h / , in order to make the values mutually comparable between different layers. The only reason we use normalised topological charge is to track changes in the magnetisation of individual layers in the system. We show in Fig. 2(a) how normalised topological charges  Q of the equilibrium state in individual layers depend on the thickness of the top layer h top . Although the topological charge of the bottom layer  Q bottom remains relatively constant, a sharp decrease of  Q top between 8 nm and 9 nm thicknesses of the top layer is evident. Consequently, we assume that depending on the thickness of the top layer, two different configurations can emerge as equilibrium states in the studied system. Subsequently, we explore these two magnetisation configurations.
To explore the state that emerges for h top ≤ 8 nm, we visualise the magnetisation obtained for h top = 5 nm in Fig. 2(b) and for h top ≥ 9 nm magnetisation configurations, we show in Fig. 2(c) the state obtained with t top = 15 nm. We show the magnetisation vector field together with the m z scalar field after we cut out a slice from the sample in order to be able to look inside and inspect its centre. We observe that for both states, vortex-like configurations with clear chiralities and polarisations emerge in individual layers. In the literature, the vortex-like www.nature.com/scientificreports www.nature.com/scientificreports/ states formed as a consequence of DM energy are also called quasi-ferromagnetic 12 , edged vortex 13 , or incomplete skyrmion 10,14,15 states. In both states we identified, the polarisation and the handedness of the vortex-like states in bottom layers are the same and related by the right hand rule -characteristic for the negative DM energy constant. On the other hand, the vortex-like states in the top layer differ for two equilibrium states we obtained.
The vortex-like state in the top layer with 5 nm thickness shown in Fig. 2(b), has the same handedness and polarisation as the thicker bottom layer. However, this is against the left hand rule which establishes the relationship between the handedness and polarisation for a positive DMI constant. Therefore, the DM energy density in the top layer is significantly higher than in the bottom layer and one can ask why the polarisation of the vortex-like state in the top layer is not reversed in order to minimise its energy. Because the energy contribution of the top layer to the total energy of the system is relatively small due to its thickness, it is energetically cheaper for the top layer to follow the magnetisation of the bottom layer (violate the left hand rule, and eventually have the higher energy density), than to reorient itself, which would result in an energetically more expensive Bloch point at the interface between two layers. Now, we look at magnetisation configuration in the system with h top = 15 nm, shown in Fig. 2(c). The chiralities of vortex-like states in both layers are the same, but their polarisations are opposite and a Bloch point emerges between them. For h top ≥ 9 nm and above, the energy contribution of the top layer to the total energy of the system becomes significant. Therefore, the top layer now does not follow the polarisation of the bottom layer, but reverses its polarisation to follow the left hand rule. It is now energetically cheaper for the system to host a Bloch point than to tolerate the top layer following the polarisation of the bottom layer and having high energy density.
The location marked by a circle in the middle of the sample in Fig. 2(c) identifies a Bloch point: a noncontinuous singularity in the magnetisation vector field where the magnetisation magnitude vanishes to zero 8,9 . Because micromagnetic models assume constant magnetisation magnitude, the precise magnetisation configuration at the Bloch point cannot be obtained using micromagnetic simulations 16 . However, it is known how to identify the signature of a Bloch point in such situations: the magnetisation direction covers any sufficiently small closed surface surrounding the Bloch point exactly once 17,18 . We illustrate this property in Fig. 2(d).

Hysteretic behaviour.
In the previous section we demonstrated the existence and stability of a Bloch point in the proposed system at zero external magnetic field. In all relaxation simulations so far, we started the relaxation from the uniform magnetisation configuration in the positive out-of-plane (+z) direction which resulted in the vortex-like state in the bottom layer with the +z polarisation. Furthermore, for h top ≥ 9 nm, top layer relaxed in the vortex-like configuration with negative −z polarisation. This resulted in a head-to-head Bloch point as we showed in Fig. 2(d). At this point, we can assume that if we started our relaxations from the uniform configuration in the −z direction, the polarisations of the bottom and the top layer would be in −z and +z directions, respectively. Consequently, a tail-to-tail Bloch point would emerge. Motivated by this assumption and knowing that uniform magnetisation configurations can be obtained by applying a sufficiently strong external magnetic field, we now explore the hysteretic behaviour of the proposed system in order to verify our assumption.
We fix the thickness of the bottom layer with D < 0 to be 20 nm as before and set the thickness of the top layer with D > 0 to be 10 nm so that we can expect a Bloch point to emerge. The diameter of the disk is 150 nm. We apply an external magnetic field in the out-of-plane +z direction and vary it between 1 T and −1 T in steps of 0.1 T. We simulate the hysteretic behaviour in the standard way by relaxing the system to an equilibrium state www.nature.com/scientificreports www.nature.com/scientificreports/ after changing the external magnetic field, and then using the resulting state as the starting point for a new energy minimisation. We show the average out-of-plane magnetisation component 〈m z 〉 as a function of an external magnetic field μ 0 H in Fig. 3(a). From this plot, we see that hysteretic behaviour of the simulated system is evident and now we investigate what are the equilibrium magnetisation states at zero external magnetic field.
We start our hysteresis simulations with 1 T external magnetic field, which saturates the sample magnetisation in the +z direction. By reducing the external magnetic field in steps of 0.1 T, we reach μ 0 H = 0 and plot the equilibrium magnetisation vector field together with m z scalar field in the xz cross section containing the disk sample centre in Fig. 3(b). We see that bottom and top layers relax in vortex-like magnetisation configurations with the same handedness, but with +z and −z polarisations, respectively. This results in a head-to-head Bloch point as we observed in the previous section. Now, we keep reducing the external magnetic field and when we reach −1 T, the sample is again fully saturated, but this time in the −z direction. Starting from this magnetisation state, we increase the field in same steps until we reach μ 0 H = 0, and plot the xz cross section of the equilibrium state in Fig. 3(c). Bottom and top layers are now hosting vortex-like states with the same handedness, but now with −z and +z polarisations, respectively. Therefore, a tail-to-tail Bloch point configuration is formed. Finally, we keep increasing an external magnetic field until we reach the 1 T value from which we started the hysteresis loop simulation. We show the evolution of the magnetisation field in the xz cross section through the entire hysteresis loop in Supplementary Video 1. We demonstrated that the studied system undergoes hysteretic behaviour and that two different types of Bloch point configurations can emerge. In addition, by applying an external magnetic field, one can manipulate the Bloch point and switch between head-to-head and tail-to-tail configurations.
Creation mechanism. Now, we explore how the Bloch point is created in the sample. We simulate the same sample we used in hysteresis simulations, but this time we do not apply an external magnetic field. We start by initialising the system in the positive out-of-plane +z direction. After that, we let the magnetisation evolve towards its equilibrium state and record the magnetisation configuration every 1 ps. We set the Gilbert damping value to the experimentally measured value 19 α = 0.28. In Fig. 4(a)   and 〈m z 〉 bottom = 〈m z 〉 top = 1 at t = 0 ps, which corresponds to the initial uniform state, all quantities evolve. In order to explore how the Bloch point is created in the sample, we show snapshots at different points in time in Fig. 4(b). We show the magnetisation vector field together with the m z scalar field in the xz cross section containing the centre of the sample.
By locating a Bloch point in the magnetisation fields (denoted in magnetisation field plots) in Fig. 4(b), we explore its creation mechanism. Starting from the uniform configuration at t = 0 ps, until approximately t = 55 ps www.nature.com/scientificreports www.nature.com/scientificreports/ the magnetisations in both layers evolve so that vortex-like states with the same handedness and polarisation are formed. At t = 58 ps we identify that the polarisation of the top layer begins to reverse and a Bloch point is formed at the top boundary of the sample. Furthermore, Bloch point propagates downwards -towards the interface between layers, until it eventually reaches its final position at t = 80 ps. We notice that the Bloch point creation mechanism is based on the polarisation (core orientation) reversal of the top layer which is the same as the reversal mechanism for skyrmionic states in confined helimagnetic nanostructures 10 . More precisely, in the process of reversal of the top layer polarisation, the Bloch point is created in the same way -at the boundary, propagating through the thickness of the sample and sitting at the interface between layers rather than being expelled out of the sample. We show the entire creation mechanism in the xz cross section of the sample in Supplementary Video 2.

Discussion
Using finite element micromagnetic simulations, we studied nanostructures composed of two layers with different handedness. We found that for certain thicknesses of individual layers, a stable Bloch point can emerge at the interface between the two layers. By applying an external magnetic field we found that the system undergoes hysteretic behaviour and two different types of Bloch point configurations exist. In the first, the magnetisation at the centre of each layer points towards the interface -the head-to-head configuration. In the second, the magnetisation points away from the interface towards the outer surface of the disk -the tail-to-tail configuration. We demonstrated the switch between these two different types of Bloch point using an external magnetic field. Finally, we simulated the time evolution in order to determine the mechanism by which the Bloch point is created, showing that the Bloch point is created at the boundary of the thinner layer and then propagates through the thickness of the sample until it reaches the boundary between layers where it sits. Now, we discuss our findings focusing on four main topics: (i) the possibility of fabrication, (ii) the suitability of micromagnetic simulations for the study of Bloch points, (iii) similarities and differences of our studied system with previous works, and (iv) implications of our results.
We based our simulations on realistic FeGe material parameters 10 for a realistically sized disk in order to encourage the experimental verification of our predictions. In addition, in all of our simulations we used a uniform magnetisation configuration as an initial magnetisation state which is feasible to obtain in an experimental setup by applying a strong field. More precisely, applying an external magnetic field of approximately 1 T, and then removing it, should allow the system to relax into the magnetisation configuration containing a Bloch point. However, before our predictions can be validated, one should be able to grow a helimagnetic thin film which consists of two layers with different chirality. More precisely, in the growth of a thin film, it is necessary to be able to switch the Dzyaloshinskii-Moriya energy constant D once a layer of sufficient thickness has been grown. Recent research by Spencer et al. 20 demonstrated that by changing the Co content y in Fe 1−y Co y Ge compound grown on Si(111), it is possible not only to modify the value of the DMI constant, but also to change its sign. This finding may contribute to the growth of a thin film system, composed of two layers with different chirality in order to fabricate our proposed system hosting a stable and manipulable Bloch point.
In this work, we used finite elements to discretise the continuous magnetisation fields. Because of the imposition of a constant magnetisation saturation in micromagnetics, this method is inaccurate to determine the Bloch point structure 16 . This is because at the Bloch point the magnetisation vanishes to zero 8,9 . However, it is well known that micromagnetics can provide us enough information to determine whether the Bloch point is present in the sample and its position, by analysing the magnetisation field on a sufficiently small closed surface surrounding the Bloch point. More precisely, it is necessary for the magnetisation direction to cover a sufficiently small closed surface surrounding the Bloch point exactly once 17,18 . Using this condition we were able to determine the existence and the position of a Bloch point. However, because of the imposed constant magnetisation saturation, we did not attempt to determine its internal structure, which was not the topic of this work. As an example of suitability of micromagnetics to determine the existence of a Bloch point, Rybakov et al. 21 predicted, using micromagnetic simulations, the existence of a chital bobber -a hybrid particle emerging at the surface of chiral magnets composed of both smooth magnetisation configuration and a Bloch point. The existence of a Bloch point as a part of chiral bobber was then experimentally confirmed in FeGe 22,23 , demonstrating the validity of the micromagnetic predictions.
To some extent similar concept was introduced by Zhang et al. 24 , where two ferromagnetic layers hosting skyrmions were antiferromagnetically coupled in order to suppress the skyrmion Hall effect, when skyrmions are driven using spin-polarised currents. However, in that system, no Bloch point can emerge between two skyrmions with different polarisations because they are physically separated by an insulating spacer of non-zero thickness. On the other hand, we achieve different polarisations (in our case, of vortex-like states) by stacking helimagnetic layers with different chirality. Because of that, there is no insulating spacer between the two layers with different chiralities and a Bloch point can be identified.
Our discovery of a stable and manipulable Bloch point in a planar magnetic system could add the Bloch point to the collection of already existing and well studied particle-like magnetic configurations that could possibly change the way we store and process data. Now, we discuss and speculate about the possibilities of how the Bloch point we explored in this work can be employed in spintronic applications. We found two different Bloch point configurations that can be switched using an external magnetic field, thus one can propose that a single information bits (0 or 1) can be encoded using the type of a Bloch point (HHBP or TTBP) -similar to the bit-patterned media. However, from the hysteresis simulations, we saw that the average out-of-plane component of magnetisation 〈m z 〉 at zero field is almost zero due to the symmetry of the magnetisation field. This poses a difficulty for reading such information bits. A possible solution is to increase the thickness of the bottom layer, which does not affect the stability of the Bloch point, but increases 〈m z 〉 at zero field. Furthermore, our initial studies show that the Bloch point does not emerge only is disk geometries, but also in a wide variety of different geometries, such as www.nature.com/scientificreports www.nature.com/scientificreports/ nanostrips. In this case, the manipulation of a Bloch point is not only limited to its type, but also its position in the sample. This manipulation can be achieved using an external excitation, for instance, external magnetic fields or spin-polarised currents. However, a Bloch point coupling to such external excitations is beyond the scope of this work and should be a topic of further research.
All results obtained in this work can be reproduced from the repository in Ref. 25 which contains micromagnetic simulation, data analysis, and plotting scripts.

Methods
In our simulations, magnetisation dynamics is governed by the Landau-Lifshitz-Gilbert (LLG) equation 26,27   The first term w ex = A[(∇m x ) 2 + (∇m y ) 2 + (∇m z ) 2 ] with material parameter A in Eq. 4 is the energy density of the symmetric exchange contribution. The unit vector field m, with Cartesian components m x , m y , and m z , represents the normalised magnetisation field M = M s m, where M s is the saturation magnetisation. The second term w dmi = Dm · (∇ × m) with material parameter D is the Dzyaloshinskii-Moriya (DM) energy density. This form of DM energy density is obtained by including Lifshitz invariants for materials of crystallographic classes T and O. The coupling of magnetisation to an external magnetic field H is defined by the Zeeman energy density term w z = −μ 0 M s H · m, with μ 0 being the magnetic constant. The w d term is the demagnetisation (magnetostatic) energy density. Because demagnetisation energy cannot be neglected in the study of FeGe helimagnetic nanostructures 10,14,19,28 we include its contribution in all simulations. We assume the simulated material is isotropic and neglect the magnetocrystalline anisotropy energy contribution. The specific boundary conditions 12 have been validated by a series of standard problem simulations proposed by Cortés-Ortuño et al. 29 . The FeGe material parameters we use are 10  We discretise the finite element mesh so that the maximum spacing between two neighbouring mesh nodes is below l max = 3 nm.
We implemented this micromagnetic model in the finite element method framework and developed a micromagnetic simulation tool Finmag 30 . For the low-level finite element operations, we use the FEniCS project 31 and for the adaptive step time integration we use Sundials/CVODE solver 32,33 .