Magnetostatic interaction between Bloch point nanospheres

Three-dimensional topological textures have become a topic of intense interest in recent years. This work uses analytical and numerical calculations to determine the magnetostatic field produced by a Bloch point (BP) singularity confined in a magnetic nanosphere. It is observed that BPs hosted in a nanosphere generate magnetic fields with quadrupolar nature. This finding is interesting because it shows the possibility of obtaining quadrupole magnetic fields with just one magnetic particle, unlike other propositions considering arrays of magnetic elements to generate this kind of field. The obtained magnetostatic field allows us to determine the interaction between two BPs as a function of the relative orientation of their polarities and the distance between them. It is shown that depending on the rotation of one BP related to the other, the magnetostatic interaction varies in strength and character, being attractive or repulsive. The obtained results reveal that the BP interaction has a complex behavior beyond topological charge-mediated interaction.


INTRODUCTION
The analysis of the properties of solitons is a cornerstone in current research regarding technological applications based on the control of magnetic quasiparticles [1][2][3][4][5].The main feature that allows the use of solitonic magnetic patterns as information carriers resides in their strong stability ensured by a topological protection [6].Therefore, the focus of applied research is on describing the fundamental properties of topological textures and their current-driven motion.Up to recent years, the interest was centered in studying solitons lying in quasi-2D-systems, such as skyrmions [7], skyrmioniums [8], biskyrmions [9], and bimerons [10][11][12][13].Nevertheless, recent advances in producing and characterizing nano and microstructures with a plethora of shapes and sizes [14][15][16][17] renewed the interest in describing the properties of three-dimensional (3D) magnetic profiles [18][19][20][21][22][23][24][25][26][27].Among the 3D magnetic quasiparticles with topological protection, we can highlight the Bloch point (BP) [28,29], which is a structure that presents a singular point at its center, where ferromagnetic order is destroyed [30][31][32].The defining property of a BP is that in a closed surface around its center, the direction of the magnetization field covers the whole solid angle an integer number of times.Theoretical and experimental works showed that these magnetic singularities appear in magnetic nanodots with perpendicular magnetic anisotopy [24], in ferromagnets during the process of vortex core reversal [33][34][35], the reversal of skyrmions in confined helimagnetic structures [36], in a bilayer of nanodots of FeGe with different chiralities [37], in modu-lated nanowires with intrinsic Dzialoshinskii-Moriya interaction [38], and in cylindrical magnetic nanowires as the center of a vortex domain wall [39][40][41].
To properly use BPs in technological applications based on their stabilization and motion, it is crucial to analyze the fundamental properties of these structures.Therefore, several interesting phenomena regarding the BPs behavior have been reported.For instance, the analysis of the magnetization resonant modes in modulated nanowires evidenced that the magnetic response to an external magnetic field can establish a strategy for detecting BPs in such systems [42].Also, from analyzing the scattering of magnons by BPs, Elías et al. [43] showed that this system is described by the same solutions of the system electron×magnetic monopole.Therefore, the non-trivial topological structure of the Bloch point manifests in the propagation of spin waves, endowing them with a gauge potential that yields the emergence of spin wave vortices [44,45] as a result of such a scattering.Another exciting property regarding the analogy between BPs and phenomena belonging to other physical contexts is the spontaneous emission of spin waves when a BP domain wall displaces in a cylindrical nanowire.In this case, when the DW velocity is above a threshold value, one can observe a Cherenkov-type DW breakdown phenomenon originating from an interaction between the spontaneously emitting SW and the BP domain wall [46].
Because the structure of a BP depends on the boundary conditions, an interesting problem, from the fundamental point of view, is the analysis of the BP properties as a function of the geometrical parameters of a magnetic system where it is confined.Therefore, due to the BP symmetry, the sphere is a natural geometry where these topological structures can appear as a metastable static configuration.Although a formal solution describing the radial BP represents the maximum energy state, hedgehog BPs appear as intermediate metastable states in spherical nanoparticles before the appearance of the more stable twisted BPs [47].In this regard, if the complete magnetostatic interaction is considered for determining the BP profile, the twist angle of a BP hosted in a sphere is approximately 105 • [48].This swirling effect originated from the magnetostatic energy leads to a non-zero 3D topological charge, which does not occur for radial BPs [49].Based on the above described, this work analyzes the fundamental properties of a magnetic sphere hosting a BP magnetization pattern.Through analytical calculations, we determine the magnetic field generated outside the sphere as a function of the BP helicity.We show that the magnetic field generated by a BP hosted in the sphere resembles that of a quadrupole.The strength of the quadrupolar field is modulated by a factor that depends on the material magnetic parameters, the polarity, and the helicity of the BP.These results allow us to analyze the interaction between two nanospheres having a BP as the magnetic state.In this case, we obtain the BP interaction as a function of the relative orientation between their polarities and the BP-to-BP distance.We conclude that the strength of the BP interaction depends on the rotation angle of one BP related to the other.Additionally, the relative orientation between the BP polarities also determines the character of the magnetostatic interaction, which can be attractive or repulsive.These obtained results reveal the complex behavior of the BP interaction, whose discussions should be beyond topological charge-mediated interaction.

