Spinmotive force due to motion of magnetic bubble arrays driven by magnetic field gradient.

Interaction between local magnetization and conduction electrons is responsible for a variety of phenomena in magnetic materials. It has been recently shown that spin current and associated electric voltage can be induced by magnetization that depends on both time and space. This effect, called spinmotive force, provides for a powerful tool for exploring the dynamics and the nature of magnetic textures, as well as a new source for electromotive force. Here we theoretically demonstrate the generation of electric voltages in magnetic bubble array systems subjected to a magnetic field gradient. It is shown by deriving expressions for the electric voltages that the present system offers a direct measure of phenomenological parameter β that describes non-adiabaticity in the current induced magnetization dynamics. This spinmotive force opens a door for new types of spintronic devices that exploit the field-gradient.

Spinmotive force (SMF) refers to the generation of spin current, which is accompanied by an electric voltage, as a result of dynamical magnetic textures in conducting ferromagnets123. This is due to the exchange coupling between conduction electrons and the local magnetization. SMF reflects the temporal- and spatial-dependence of the local magnetization456789, and thus it offers a powerful method to probe and explore the dynamics and the nature of magnetic textures. In addition, SMF can be a new source of electromotive force, directly converting the magnetic energy into the electric energy of conduction electrons. While the classical electromagnetism tells us that the conventional inductive electromotive force requires a time-varying magnetic flux, it has been reported that an electromotive force can be generated by a static and uniform magnetic field via the SMF mechanism101112.

In the last few years, more attention has been focusing on topologically nontrivial magnetic structures such as magnetic vortices in soft ferromagnetic nanodiscs1314 and skyrmion lattices in chiral magnetic thin films1516. The SMF offers some insights into and gains benefit from such magnetic systems; the polarity of a moving magnetic vortex core can be electrically detected1718, and arbitrarily-large ac SMF was predicted due to skyrmion lattice motion19. To the best of our knowledge, however, there has been no work on SMF induced in systems that contain magnetic bubble domains. Magnetic bubbles are observed in ferromagnetic films with out-of-plane anisotropy as spot-like closed domains, where the magnetization is oriented to the opposite direction to the one outside the bubbles. The structures of magnetic bubbles are similar to those of skyrmions and vortices in the sense that bubbles carry a topological number called skyrmion number. Since their first observation in the 1960s, magnetic bubbles have been showing distinctive and interesting behaviors20212223.

In this paper, SMF due to the motion of magnetic bubble arrays is theoretically investigated based on the steady-motion model. As the bubble motion may be induced by a spatially-varying magnetic field, our work reveals that a magnetic field gradient can be exploited to generate a spin current and an associated electric voltage. By deriving expressions for the electric voltages, we demonstrate their cumulative nature, i.e., they can be proportional to the number of involved bubbles. An important fact to be stressed is that the present system can confirm the SMF originating from the non-adiabatic dynamics of the electron spin, leading to a direct measurement of the controversial, so-called, β parameter.

Results

Steady-motion model for magnetic bubble

Let us begin by reviewing the collective-coordinate model for magnetic bubble dynamics. We consider a cylindrical bubble domain with radius R in a thin film. This type of magnetic structure can be stabilized by applying bias magnetic field in the out-of-plane with appropriate magnitude. The distribution of the magnetization direction m = (sin θ cos ψ, sin θ sin ψ, cos θ) is assumed to be two dimensional and described as20 [see Fig. 1]where (r, χ, z) is the cylindrical coordinate measured from the bubble center, Δ is the wall width parameter, ψ0 is a constant, Q is the topological parameter defined asand S is the winding number:where is the contour integral taken counterclockwise around the circumference. The domain wall separating the bubble and the outside, in general, can contain vertical Bloch lines, i.e., there are many possibilities in the way of distributing the azimuthal angle ψ along the perimeter. In Eq. (2) we assumed the linear dependence of ψ on χ, as we focus on the following two cases. First, when the magnetization rotates one full turn around the wall of the bubble with no Bloch line, ψ = χ ± π/2, that is, S = 1 and ψ0 = ± π/2, with the sign + (−) corresponding to left (right) handed chirality. Eq. (2) is also a good approximation when the Bloch lines are packed so closely that the distance between the adjacent Bloch lines is comparable to the wall width, i.e., . If one considers a small number of Bloch lines, the distribution of ψ would be no longer as simple as Eq. (2).

Figure 1

Schematic of magnetic bubble structures in a thin film.

