Nanometric alternating magnetic field generator

In this work we introduce an alternating magnetic field generator in a cylindrical nanostructure. This field appears due to the rotation of a magnetic domain wall located at some position, generating a magnetic region that varies its direction of magnetization alternately, thus inducing an alternating magnetic flux in its vicinity. This phenomenon occurs due to the competition between a spin-polarized current and a magnetic field, which allows to control both the angular velocity and the pinning position of the domain wall. As proof of concept, we study the particular case of a diameter-modulated nanowire with a spin-polarized current along its axis and the demagnetizing field produced by its modulation. This inhomogeneous field allows one to control the angular velocity of the domain wall as a function of its position along the nanowire allowing frequencies in the GHz range to be achieved. This generator could be used in telecommunications for devices in the range of radiofrequencies or, following Faraday’s induction law, could also induce an electromotive force and be used as a movable alternate voltage source in future nanodevices.

field, and a rotation in the plane perpendicular to the main axis of the nanostructure, which periodically varies the magnetization of the nanostructure, generating an alternating magnetic field 27 . Both motions occur simultaneously while the nanostructure continues under the influence of a magnetic field or current (see Fig. 1b).
Considering the work previously mentioned, in this work we present a time-stable nanometric alternating magnetic field generator, which is capable of maintaining the rotation of the DW in a fixed position, thus permitting the in-situ knowledge of the source of field variation. In this way, we will establish the relation that must be satisfied by the spin-polarized currents and the magnetic fields to obtain the phenomenon described previously. As proof of concept, we will show that by means of a spin-polarized current and a demagnetizing field, created by the geometric change of a modulated diameter nanowire, it is possible to fix the position of a TDW, and generate a stable rotation of this DW.

Analytical Model
Initially, we will consider a nanowire with cylindrical symmetry containing a TDW. We will also consider that the DW is simultaneously under the action of a magnetic field and a spin-polarized current. In this way, the dynamics of the TDW can be estimated by the Landau-Liftshitz-Gilbert equation, with the inclusion of additional STT terms 28,29 : where M is the local magnetization, H eff the effective field, M s the saturation magnetization, γ the gyromagnetic ratio, α the Gilbert damping factor, and β the non-adiabatic spin-transfer parameter. The vector u is defined as where g is the Landé factor, μ B the Bohr magneton, e the electron charge, J the current pin , the TDW stops its linear movement, however it keeps rotating with an angular velocity φ  pin in an equilibrium position z pin .
Scientific RepoRts | 7: 4736 | DOI:10.1038/s41598-017-05026-4 density and P the polarization rate of the current (u is positive for P > 0). For simplicity, we will use ρ = u J , where ρ is a constant of proportionality.
As a starting point, we will use the equations proposed by Yan et al. 27 , who adapted the analytic model proposed by Mougin et al. 25 . These equations describe the velocity of a single 180° DW in a magnetic nanostructure of reduced dimensions with cylindrical symmetry under the action of an external applied magnetic field and/or an electrical current. The linear velocity of the DW is proportional to θ  , while the angular velocity of rotation is φ  . When the TDW is driven just by a static field applied along the z direction (H = H z , u = 0), the velocities are given by θ αφ , and the velocities of the magnetization due to the current-driven case where wc denotes the center of the TDW (θ = π/2), θ(z) is the profile of the wall and u z is the vector in the z-direction defined previously. Following from these ideas, we explicitly write the equations that consider the combined action of an external magnetic field and a spin-polarized current, θ γα α αβ ρ In these equations, Δ corresponds to the characteristic wall width, while u z = ρJ z . From the above equations it can be deduced that a TDW can remain rotating in a fixed position if the condition of θ =  0 occurs. Then, to stop the propagation of the TDW, it is required that the magnetic field and the electron flow have opposite directions, satisfying the following relation (see Fig. 1c): Under these circumstances, combining Eqs (2) and (3), we obtain the angular velocity, φ  J ( ) pin z , when the TDW remains pinned, pin z z 2 Figure 2 shows the magnetic field H z pin (obtained from equation 3) and the angular velocity φ  pin (obtained from equation 4) as a function of the current density, with values typically used in experiments. From this figure we can observe both the magnetic field and the angular velocity of rotation increase linearly with the current density. In addition, it is interesting to note that this angular velocity of the pinned TDW reaches GHz, which allows this configuration to have applications in the area of communications.