MAGNETIC FIELD OF A BP
This section presents the analysis of the magnetostatic field generated by nanospheres whose dimensions allow the nucleation and stabilization of a BP singularity as a metastable configuration.The performed calculations are obtained in the framework of the micromagnetism approach, where the magnetization, M, can be considered as a continuous function depending on the position inside the magnetic body.The considered system consists of a magnetic sphere with radius R made of a material having saturation magnetization M s .The normalized magnetization field can be written in spherical polar coordinates as M (r)/M s ≡ m = (sin Θ(r) cos Φ(r), sin Θ(r) sin Φ(r), cos Θ(r)).Under this framework, the magnetic profile of a BP configuration can be modeled with the ansatz [47] Θ where θ and φ are the standard polar and azimuthal angles describing the spherical coordinates, p = ±1 stands for the BP polarization, q ∈ Z corresponds to the winding number, and γ is the helicity.Although γ has slight variations along the radial position inside the sphere [49], we consider it a fixed value in the sphere volume.Thus, in this parametrization, the (2D) topological charge is Q = pq.It is worth noticing that this work focuses on determining the properties of topological magnetic textures with q = 1 (BP).The properties of their q = −1 counterparts (Anti BP) will be considered in future works To better describe the magnetic textures determined from Eqs.
(1), we present some types of BPs with different helicities and polarities in Fig. (1).One can notice that while p determines the direction in which the magnetic moments in the poles of the sphere point, γ determines the direction of the magnetic vector field in the sphere equator.The analysis of Fig. 1 evidences that, as expected, the magnetization vector field of a BP obeys the hairy ball theorem [50], which states that any continuous tangent vector field on the sphere must have a vanishing point.
FIG. 1. Magnetization vector field of BP singularities confined in a nanosphere for different values of p and γ.
One of the main objectives of the present work is determining the magnetostatic field H d , generated by a BP nucleated in a ferromagnetic nanosphere in the whole space.Indeed, although the calculation of H in d inside the sphere has already been performed [47,48], the magnetic field, H out d , generated by the BP in the region outside the sphere has not been addressed in detail.Determining the specific profile of H out d is important to better understand the fundamental physics behind BP structures and to analyze the magnetostatic interaction between two nanospheres nucleating BPs as metastable states.
In the absence of current densities, the micromagnetism approach allows us to obtain the magnetic field as where n is the unitary vector normal to the surface, the first and second integrals in the above equation are evaluated along the sphere volume (V ) and external surface area (S), respectively.The magnetostatic potential given in Eq. ( 3) can be fully calculated by expanding the Green's function G(r, r ) = |r − r | −1 into the spherical harmonic basis (see details in the Appendix ).Therefore, after some algebra, we obtain the magnetostatic potential generated by an isolated BP confined in a nanometric sphere, given by and where the superscripts in and out are associated with the magnetostatic potential inside and outside the sphere, respectively.In this context, from the definition of the demagnetizing field, we have that The substitution of Eqs. ( 4) and ( 5) in Eq. ( 6) yields the magnetostatic field inside and outside of a sphere hosting a BP, respectively given by and One can notice that the value of γ plays different roles in the magnetic field inside and outside the sphere.That is, in the region r > R, γ contributes to a global factor to the magnetostatic field, which vanishes when the BP is spherically symmetric (γ = 0 or π for p = 1 or p = −1, respectively).In such a case, the dipolar term in a multipolar expansion is expected to be null.As a consequence, hedgehog-like BPs can be only stabilized in the absence of dipolar energy, since the latter is proportional to ∝ m • H d .On the other hand, because γ does not contribute with a global factor in the magnetostatic field inside the sphere, even spherically symmetric BPs generate a dipolar field inside the nanoparticle.Indeed, a hedgehog BP pointing radially outside (p = 1 and γ = 0) generates a constant magnetostatic field pointing radially inside the sphere.This demagnetizing field yields the magnetic moments to deviate from the radial direction, and adopts a quasi-tangential configuration in the nanosphere equator [47,48].
The vector field of the magnetostatic field outside the sphere is depicted in Fig. 2, for p = 1 and γ = π/2.From the analysis of the presented vector field and Eq. ( 8), one notices the emergence of a magnetic field with quadrupolar nature, whose strength decreases with γ.This result is very interesting since quadrupole magnetic fields can be used, for instance, for trapping atoms in radio-frequency dressing experimental settings [52], designing magnetic tweezers to exert forces on magnetic particles [53], and changing the properties of ferrofluid flow and heat transfer [54].The quadrupolar fields are generally obtained by two ferromagnetic bars parallel to each other, with the north pole of one next to the south of the other.The same field profile can be generated by two properly spaced coils with currents in opposite directions.Another interesting configuration of magnetic elements to obtain a quadrupolar field is considering four pole tips, with two opposing magnetic north poles and two opposing magnetic south poles.At the micro and nanoscale size, this last configuration can be observed in square spin ice systems [55][56][57], and magnetic tweezers [53].Therefore, the magnetic system considered in our work overcomes the necessity of using two or more magnetic elements to generate a quadrupole magnetic field.