Red and Blue arrows in the upper figure and black ones in the bottom indicate the magnetization.

When a magnetic field is applied in the z direction with its magnitude varying in the x-y plane, a bubble is driven to move in the film seeking positions with lower Zeeman energy. In the following the steady motion of the bubble under the constant gradient ∇H is assumed, i.e., during its motion with constant velocity v the bubble stays rigidly cylindrical with constant radius R and the ψ-distribution does not change with respect to the coordinate frame moving with the bubble:

Assuming Eq. (5) and that the magnetization dynamics obeys Landau-Lifshitz-Gilbert equation, the equation of motion for the bubble is given by20 (the derivation is shown in Methods)where γ is the gyromagnetic ratio and α is the dimensionless Gilbert damping parameter. represents the unit vector along the z-direction. An assumption that was made when deriving Eq. (6) is that the spin transfer with the conduction electrons are negligible. Remarkably, the bubble is deflected away from the direction of the field gradient at an angle that is determined by the material parameters (see Methods for the detail). The extension of the above discussion to multiple-bubble problem is straightforward.

Spinmotive force due to bubble motion

Let us examine the SMF induced by the steady motion of bubble indicated by Eq. (5). We assume that a conduction electron in the ferromagnetic film is described by a one-body Hamiltonianwhere m is the electron's mass. The second term represents the exchange interaction between the electron spin and the magnetization, with Jex being the exchange coupling energy. According to theory of SMF456, dynamical magnetization exerts an effective electric field on the electrons via the exchange coupling, which is called spin electric field since its sign depends on the electron spin (see Methods):with

The upper (lower) signs in Eq. (8) correspond to the electron with majority (minority) spin. and are referred to as adiabatic and non-adiabatic spin electric fields, respectively, as β = ħ/2Jexτsf is the dimensionless parameter describing the non-adiabaticity in the electron spin dynamics5689, with τsf the relaxation time for the electron spin flip. The spin electric fields (9) and (10) require both time and spatial dependences of the magnetization, and this condition is satisfied around the perimeter of moving bubbles.

The spin electric field (8) induces a spin current and a charge current in the sample, with the electric conductivity for the majority (minority) electrons. These currents generate by-products such as the charge redistribution, the spin accumulation, and the charge/spin diffusion current. In an open circuit system, the electric field ind = −∇ϕ−∂/∂t appears to keep the total charge current zero, i.e., , where ϕ and are the electromagnetic scalar and vector potentials. An electric voltage between two given points and , which is given by the difference in the electric potential ϕ() − ϕ(), enables one to detect the spinmotive force electrically. To determine the gauge potentials one has to fix the gauge, and here let us adopt the Coulomb gauge, ∇ · = 0. From the above equation for the open circuit condition, one obtains the Poisson equationHere is the conservative electric field induced by the electric potential distribution, where is the spin polarization of the conduction electrons. In the above argument, we have neglected the contribution from the diffusive current to the total charge current as we have metals in our mind as the samples; once the effective U(1) electric field is given, the problem of computing the electric voltage induced by the electric field falls within the classical electromagnetism and the established transport theory, and it is known that in metals the induced diffusion potential is mostly negligible compared to the electric potential, unlike in semiconductors. Technically, the diffusive current can be taken into account by replacing the electric potential −eϕ by the electrochemical potential , where is the Fermi energy. ( may be dependent on the space and the spin electric field in complex ways.)

Eq. (11) can be applied to systems with arbitrary sample geometry and magnetic texture. In Fig. 2, we show an example of electric potential distribution by numerically solving the Poisson equation (11) with spin electric fields (9) and (10), where the steady motion of eight identical bubbles in a square-shape thin film is assumed [see Methods for the numerics]. Here the coordinate system is set that the bubble array flows along the x direction. It is seen that the potential drop occurs at the position of the bubbles, as is expected. Notice that the adiabatic field gives rise to the net potential drops only in the y direction (perpendicular to the bubble flow), while the non-adiabatic one only to the x direction (along the bubble flow), indicating that in this setup the two contributions can be separately identified by longitudinal and perpendicular voltage measurements. While S = 1 is assumed in Fig. 2, qualitatively much the same profiles are also obtained in the case of |S| = R/Δ, but with height of each potential drop being different (not shown). The S dependence of the electric voltage will be discussed later.

Figure 2