Micromagnetic Simulations
As proof of concept, we will consider a demagnetizing field created by the geometric change of a cylindrical nanowire with modulated diameter, by means of which it is possible to fix the position of a TDW, and this position generates a rotation of the DW stable in time. For this, we performed numerical simulations using a micromagnetic simulator 30 extended to consider also the dynamics of the magnetization induced by a spin-polarized current 12 . Our starting point is a cylindrical diameter modulated Ni 80 Fe 20 nanowire with a total length L = 2030 nm. The nanowire is composed of three segments: the segment A has a diameter d A = 10 nm and a length L A = 1000 nm, the segment B has a diameter d B = 40 nm and a length L B = 1000 nm, and between both segments, the diameter decreases linearly from 40 to 10 nm over a 30 nm length, as shown in Fig. 3a.
We have used a saturation magnetization of M S = 800 × 10 3 A/m, a stiffness constant A = 13 × 10 −12 J/m and cell sizes of 1.25 × 1.25 × 2.5 nm 3 . Besides we have neglected the crystalline anisotropies since the nanowires are polycrystalline. Additionally, the current density varies with the transverse area of the nanostructure and, from now on, when we talk about simulations with electric current, we will refer to the current density of segment A. Finally, we use α = 0.01, β = 0.04, and P = 0.4 26,29 . Figure 3b shows the initial configuration of the system, where a TDW lies near the free end of segment A, and due to this it minimizes energy versus other types of DWs 3 . Then, at t = 0 s, an electric current is applied with current densities J ranging between 1 × 10 12 A/m 2 and 1 × 10 13 A/m 2 , as typically used in experiments 19,[31][32][33] . It is important to consider that these high current densities can only be applied in pulses which typically last less than 100 ns. In our case, and as we will show later (see Fig. 4a), the angular velocity will generally converge for times less than 40 ns, so that if a 100 ns pulse is applied, the system will exhibit a constant angular velocity for about 60 ns. The current densities used in our simulations produce an Oersted field that can be neglected 29 , and has therefore not been considered in our simulations. An electric current along −z displaces the TDW towards the +z direction, i.e., in the electron flow direction (from segment A to segment B).
The advantage of using a diameter-modulated nanowire is to replace the source of the external magnetic field by the demagnetizing field created naturally by the modulation and segment B. As can be seen in Fig. 3c, the demagnetizing field H d in the thin nanowire region is inhomogeneous and it is oriented in the −z direction, which could allow one to slow the propagation of the TDW with direction +z, driven by an electric current. Due to the non-homogeneity of the demagnetizing field along the wire for certain values of current density, it is possible to find an equilibrium position z cr satisfying Eq. (3) but not suppressing its rotation, where = H z H ( ) d cr z pin . From our simulations, we can obtain that at z = 1000 nm, the demagnetizing field H d reaches values close to 0.5 T in the −z direction, and decreases rapidly along the segment A. Thus, when the TDW is far from the modulation, the effect of the demagnetizing field on the DW is negligible. However, when the TDW is near the interface, the influence of the demagnetizing field becomes important and both θ  and φ  velocities are affected.
From the results obtained one may observe a linear relationship between the speed of propagation of the TDW in the segment A and the current density. Once the TDW approaches and reaches the modulation, the DW strongly slows its propagation, increasing its angular velocity and two possible results are generated. First, if the current is sufficiently high, J > 2 × 10 12 A/m 2 , the TDW will continue its movement towards and along the segment B, significantly reducing its velocity, which will be proportional to the decrease in current density due to the change in diameter from d A = 10 nm to d B = 40 nm. On the other hand, for low current densities, J ≤ 2 × 10 12 A/m, the TDW will stop its displacement at a point near the modulation due to the contribution of the demagnetizing field when it is equal to that of the spin-polarized current, thus finding the equilibrium position z cr satisfying Eq. (3). Figure 4a shows the behavior of the angular velocity φ  pin of the TDW as a function of time for the regime of a low current density. We find that for J = 1 × 10 12 A/m 2 , the TDW stops its motion at z cr = 996.5 nm, but it continues rotating to reach a maximum angular velocity of 48.2 rad/ns. Then, the angular velocity varies periodically over time until it converges to a stable value of 37.2 rad/ns, which is equivalent to a frequency of 5.92 GHz, in the microwave range. At this point it is important to mention that we have performed micromagnetic simulations for different values of the non-adiabatic spin-transfer parameter β (not shown here), and the frequency variation is less than 3%. From Eq. (4), and considering a current density of J = 1 × 10 12 A/m 2 , we have obtained an angular velocity φ = .
 36 6rad/ns pin that is in good agreement with the value obtained from the micromagnetic simulations.