INTERACTION BETWEEN TWO MAGNETIC SPHERES HOSTING BLOCH POINT
Once we have established the particular shape of the magnetostatic field generated by a Bloch Point confined in a nanosphere, we can explore how two BP's interact through the magnetostatic potential.Let us consider the interaction between two BP's hosted in spherical nanoparticles with the same radii R 1 = R 2 = R, separated by a distance d.We also consider that the magnetic spheres are not subject to any other force than that produced by the interaction between the magnetization field of one of them and the magnetic field produced by the other.Under these assumptions, the total magnetostatic energy E m of the system is given by where M (r) = M 1 (r) + M 2 (r) and H is the magnetic field produced by the magnetic spheres.The situation can become very complex if the magnetic field is strong enough to modify the structure of the BPs.In this case, finding a new equilibrium BPs configuration would be necessary to minimize the system's energy.Due to this difficulty and the fact that the BP structure is mainly dominated by its self-dipolar field, we assume that the BPs structure do not change.This assumption can be justified by considering hard ferromagnetic materials in such a way that weak external fields cannot modify their configuration.In this approximation, Eq. ( 9) can be simplified to where the subindices 1 and 2 relates with the two BPs nucleated, respectively, in the blue and red spheres represented in Fig. 3.
From the reciprocity theorem [51], we can state that both terms in the right-hand of the previous equation are equal.Therefore, we can rewrite Eq. ( 10) as 3. Schematic representation of the idealized interaction between two nanospheres with a BP as a metastable state.
In a), we assume that the blue BP (sphere 1), localized in the frame O, generates a magnetostatic field that allows the interaction with the red BP (sphere 2), localized in the frame O .The vector d represents the separation between O and O , forming an angle α with the Y -axis.The BP in O is allowed to perform rigid rotations in response to the interaction energy, which are represented in b) and c) for two limiting cases.
In above equation, we consider that the BP stabilized in the sphere denoted by 1 generates a magnetostatic field H 1 in the region where the BP 2 , which has magnetization field 11) reads It is worth noticing that Eq. ( 13) is integrated out in the frame where the BP 2 is located.Therefore, we need to express the potential U d1 , generated by the BP 1 , from the frame of the sphere 2 (see Fig. 3).We have then determined a relationship between the two reference systems by assuming that, in the most general scenario, the two reference frames, O and O , are located in the center of the spheres 1 and 2, respectively.Additionally, we consider that the BPs axis coincides with the Z and Zaxis in each case.The reference frames are separated by a vector d, pointing from O to O , as depicted in Fig. 3. Also, sphere 2 (O ) is free to perform rigid rotations in response to interacting with the other BP.In this context, after some mathematical manipulations and considerations (see Appendix for further details), we obtain the potential U d1 in the frame O , given by where c 1 = M s1 (p 1 −cos γ 1 )R 4  1 /12 and we introduced the rotation angles ξ and ψ, corresponding to two sucessive rotations around X and Y , respectively.
A general result for the integrals presented in Eq. ( 13) becomes quite complicated.However, we can obtain interesting results by exploring some particular cases of interaction.Thus, we consider that the BP 1 (which generates de magnetostatic field) is characterized by γ 1 = arccos (−p 1 /4) (we recall the fact that γ = 0 produces a null interaction energy) and fix its polarity p 1 pointing in the Z−axis direction.Similarly, for the BP 2 , we adopt γ 2 = arccos (−p 2 /4) and consider that p 2 = ±1.It is worth noting that because we assume that the interaction between the BPs does not change their helicities significantly, we use the γ value for a BP free of external interactions [48].Finally, we fix the position of BP 1 and vary the position of BP 2 along the plane yz.In this way, we explore three representative cases for the interaction between the BPs as a function of d = (0, d cos α, d sin α), where α is the angle between d and the axis-Y (see Fig. 3).Finally, we restrict our analysis to cases where the BPs interaction is appreciable, i. e., the separation between them is |d| = R 1 + R 2 + , being a small and positive parameter.The main idea is to describe how the magnetization field of the BP 2 changes in the presence of the magnetostatic potential U d1 .Therefore, we explore possible rigid rotations of the BP structure, which can be interpreted as variations of the direction to which the BP 2 magnetic moments in the sphere poles (Z -axis) point.Thus, from determining the interaction energy as a function of the rotation angles ψ and ξ, we obtain the respective rotation angle minimizes the energy or, equivalently, the direction to which Z -axis points in order to minimize the system's energy.From now on, we numerically solve Eq. ( 13) to determine the system's energy.In the obtained results, we adopt arbitrary unities, that is, R = 1, µ 0 = 1, M s = 1, = 0.1.Therefore, all the presented numerical results give us the qualitative behavior of the analyzed system.sults are depicted in Fig. 4, where we present the magnetostatic energy (in arbitrary unities) as a function of the rotation angles ψ and ξ that determine the final direction to which the Z -axis points.Fig. 4-(a) and (b) depict to the magnetostatic energy when the BP 2 has a polarity p 2 = +1 and p 2 = −1, respectively.An opposite behavior can be observed from the analysis of the obtained results.From Fig. 4-(a) we can see that the minimum energy is degenerated in four pairs of rotation angles {ψ, ξ}, namely {(0, 0), (π, 0), (0, π), (π, π)}.In this case, we can state that the BPs tend to align the magnetic moments in the sphere poles in the same direction ( Ẑ = ± Ẑ).Additionally, the minimum energy is independent of the helicity orientation (clockwise or anticlockwise).Indeed, establishes an energy barrier between the distinct magnetic ground states for the considered system.In contrast, when p 2 = −1, we obtain an opposite behavior.That is, the system minimizes its energy for ξ = π/2, whatever the rotation angle ψ, and the maximum values for the magnetostatic energy degenerate in the four pairs of rotation angles which minimize the energy in the case p = +1.The magnetization profile of the BP 2 in the minimum energy configuration is depicted in Fig. 5-(b).In Fig. 5-(c), we present the magnetization profile of BP 1 just for comparison with the final magnetic states of the BP 2 .
Case 2: d = (0, d, d) We continue our analysis on the magnetostatic interaction between two BPs hosted in nanospheres by considering that they lie in the Y Z-plane, and their relative position forms an angle α = 45 • .Our results evidence the same qualitative behavior as compared with the case of α = 0 • .That is, there are certain rotation angles where the minimum energy is reached, depending on the BP polarity and the specific orientations between the axes Z and Z .Nevertheless, due to the rotational symmetry of the BPs, the configuration that minimizes the energy when α = 45 • is infinitely degenerated, as evidenced from the results presented in Fig. 6, which shows the magnetostatic energy as a function of the rotation angles ψ and ξ.For each rotation ψ, around X-axis, there is an associated rotation ξ, around Y -axis, that leads the BPs to assume the minimum energy configuration.Indeed, the analysis of Fig. 6-(a) reveals that for p 2 = 1, there are infinite pairs of angles {ψ, ξ} that minimize the energy.For instance, the minimum energy configuration of the BP 2 for the case ψ = 0 is shown in Fig. 7-(a).Indeed, only one rotation around the Y -axis (ξ = π/2) is allowed since any other value of ξ would yield higher energy.Additionally, two pairs of rotation angles produce the configuration with maximum energy, e.g., (ψ, ξ) ∼ (0, 0) and (ψ, ξ) ∼ (π, π), which are separated by the energy valley formed by the infinitely degenerated minimum energy states.If the BP 2 polarity is p 2 = −1, one observes a similar behavior as the case where p 2 = 1 and α = 0. Indeed, because the inversion of the polarity leads to the opposite behavior in the magnetostatic interaction between the BPs, two pairs of rotation angles minimize the energy, as evidenced by the results presented in Fig. 6-(b).In this case, the pairs (ψ, ξ) ∼ (0, 0) and (ψ, ξ) ∼ (π, π) represent the configuration with minimum energy, which are separated by the energy barrier for the infinite number of pairs (ψ, ξ).Nevertheless, it is worth noticing that although the points (ψ, ξ) ∼ (π, 0) and (ψ, ξ) ∼ (0, π) do not correspond to the minimum energy configuration, they can also be considered metastable states due to the pronounced energy barrier at the lines ξ ∼ π/2.The possible configurations of BP 2 having minimum energy or metastable states are presented in Fig. 7-(b).Finally, we consider the interaction between two BPs lying in the plane Y Z, with their relative position forming the angle α = 90 • .Under this assumption, the BP 2 is located at the Z-axis, just above the BP 1 .We analyze the behavior of the magnetostatic energy as a function of possible BP 2 rotations around the axis X and Y .The obtained results are presented in Fig 8 -(a), from which one can notice that, for p = +1, the configuration having minimum energy is precisely the same as the case α = 0 • .That is, the minimum energy is degenerated in four pairs of rotation angles {ψ, ξ}, which can be reduced to only two rotations because the minimum energy is independent of the helicity orientation, as represented in Fig. 9-(a) in ξ = π essentially changes the sign of γ.These minimum energy states are separated by an energy barrier whose maximum values are associated to an infinite number of possible rotations giver by (ψ, π/2) and (π/2, ξ), as depicted in Fig. 8-(b).Therefore, if ψ = π/2, the energy is independent of the rotation angle ξ.The same behavior is valid if we fix ξ = π/2 and change the values of ψ.Two examples of the BP 2 configuration ((0, π/2) and (π/2, 0)) producing the minimuum energy state are presented in Fig. 9-(b).