The distributions of electric potential ϕ induced in a thin film due to (a) the adiabatic field (9) and (b) the non-adiabatic field (10), calculated by solving the Poisson equation (11) numerically, where the steady motion of eight identical bubbles along the x direction is assumed.

The profile of each bubble is given by Eqs. (1) and (2), with R = 50 nm, Δ = 2 nm, Q = 1, S = 1 and ψ0 = π/2. The other parameters assumed here are γ = 1.76 × 1011 T−1 s−1, α = 0.02, , P = 0.5, and the side lengths are 900 nm. The field gradient makes an angle 75.9° to the x axis with its strength |R∇H| = 10 Oe.

Next, let us adopt a “quasi-one-dimensional” approximation for the electric voltage and derive its analytic expression. We limit ourselves to a rectangular thin film where bubbles move along either of two sides of sample, as is the case in Fig. 2. Assuming a bubble moving in the positive x direction, we estimate the electric voltage V along the motion of the bubble by where L( is the side length of the sample along the x(y) direction. The y-integral operation and the appearance of the factor 1/L are for spatial-averaging of the electric potential along the y-axis. Eqs. (1), (2), (5), (6) and (8) have been used to carry out the integrals. Eq. (12) gives the exact solution for V in one-dimensional systems: . In two-dimensional systems as is the present case, Eq. (12) is still a good approximation when effects on V from ∂ϕ/∂y and the detailed profile of ∂ϕ/∂x can be ignored. This is often the case if there are some appropriate symmetries in the system, and the divergence of the electric field, i.e., the moving bubble, keeps away from the electrodes. Similar concept for the electric voltage V measured in the perpendicular direction to the motion of the bubble leads to

The dc electric voltages appear both in the x and y directions, being proportional to the field gradient for both magnetic configurations S = 1 and . It is seen from Eq. (12) that V for S = 1 is larger than that for |S| = R/Δ under the same applied field, because R > Δ. With the parameters shown in Fig. 2, V for S = 1 is about an order of magnitude greater compared to that for |S| = R/Δ. On the other hand, V has little dependence on |S|.

When there are multiple bubbles, the net electric potential distribution is given by the super-position of all individual electric potentials produced by each bubble. If N identical bubbles move in the x direction, the simplest extensions, and , of the above expressions for the electric voltages, V and V, respectively, may be

Figure 3 compares Eqs. (14) with the electric voltages obtained by numerically solving the Poisson equation (11), showing good agreement between them.

Figure 3

The field dependence of dc electric voltages induced by the steady motion of arrays of identical bubbles.

Three different configurations are examined, which are indicated in the top panel. The dynamics of each bubble is described by Eqs. (1), (2) and (5) with S = 1 [see the caption of Fig. 2 for the other parameters assumed here]. (a) obtained by Eq. (14) and by numerically solving the Poisson equation (11) correspond to solid lines and open symbols, respectively, and (b) similarly for . The electric voltages are proportional to the field gradient, and depends on the number of bubbles and the sample geometry. See also Discussion in the main text.

Discussion

Eqs. (12), (13) and (14) indicate key features of this SMF. is proportional to N/L, and thus, roughly speaking, depends on the “density” of bubbles along the y axis and the “number” of bubbles along the x axis, which contribute to, respectively, the height of each potential drop and the number of occurrence of potential drop along the x axis. The dependence of on the number of bubbles and the sample geometry is demonstrated in Fig. 3(a). Comparing the configurations i and ii, which share the same sample shape, the slope of V is twice larger for ii because the configuration ii contains twice as many bubbles, i.e., the sites where the potential drop occurs. It is not necessarily that the same value of N leads to the same magnitude of electric voltage; the configuration iii provides twice larger than ii does under the same applied field because of the difference in the factor 1/L. Similar discussion is applied for , see Fig. 3(b). Eq. (14) indicates that one may control the dc electric voltages by adjusting the sample geometry and the number of bubbles.

Remarkably, only the non-adiabatic field contributes to , while only the adiabatic one to . Since in most magnetic materials β is believed to be smaller than unity and thus , it is difficult to identify the contribution of the non-adiabatic field in the conventional systems. In fact, there has been no experimental confirmation of its effects. Now, by the measurement of one can detect free from the larger adiabatic contribution. This can lead to direct and unambiguous measurement of the phenomenological parameter β; as both V and V contain in their expressions P, which may be another uncertain material parameter, one can be free from P by measuring the ratio of V and V.

In the present study, we assumed the specific and simple profile of magnetization, Eqs. (1), (2) and (5), to derive the analytic expressions for the electric voltages. Although the investigation of the effects of the bubble shape distortion, local disorder, and etc. must be interesting, it will require systematic study based on numerical approaches to the magnetization dynamics, which is out of reach of this paper. The assumption made in the present study is reasonable when the gradient of applied field is sufficiently moderate22, and our analytical model will work well in this field range.

In conclusion, We have shown for the first time that a magnetic field gradient can generate a SMF, i.e., spin currents and associated electric voltages, by driving the motion of magnetic bubbles. Based on the steady-motion model, expressions for the dc electric voltages in longitudinal and perpendicular to the bubble motion are derived, which turned out to be controllable by tuning the sample geometry and the number of involved bubbles. An important implication of our result is that the present system offers an experimental determination of the phenomenological parameter β that describes non-adiabaticity in the electron spin dynamics. This SMF can lead to a new route for basic study of the electron-magnetization interaction as well as a new concept in spintronic devices, exploiting the gradient of magnetic fields.

Methods

Steady-motion of bubble

In order to get to the equation of motion for the bubble, consider the increment of the stored magnetic energy , with the magnetic energy density w, due to the variations δθ and δψ:where h is the film thickness and MS is the saturation magnetization. δw/δθ and δw/δψ have been expressed in terms of ∂θ/∂t and ∂ψ/∂t by using Landau-Lifshitz-Gilbert equation of motion without spin-transfer-torque effect, and the time derivatives have been further expressed in terms of r and χ assuming Eqs. (1), (2) and (5). The change δU in the internal energy is supposed to be balanced by the external pressure on the bubble due to ∇H:where denotes the position of the center of the bubble. Rearranging the above equation, Eq. (6) is obtained. The net mobility of the bubble is found by solving Eq. (6) for || in terms of |R∇H| as

The angle ρ of deflection of the bubble away from the field gradient may be defined asWith S = 1, Q = 1, α = 0.02, Δ = 2 nm and R = 50 nm, corresponding to the calculation in Fig. 2, one obtains .

Derivation of spin electric fields

Under the Hamiltonian (7), the Heisenberg equation of motion for the electron is given bywith and denoting the operators for the electron's position and the force acting on the electron, respectively. The actual motion of the electron is obtained by determining the expectation value 〈〉↑↓ of the electron spin with majority (↑) and minority (↓) states. Notice that · ∇ = 0 and thus the component of 〈〉↑↓ that is (anti-)parallel to does not contribute to the force. Let us decompose the electron spin as , where the upper (lower) sign corresponds to the majority (minority) spin and δ↑↓ represents a slight deviation from . The expectation value of the force is written as ; what causes non-zero force due to the exchange coupling is a misalignment between the electron spin and the magnetization. Assume that the electron spin dynamics is described by

The first term on the right-hand side represents the Larmor precession about the magnetization, and the damping motion toward the magnetization is phenomenologically introduced by the second term, which describes the non-adiabaticity in the electron spin dynamics, with τsf the relaxation time for the electron spin flip. By substituting the above expression for 〈〉↑↓ into Eq. (20), δ↑↓ is expressed in terms of , by which one obtains with given by Eq. (8).

We have considered an open circuit condition. More generally Eq. (20) should include the divergence of the spin current carried by conduction electrons, leading to appearance of the spin magnetic field. While we have followed Ref. 9 here, essentially the same result was obtained by a different approach where Onsager's reciprocal relation is taken into account56.

Numerical approach to the Poisson equation

In the numerical calculations, the sample is divided into number of meshes, and a single magnetization vector is assigned to each mesh, where i is the index of the meshes. The time evolution of {(t)} (i = 1, 2, …, Nm), where Nm is the number of the meshes, is given based on Eqs. (1), (2), and (5). As the spin electric field is updated at each mesh by Eq. (8), the induced electric potential distribution {ϕ(t)} satisfies the Poisson equation [see Eq. (11)]

In two-dimensional discrete systems, the above equation is equivalent towhere (i, j) are the indices for the meshes in two dimension, and ΔxΔy is the area of the mesh. The potential distribution is obtained by solving the above equation self-consistently with Neumann boundary condition, where the spatial derivative of the electric potential is zero at the sample edge.

Author Contributions

Y.Y. carried out theory calculation and wrote the manuscript. S.H. analyzed the data. J.I. and S.M. conceived the problem. J.I., S.M. and J.S. supervised the study. All authors discussed the results and commented on the manuscript.