Moreover, we find that for J = 2 × 10 12 A/m 2 the TDW stops its motion into the modulation at z cr = 1002 nm, but it continues rotating to reach an average angular velocity of 61 rad/ns that is equivalent to a frequency of 9.78 GHz, with periodic variations between 59 rad/ns and 63 rad/ns. This frequency variation is due to the overlap between the TDW and the spin waves generated when the TDW enters the modulated zone. For higher currents densities a more intense magnetic field is necessary to stop the TDW. However, the modulation zone generates a magnetic field whose maximum value is −0.5 T, which is not enough to stop the TDW when it moves due to the application of current densities J > 2 × 10 12 A/m 2 .
The behavior described above not only refers to the geometric parameters chosen for the simulated nanostructure, but it will in principle work for any nanostructure that satisfies equation 4. In fact, if a conductive circuit is placed near the nanostructure, then an electromotive force will be induced because of the change in the magnetic flux created by the rotation of the TDW, as shown in Fig. 4b. In this way, we have proposed this nanostructure as a nanometric alternating magnetic field generator of controlled frequency and position 34 .
In the literature there are several studies that show that the width of the domain walls increases while the magnetization decreases with an increase in temperature [35][36][37] . The model used in this article considers that there is a magnetic field that stops the displacement of a domain wall. In this case, the magnetic field is assumed to be approximately constant within the domain wall. However, if the temperature increases considerably the domain wall width, probably the magnetic field will no longer be constant inside the wall. Furthermore, increasing the temperature decreases the magnetization of the nanowire, which implies that the demagnetizing field will decrease and probably will not be able to stop the domain wall.
On the other hand, there are numerous articles [38][39][40][41][42] which indicate that a current density in the order of 10 12 A/ m 2 generates a temperature of about 0-100 K for each ns of pulse over a nanowire. Of course, the temperature increase will depend on the geometric and magnetic parameters, the magnitude and time of the pulse, and the substrate where the nanowire will be supported. For all of the above, our model is valid for low temperatures, and therefore, in order to operate the device presented in this article, suitable conditions must be sought so that the applied current density does not considerably increase the temperature of the device. A good thermal conductor as a substrate is an alternative to achieve a lower average temperature. For example, Fangohr et al. 39 showed that a Permalloy nanowire with notch placed on a silicon wafer exhibits a maximum temperature increase of 17 K when a 20 ns pulse (current density at the wire ends is 10 12 A/m 2 ) is applied, conditions very similar to those exhibited by our nanogenerator. In addition, Parkin et al. 19 showed that by choosing a suitable substrate it is possible to move magnetic domain walls using pulses between 20 and 100 ns with a current density in the order of 10 12 A/m 2 , same order of magnitude used in our investigation.

Conclusions
In conclusion, we have demonstrated through analytical calculations and micromagnetic simulations that the competition between a spin-polarized current and a magnetic field allows to stop the translation movement of a TDW traveling along a cylindrical geometry, but without stopping its angular movement, generating a magnetic region that varies its direction of magnetization alternately, thus inducing an alternating magnetic flux. From a theoretical point of view, we have obtained analytical expressions for both the magnetic field H z pin necessary to stop the translation movement of the DW (Eq. (3)) and its angular velocity φ  pin (Eq. (4)). Interestingly, both expressions increase linearly with the current density. Furthermore, as proof of concept, we have investigated the propagation of a TDW driven by a spin-polarized current in a cylindrical diameter modulated Ni 80 Fe 20 nanowire by means of micromagnetic simulations. It is observed that once the TDW approaches and reaches the region of change of diameter, the DW strongly decelerates its propagation, completely stopping its displacement for certain current densities. Finally, we obtained that the angular velocity of the pinned TDW is in the microwave range, which allows these structures to have applications in communications.