Micromagnetic Simulations of Submicron Vortex Structures for the Detection of Superparamagnetic Labels

We present a numerical investigation on the detection of superparamagnetic labels using a giant magnetoresistance (GMR) vortex structure. For this purpose, the Landau–Lifshitz–Gilbert equation was solved numerically applying an external z-field for the activation of the superparamagnetic label. Initially, the free layer’s magnetization change due to the stray field of the label is simulated. The electric response of the GMR sensor is calculated by applying a self-consistent spin-diffusion model to the precomputed magnetization configurations. It is shown that the soft-magnetic free layer reacts on the stray field of the label by shifting the magnetic vortex orthogonally to the shift direction of the label. As a consequence, the electric potential of the GMR sensor changes significantly for label shifts parallel or antiparallel to the pinning of the fixed layer. Depending on the label size and its distance to the sensor, the GMR sensor responds, changing the electric potential from 26.6 mV to 28.3 mV.


Introduction
Magnetoresistive sensors are part of many technical applications and are already implemented in several industrial sectors. Due to their high sensitivity and integration capability, magnetoresistive sensors are also a promising technology in the field of flexible electronics [1], human-computer interactions [2] and biomedicine [3]. In the case of biomedical applications, the detection of magnetic signals can be challenging since human organs generate very small magnetic fields ranging from nT to fT [4]. Therefore, magnetoresistive sensors can be combined with bio-functionalized labels to enhance magnetic interaction. Excited by an external field, these labels generate magnetic stray fields, which can be detected electronically by a magnetoresistive sensor element [5].
In contrast with other technologies, this provides a multitude of advantages. Compared to the detection of fluorescence signals induced by light absorption of biological samples, no interference appears with most biological samples when using magnetic stray fields [5]. Additionally, the strong coupling of the magnetic stray field in magnetically labelled samples with an external magnet offers the opportunity for efficient label separation and transportation [6]. Since magnetic field sensors are already implemented in commercial motion-control applications, CMOS techniques have been Magnetic vortices arise in micron-sized, nanometer-thick, circular ferromagnetic disks, where the lateral extension is much larger than the thickness of the layer. In the absence of a magnetic field, the magnetization curls up in-plane along the edges to minimize magnetostatic energy, leading to a flux closure state with an out of plane magnetization in the center of the disk (vortex core) [12]. In the presence of an external magnetic field, the parallel and antiparallel oriented magnetic domains increase and decrease equally, and the magnetic vortex moves proportionally to the external magnetic field until it annihilates at the annihilation field Ha (Figure 1b). This change in the layer's magnetization results in a change in electrical resistance if the vortex layer represents the free layer of a GMR spin-valve stack [10].
In principle, ferromagnetic and superparamagnetic labels can generate detectable magnetic fields. As the size of several important biomarkers is assigned to the nanometer range, the size of magnetic labels to be detected should be comparable [13]. In this range, most labels show superparamagnetic properties reaching critical diameters smaller than 100 nm [14]. In contrast to superparamagnetic labels, ferromagnetic labels exhibit a nonvanishing magnetization, even in the absence of a magnetic field. As a consequence, ferromagnetic labels are attracted to each other and Magnetic vortices arise in micron-sized, nanometer-thick, circular ferromagnetic disks, where the lateral extension is much larger than the thickness of the layer. In the absence of a magnetic field, the magnetization curls up in-plane along the edges to minimize magnetostatic energy, leading to a flux closure state with an out of plane magnetization in the center of the disk (vortex core) [12]. In the presence of an external magnetic field, the parallel and antiparallel oriented magnetic domains increase and decrease equally, and the magnetic vortex moves proportionally to the external magnetic field until it annihilates at the annihilation field H a (Figure 1b). This change in the layer's magnetization results in a change in electrical resistance if the vortex layer represents the free layer of a GMR spin-valve stack [10].
In principle, ferromagnetic and superparamagnetic labels can generate detectable magnetic fields. As the size of several important biomarkers is assigned to the nanometer range, the size of magnetic labels to be detected should be comparable [13]. In this range, most labels show superparamagnetic properties reaching critical diameters smaller than 100 nm [14]. In contrast to superparamagnetic labels, ferromagnetic labels exhibit a nonvanishing magnetization, even in the absence of a magnetic field. As a consequence, ferromagnetic labels are attracted to each other and tend to aggregate. This limits the performance of a sensor and leads to potential health issues due to embolism if an in vivo sensing concept is contemplated [13]. By using superparamagnetic labels, this can be avoided. Superparamagnetic properties arise in small ferromagnetic particles below a critical size that depends on the material properties. In this case, labels will only provide a magnetic moment if a small external field is applied. In the absence of the external field and above the blocking temperature, the spins of the label flip randomly, and the averaged magnetic moment disappears [14]. Hence, the magnetic moment of a superparamagnetic label can be activated by the time the label reaches the detector using an external field (magnetizing field) perpendicular to the sensor surface [15].
In developing a sensor system for the detection of a single label, micromagnetic simulations help understand the interaction of the sensor and the label. Here, the 3D finite element magnum.fe tool with the self-consistent feature as published in [16] is used to introduce a spherical label directly to the simulation domain and vary its diameter, distance, horizontal and vertical position to the surface of a simplified GMR element. By solving the Landau-Lifshitz-Gilbert equation numerically, the changes in the sensor's magnetization generated by the stray field of the label are determined and used to simulate the GMR response in a post-processing procedure via the self-consistent micromagnetic model [16].