Magnetostatic energy as a function of the distance
The above-analyzed cases present the same qualitative behavior regarding the magnetostatic interaction as a function of the rotations of BP 2 concerning BP 1 .From this perspective, we study the nature of these interactions if they are attractive or repulsive.We have then determined the magnetostatic energy of the system composed of the two BPs as a function of the distance between them.To clarify the role of the relative orientation between the BPs on the interaction, we present the results for the magnetostatic energy when the BPs are oriented to reach the minimum and maximum energy state for the cases α = 0 • with p 2 = +1 (Fig. 10-(a)), and α = 90 • with p 2 = −1 (Fig. 10-(b)).Although we present only these two particular cases, the same qualitative behavior is observed in all analyzed cases under the respective combinations of α and p 2 = ±1, yielding maximum and minimum energy configurations.Our results show that the magnetic configurations presenting minimum and maximum energies have attractive and repulsive natures.Indeed, the analysis of Figs.10-(a) and 10-(b) reveals that the energy increases with the nanospheres distance when the relative orientation of the BPs leads to the minimum energy configuration, evidencing an attractive interaction.On the other hand, a decrease in the magnetostatic energy is observed when the BPs are oriented at an angle defining the maximum energy configuration, revealing, in this case, the repulsive nature of the BPs interaction.
Because the relative rotation between the BPs changes the nature of the interaction from attractive to repulsive, there should be a set of specific rotation angles for which the BPs do not interact.To corroborate this statement, we calculate the magnetostatic energy numerically as a function of ξ for a fixed ψ and distance between the nanospheres.The obtained results are depicted in Figs.10-(c) and 10-(d) respectively for α = 0 • and α = 90 • , where we fix ψ = 0 • .In specific, for the case α = 0 • with p 2 = +1 the system becomes noninteracting when ξ ∼ 41 • or ξ ∼ 138 • .Similarly, in the case α = 90 • with p 2 = −1, the BPs do not interact when ξ ∼ 54 • or ξ ∼ 125 • .As expected, the qualitative behavior for both cases is similar, but the specific angles for which the interaction vanishes depend on α.

DISCUSSION AND CONCLUSIONS
The knowledge of the magnetostatic field generated by magnetic nanoparticles is an important issue to consider when analyzing the interaction between magnetic particles.In this work, we analyzed the magnetic field of nanospheres hosting a BP as a metastable state.Our results evidenced the quadrupolar nature of the magnetostatic field generated by a BP.Furthermore, the strength of this field depends on the BP helicity, vanishing for a hedgehog BP pattern.This result is interesting from both perspectives, fundamental and applied physics.Indeed, the magnetic system considered in our work overcomes the necessity of using two or more magnetic elements to generate a quadrupole magnetic field, which can give a new breath for technological propositions regarding quadrupole magnetic fields [52][53][54].
The magnetostatic interaction between two nanospheres with a BP metastable state is also determined.Our results evidenced that the interaction between BPs depends on the relative position of one BP regarding the other.Additionally, the relative orientation between the BP axis determines the strength of the magnetostatic energy.Depending on the BPs polarities, there are relative orientations of the BPs magnetic moments for which the energy is minimized.Nevertheless, the energy minima are separated by relative rotation angles for which the energy is a maximum.In this context, although the interaction between planar vortices (2D counterparts of BPs) can be analyzed from the perspective of their topological charges [58], we can conclude that the interaction between BPs should consider arguments beyond the approximation of topological charge-mediated interaction.Indeed, BPs with the same topological charge (p 1 = p 2 = +1) can interact with attractive or repulsive potential, depending on the relative rotation angle between them.The same behavior is valid for BPs with opposite topological charges (p 1 = −p 2 ).Finally, because the nature of the interaction can change from attractive to repulsive, we showed that there are relative orientations between the BPs for which they do not interact.These results evidence that nanospheres with BPs as magnetic states should present an exciting interaction dynamic, highly depending on the initial conditions.Although the results presented in this work concern the magnetostatic interaction of BPs, they give some insights in possible ways to understand the complex behavior of three-dimensional topological textures.