Simulation Methods and Parameters
Basically, a GMR spin-valve system is fabricated as an array including several hundreds of single GMR elements. In a basic setup, each element consists of a magnetic free layer and a magnetic fixed layer separated by a non-magnetic conducting layer ( Figure 1a). The fixed layer is usually designed to have a stable and homogeneous magnetization configuration, which is typically achieved by an antiferromagnetic exchange bias layer [17]. On the other hand, the free layer magnetization is sensitive to changes in the external field as introduced by the magnetic labels. The relative change in fixed and free layer magnetization is reflected in its GMR signal.

Simulation Model
In order to model the free layer's magnetization rotation due to the label's stray field, the Landau-Lifshitz-Gilbert equation describes the magnetization dynamics of a three-dimensional magnetic domain [18]: It connects the precessional damped motion of the normalized magnetization m to an effective magnetic field h eff , where γ and α are the gyromagnetic ratio and the damping constant, respectively. According to the problem setting, the effective field term h eff includes the demagnetization field h demag and the exchange field h ex describing the dipole-dipole interaction and the quantum-mechanical effect of the exchange interaction. In order to simulate the activation of the superparamagnetic label, the effective field is completed by a Zeeman-field contribution h zeeman [19].
For the numerical solution of the Landau-Lifshitz-Gilbert equation, the 3D finite element method tool magnum.fe is used [20]. The terms of Equation (2) are defined as reported in [20] as well as the corresponding boundary conditions. During the generation of the finite element mesh, the label is placed close to the surface of the soft-magnetic free layer to determine how the layer's magnetization is affected by the stray field of the label. Subsequently, this free-layer configuration is used to calculate the GMR response of the sensor by solving the self-consistent spin-diffusion model.
Within the spin-diffusion approach, the spin accumulation s is introduced characterizing the excess of one spin orientation due to the magnetization direction of a ferromagnetic material. Its diffusion through the layer stack is defined by the equation of motion in [16], which provokes the spin current j s . The spin accumulation s is assumed to be always in equilibrium, since it undergoes temporal changes two orders of magnitude smaller compared to the dynamics of the magnetization, which leads to Equation (3) [16].
Sensors 2020, 20, 5819 4 of 13 The spin current j s is defined by Equation (4) coupled to the charge current j e in Equation (5) via the electric potential u using E = −∇u.
As reported in [16], the system of Equations (3)-(5) is solved in a self-consistent fashion by introducing boundary conditions for the electric potential u, the spin accumulation s and the electric current j e . For this purpose, Equation (3) is complemented by the source equation for the electric current j e [16].
∇ · j e = 0 (6) By using this model, the electric potential u of a GMR element can be directly extracted from the simulation characterizing the GMR effect for an applied electric current. The respective material parameters of Equations (3)-(5) are specified in the following section.
In general, superparamagnetic properties can be observed below a particle diameter of 100 nm, where a magnetizing field h zeeman of a few mT provokes a strong saturation magnetization M s,p in the range of 20 emu/g to 90 emu/g [30,31]. In the simulation, the label's saturation magnetization M s,p is assumed to be 312 kA/m, which corresponds to 60 emu/g using a material density of 5.2 g/cm 3 [32]. Following [30,31], a magnetizing z-field µ 0 h zeeman of 100 mT is applied, and the label's used exchange constant A ex,p of 10 pJ/m is extracted from [33].
Both the diameter of the label (d p ) as well as its lateral and vertical position relative to the sensor's free layer are varied. While d p receives values between 20 nm and 100 nm, the vertical position ∆z is varied between 10 nm and 100 nm ( Figure 2a). The horizontal and lateral position is defined as the label shift in the x-or y-direction (∆x/∆y) starting from the center of the circular free layer ( Figure 2b). It is changed in the negative and positive x-and y-direction from 25 nm to 50 nm and 75 nm.
According to [34], the detection of a single magnetized label is possible as long as the size of the sensor structure is similar to the size of the label. Therefore, a sensor with a diameter in the order of 100 nm is needed to detect the specified magnetite labels. Typically, the free layer of a GMR element consists of a soft ferromagnetic material [35]. The formation of vortices in such submicron layers was theoretically and experimentally investigated in [36][37][38][39], showing that the formation of a vortex state highly depends on the aspect ratio of the layer. Ensuring a stable vortex appearance, this was complemented by Equation (3), which connects the layer thickness t L to the layer diameter d L and the exchange length of the material (l ex,L = Using a CoFe alloy (M s,L = 2.0 T [40], A ex,L = 15 pJ/m [10]) as the free layer material, an exchange length l ex of 3.1 nm is obtained. With the objective of d L ≈ 100 nm, t L has to be > 8.5 nm, and a layer thickness of t L = 10 nm will lead to a stable vortex structure. This is confirmed by the simulation result displayed in Figure 3, where the vortex state occurs as the CoFe layers ground state after 6 ns. According to [34], the detection of a single magnetized label is possible as long as the size of the sensor structure is similar to the size of the label. Therefore, a sensor with a diameter in the order of 100 nm is needed to detect the specified magnetite labels. Typically, the free layer of a GMR element consists of a soft ferromagnetic material [35]. The formation of vortices in such submicron layers was theoretically and experimentally investigated in [36][37][38][39], showing that the formation of a vortex state highly depends on the aspect ratio of the layer. Ensuring a stable vortex appearance, this was complemented by Equation (3), which connects the layer thickness tL to the layer diameter dL and the exchange length of the material (lex,L = √2Aex,L/μ0Ms,L²) [40].
Using a CoFe alloy (Ms,L = 2.0 T [40], Aex,L = 15 pJ/m [10]) as the free layer material, an exchange length lex of 3.1 nm is obtained. With the objective of dL ≈ 100 nm, tL has to be > 8.5 nm, and a layer thickness of tL = 10 nm will lead to a stable vortex structure. This is confirmed by the simulation result displayed in Figure 3, where the vortex state occurs as the CoFe layers ground state after 6 ns.  Considering the exchange lengths of the label and the free layer (lex,p = 12.8 nm and lex,L = 3.1 nm), a tetrahedron mesh is generated using an element size of 3 nm. The damping constants αL and αp are set to 1.0, since only the steady state of the free layer's magnetization will be used to compute the GMR effect within the self-consistent spin-diffusion model.
The single GMR element consists of the CoFe free layer [41], a 2 nm thick Cu layer and a 5 nm  Considering the exchange lengths of the label and the free layer (l ex,p = 12.8 nm and l ex,L = 3.1 nm), a tetrahedron mesh is generated using an element size of 3 nm. The damping constants α L and α p are set to 1.0, since only the steady state of the free layer's magnetization will be used to compute the GMR effect within the self-consistent spin-diffusion model.
The single GMR element consists of the CoFe free layer [41], a 2 nm thick Cu layer and a 5 nm thick CoFeB fixed layer [35,42]. In order to satisfy Equations (5)- (7), further material parameters are required, namely, the diffusion constant D 0 , the conductivity C 0 , the spin-flip relaxation time τ sf , the coupling strength J and the polarization parameters β and β' [43]. Except C 0 , these material parameters are taken from [16], while the conductivities of CoFe, CoFeB and Cu are assumed to be 2.5 MA/Vm [44], 0.3 MA/Vm [45] and 30 MA/Vm [46], respectively. Lastly, the fixed layer's magnetization is pinned in the x-direction and an electric current of I = 8 mA is applied in the z-direction.

Change in the Free Layer's Magnetization
During the simulation, the magnetization direction of the label becomes parallel to the external z-field (Figure 4b). If the label is located in the center above the free layer, the effective stray field of the label acting on the free layer's x/y plane will be balanced in all spatial directions, and the magnetic vortex will not move. Therefore, no significant magnetization change will occur in the x-or y-direction.

Change in the Free Layer's Magnetization
During the simulation, the magnetization direction of the label becomes parallel to the external z-field (Figure 4b). If the label is located in the center above the free layer, the effective stray field of the label acting on the free layer's x/y plane will be balanced in all spatial directions, and the magnetic vortex will not move. Therefore, no significant magnetization change will occur in the x-or ydirection. Due to the magnetizing external z-field, the z-component of the free layer increases compared to the setup without a label from mz/Ms,L = 0.01 to mz/Ms,L = 0.1. Additionally, the z-component is enhanced by the label's stray field in the z-direction, which is strongest for large labels at short distances (Figure 4a).
By shifting the label in the x-or y-direction, the magnetization of the free layer changes significantly within the x/y plane. For example, an x-shift (∆x = 25 nm) of a label with a diameter of dp = 80 nm placed in a height of ∆z = 20 nm leads to an increase in the relative averaged xmagnetization from 0.00 to 0.08. Accordingly, the rising spin alignment in the x-direction leads to an orthogonal movement of the vortex core, which slides in the positive y-direction, as illustrated in Figure 5. This can be observed vice versa for a label shift in the y-direction, since the magnetic stray field of the fixed layer is negligible due to the use of a synthetic antiferromagnet. Due to the magnetizing external z-field, the z-component of the free layer increases compared to the setup without a label from m z /M s,L = 0.01 to m z /M s,L = 0.1. Additionally, the z-component is enhanced by the label's stray field in the z-direction, which is strongest for large labels at short distances (Figure 4a).
By shifting the label in the x-or y-direction, the magnetization of the free layer changes significantly within the x/y plane. For example, an x-shift (∆x = 25 nm) of a label with a diameter of d p = 80 nm placed in a height of ∆z = 20 nm leads to an increase in the relative averaged x-magnetization from 0.00 to 0.08. Accordingly, the rising spin alignment in the x-direction leads to an orthogonal movement of the vortex core, which slides in the positive y-direction, as illustrated in Figure 5. This can be observed vice versa for a label shift in the y-direction, since the magnetic stray field of the fixed layer is negligible due to the use of a synthetic antiferromagnet.
The orthogonal counterclockwise movement of the magnetic vortex core compared to the label's shift direction is initiated by the spins surrounding the magnetic vortex core. At the edge of the core, the magnetization vector consists of magnetization components in the z-and x-/y-direction, Sensors 2020, 20, 5819 7 of 13 depending on the respective magnetic domain. Without shifting the label, its stray field interacts with the out-of-plane vectors around the vortex core with perfect symmetry. Hence, the magnetic vortex remains in the center of the free layer.
to the setup without a label from mz/Ms,L = 0.01 to mz/Ms,L = 0.1. Additionally, the z-component is enhanced by the label's stray field in the z-direction, which is strongest for large labels at short distances (Figure 4a).
By shifting the label in the x-or y-direction, the magnetization of the free layer changes significantly within the x/y plane. For example, an x-shift (∆x = 25 nm) of a label with a diameter of dp = 80 nm placed in a height of ∆z = 20 nm leads to an increase in the relative averaged xmagnetization from 0.00 to 0.08. Accordingly, the rising spin alignment in the x-direction leads to an orthogonal movement of the vortex core, which slides in the positive y-direction, as illustrated in Figure 5. This can be observed vice versa for a label shift in the y-direction, since the magnetic stray field of the fixed layer is negligible due to the use of a synthetic antiferromagnet. Figure 5. Sketch of the orthogonal movement of the magnetic vortex core with respect to the shift in the label. The shift in the magnetic vortex core leads to resistance changes across the GMR sensor stack, reaching its maximum resistance for label shifts in the x-direction. By shifting the label in the x-direction (Figure 6), the magnetization in the label center tilts in the x-direction. This is initiated by the stray field of the vortex and the pull-back mechanism of the z-component in the vortex center. The resulting x-component of the label bottom side enhances the x-magnetization of the free layer surface, similar to the increase in the z-component in center position (see Figure 4). This enhancement naturally shifts the vortex core in the y-direction, increasing the free layer's net magnetization in the x-direction. The orthogonal counterclockwise movement of the magnetic vortex core compared to the label's shift direction is initiated by the spins surrounding the magnetic vortex core. At the edge of the core, the magnetization vector consists of magnetization components in the z-and x-/y-direction, depending on the respective magnetic domain. Without shifting the label, its stray field interacts with the out-of-plane vectors around the vortex core with perfect symmetry. Hence, the magnetic vortex remains in the center of the free layer.
By shifting the label in the x-direction (Figure 6), the magnetization in the label center tilts in the x-direction. This is initiated by the stray field of the vortex and the pull-back mechanism of the zcomponent in the vortex center. The resulting x-component of the label bottom side enhances the xmagnetization of the free layer surface, similar to the increase in the z-component in center position (see Figure 4). This enhancement naturally shifts the vortex core in the y-direction, increasing the free layer's net magnetization in the x-direction. The amount of the magnetization change within the free layer can be manipulated by changing the diameter of the label and its distance to the layer surface, which is shown in Figure 7. Large particles and short distances generate relative magnetization changes up to 0.16 in the x-and ydirection, respectively.
It is noticed that the local magnetization change for a certain particle diameter dp depends on the label's shift length and the direction of the shift. This transfer curve shows properties of a sinusoidal function (Figure 7), where a label shift of 0.5dL (50 nm) provokes the highest interaction. At this position, the label's stray field urges the spins strongly to an alignment in the direction of the label's shift. Towards the center, this force is reduced, and the label's stray field acts more and more in all spatial directions on the free layer. The amount of the magnetization change within the free layer can be manipulated by changing the diameter of the label and its distance to the layer surface, which is shown in Figure 7. Large particles and short distances generate relative magnetization changes up to 0.16 in the x-and y-direction, respectively.
It is noticed that the local magnetization change for a certain particle diameter d p depends on the label's shift length and the direction of the shift. This transfer curve shows properties of a sinusoidal Sensors 2020, 20, 5819 8 of 13 function (Figure 7), where a label shift of 0.5d L (50 nm) provokes the highest interaction. At this position, the label's stray field urges the spins strongly to an alignment in the direction of the label's shift. Towards the center, this force is reduced, and the label's stray field acts more and more in all spatial directions on the free layer.
At the edges, this force is also reduced since the label slowly disappears from a position above the free layer's surface. Thus, the magnetization change decreases towards the edges.
Sensors 2020, 20, x FOR PEER REVIEW 9 of 14 Figure 7. Change in averaged x-magnetization due to the x-shift in the magnetized label. The y-shift in the label leads to an equal distribution of the averaged y-magnetization. The inset shows the label diameters dp and heights ∆z used for the simulation.

Resistance Changes due to the Label's Stray Field
In this section, the GMR from the magnetization patterns in the presence of the label is calculated following the procedure outlined in Section 2.2. Figures 8 and 9 show the results of the potential change for the respective label shifts. According to the findings of the magnetic simulations, a label shift in the ±x-direction leads to a shift in the vortex core in the ±y-direction. Due to more parallel or antiparallel spin alignment between the fixed and the free layer, the electric potential is reduced or enhanced by 0.8 mV (at position 0.5dL) compared to the setup without or centered label. The amount of the potential change depends on particle size and distance (Figure 8). For the displacement in the ±y-direction (Figure 9), the distribution of the potential response is similar, but the amount of the potential change decreases significantly, as the fixed layer is aligned in the x-direction. At position 0.5dL, the electric potential is reduced or enhanced by 0.1 mV compared to the setup without a label. Due to the displacement in the ±y-direction, more spins align in the ±y-direction, and the vortex core slides in the ±x-direction. As a consequence, more spins align orthogonally to the pinning direction of the fixed layer, and the GMR is affected only weakly.

Resistance Changes due to the Label's Stray Field
In this section, the GMR from the magnetization patterns in the presence of the label is calculated following the procedure outlined in Section 2.2. Figures 8 and 9 show the results of the potential change for the respective label shifts. According to the findings of the magnetic simulations, a label shift in the ±x-direction leads to a shift in the vortex core in the ±y-direction. Due to more parallel or antiparallel spin alignment between the fixed and the free layer, the electric potential is reduced or enhanced by 0.8 mV (at position 0.5d L ) compared to the setup without or centered label. The amount of the potential change depends on particle size and distance (Figure 8). For the displacement in the ±y-direction (Figure 9), the distribution of the potential response is similar, but the amount of the potential change decreases significantly, as the fixed layer is aligned in the x-direction. At position 0.5d L , the electric potential is reduced or enhanced by 0.1 mV compared to the setup without a label. Due to the displacement in the ±y-direction, more spins align in the ±y-direction, and the vortex core slides in the ±x-direction. As a consequence, more spins align orthogonally to the pinning direction of the fixed layer, and the GMR is affected only weakly.
Without label displacement in the x-or y-direction, only the z-component of the free layer (and thus the magnetization of the small vortex core) is changed (Figure 4a). This magnetization change in a small region compared to the total surface of the free layer does not affect the GMR. This is confirmed by calculations and illustrated in Figure 10, where the center of the x/y plane shows an electric potential equal to the potential of the no label setup (u = 27.42 mV). Furthermore, it is obvious that the largest sensitivity is achieved by label shifts parallel to the magnetization of the fixed layer, in our case the x-coordinate. By using a sensor array with neighboring sensor elements, including orthogonally aligned fixed layers, similar sensitivity in the y-direction can be obtained.
±y-direction (Figure 9), the distribution of the potential response is similar, but the amount of the potential change decreases significantly, as the fixed layer is aligned in the x-direction. At position 0.5dL, the electric potential is reduced or enhanced by 0.1 mV compared to the setup without a label. Due to the displacement in the ±y-direction, more spins align in the ±y-direction, and the vortex core slides in the ±x-direction. As a consequence, more spins align orthogonally to the pinning direction of the fixed layer, and the GMR is affected only weakly.  Without label displacement in the x-or y-direction, only the z-component of the free layer (and thus the magnetization of the small vortex core) is changed (Figure 4a). This magnetization change in a small region compared to the total surface of the free layer does not affect the GMR. This is confirmed by calculations and illustrated in Figure 10, where the center of the x/y plane shows an electric potential equal to the potential of the no label setup (u = 27.42 mV). Furthermore, it is obvious that the largest sensitivity is achieved by label shifts parallel to the magnetization of the fixed layer, in our case the x-coordinate. By using a sensor array with neighboring sensor elements, including orthogonally aligned fixed layers, similar sensitivity in the y-direction can be obtained. Without label displacement in the x-or y-direction, only the z-component of the free layer (and thus the magnetization of the small vortex core) is changed (Figure 4a). This magnetization change in a small region compared to the total surface of the free layer does not affect the GMR. This is confirmed by calculations and illustrated in Figure 10, where the center of the x/y plane shows an electric potential equal to the potential of the no label setup (u = 27.42 mV). Furthermore, it is obvious that the largest sensitivity is achieved by label shifts parallel to the magnetization of the fixed layer, in our case the x-coordinate. By using a sensor array with neighboring sensor elements, including orthogonally aligned fixed layers, similar sensitivity in the y-direction can be obtained.

Discussion
It is apparent that distinct GMR potential changes arise provided that the particle diameter is similar to the diameter of the sensor's free layer. This was already theoretically reported in [34], where also a minimum distance of the label to the sensor surface of 0.2d p /2 was required. We met this criterion in the simulations as well, starting at a minimum height of 10 nm for the largest label with a diameter of 100 nm. Furthermore, it is pointed out that the x-displacement parallel to the fixed layer magnetization leads to larger potential changes compared to the displacement in the y-direction. Figure 8 shows a potential difference up to 0.8 mV between the setup without a label and the setup where the label is placed in the 0.5d L position. This simulated change in the sensor's response is within the range of analytical models and experimental work. As reported in [47], a resistance change of 20 mΩ has been generated by a spin-valve stripe detecting a single magnetite label with a diameter of 16 nm. In addition, [48,49] investigated the detection of multiple labels and showed that resistance and potential changes of 100 mΩ and 90 mV arise using GMR spin-valve stripes with a size of a few micrometers. According to [47], the transfer curve of the vortex GMR element (Figure 8) differs from the transfer curve of a GMR spin-valve stripe. For the circular vortex structure, the position of the highest potential drop or enhancement could be found at 0.5 d L in the x-direction, whereas the center of the spin valve stripe was the position of highest interaction.
Noise is a limiting factor regarding the performance of the sensor system, since it blurs the sensor response. As reported in [10], the phase noise vanishes by the use of magnetic vortex structures. Therefore, only the thermal noise and the 1/f-noise contribute to the total noise of the GMR element [11]. Since the detection of a biomolecule is an application in the low frequency range, the influence of the 1/f noise becomes dominant, and the noise increases for small sensors [50]. For an estimate of the sensor noise in our case, the model reported in [11] can be consulted. Assuming a linear decay for the normalized Hooge parameter with an increasing aspect ratio of the free layer [11], the noise parameter S v can be approximated to S v ≈ 22 nV/

√
Hz using a frequency of f = 10 Hz. For this purpose, a connection of ten elements in series was assumed. Comparing this approximation with the results of [11] shows that S v increases by a factor 4. The square root of the S v 2 integration yields to a noise amplitude of u ≈ 104 nV, compared to an electric potential for the no label setup of u ≈ 27.4 mV, and a response just below 1 mV, which indicates an acceptable relation between the noise and the response of the simulated sensor structure.

Conclusions
In conclusion, we have used micromagnetic simulations to study the response of a GMR vortex structure to the stray field of a magnetized label.
If the label is centered above the vortex, only the average out-of-plane magnetization of the free layer will change from m z /M s,L = 0.01 to m z /M s,L = 0.1 depending on the label distance and the label diameter. For label shifts away from the sensor center, the orientation of the free layer's average stray field is affected by the shift in the label. If the label that is magnetized in the z-direction will be moved in the x-direction, the average stray field produced by the label in the region of the free layer points in the x-direction. Hence, the x-directed magnetic domain in the vortex increases, which in turn leads to the movement of the vortex core in the y-direction. Based on this mechanism, the vortex core therefore always shifts in the orthogonal direction to the label shift.
In our specific setup, a parallel or antiparallel shift in the label along the pinning axis of the fixed layer, results in a modulation of the electric potential between 26.6 mV up to 28.3 mV. This depends on the distance and the diameter of the label, reaching a maximum for small distances (10 nm to 100 nm) and label diameters similar to the sensor diameter (100 nm).
The amount of the potential change is in agreement with already reported results of analytical modeling and experimental work [47][48][49], which was performed for spin valve stripes.
In summary, a submicron GMR element using a vortex free layer has the potential as a sensor for nanoscale label detection with direct electronic readout and low noise operation in lab-on-a-chip