Fig. 5 -FIG. 6 .FIG. 7 .
FIG. 6. Magnetostatic energy of a system composed by two BPs interacting along an axis defined by α = 45 • as a function of the rotation angles ψ and ξ.

FIG. 8 .
FIG. 8. Magnetostatic energy of a system composed by two BPs interacting along an axis defined by α = 90 • as a function of the rotation angles ψ and ξ

FIG. 9 .
FIG. 9. Final configuration of the BP2 for the different rotation angles that minimize the interaction energy in the case α = 90 • .Panels (a) and (b) correspond respectively to the cases where p2 = 1 and p2 = −1.For comparison, the configuration of BP1 was included in (c).

FIG. 10 .
FIG. 10.Magnetostatic energy as a function of the BP distance.Panel (a) presents the results for α = 0 • and p2 = +1, while panel (b) depicts the case α = 90 • and p2 = −1.In both cases, open symbols correspond to the energy when the BPs are oriented to reach the maximum energy state, while filled symbols correspond to the energy when the BPs are oriented to reach the minimum energy configuration.Panel (c) depicts the behavior of the energy as a function of ξ for a fix ψ = 0 in the case α = 0 • with p2 = +1.Panel (d) shows the corresponding energy as a function of ξ for ψ = 0 in the case α = 90 • and p2 = −1.