• Views 495
• Citations 0
• ePub 9
• PDF 286
`Journal of Applied MathematicsVolume 2015, Article ID 317947, 5 pageshttp://dx.doi.org/10.1155/2015/317947`
Research Article

## On a Nonlocal Damping Model in Ferromagnetism

1MIA CNRS EA 3165, Université de La Rochelle, 17000 La Rochelle, France
2FST Errachidia, Laboratoire M2I, Equipe MAMCS, Université My Ismaïl, BP 509, Boutalamine, 52000 Errachidia, Morocco

Received 15 April 2015; Accepted 12 July 2015

Copyright © 2015 M. Moumni and M. Tilioua. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We consider a mathematical model describing nonlocal damping in magnetization dynamics. The model consists of a modified form of the Landau-Lifshitz-Gilbert (LLG) equation for the evolution of the magnetization vector in a rigid ferromagnet. We give a global existence result and characterize the long time behaviour of the obtained solutions. The sensitivity of the model with respect to large and small nonlocal damping parameters is also discussed.

#### 1. Introduction and Preliminary Results

In this work, we are dealing with a mathematical model arising in the theory of magnetization dynamics with nonlocal damping. We will consider the model proposed by Nembach et al. [1]. It is given by a modified LLG equation. The modification lies in the presence of a second-order space derivative of magnetization in the effective field. To describe the model equations, we consider a bounded and regular open set of . The generic point of is denoted by . We assume that a ferromagnetic material occupies the domain .

In what follows, represents the unit sphere of . The magnetization field of the ferromagnetic material which belongs to almost everywhere is denoted by . Its evolution is governed by the following modified LLG equation (see [1, 2]):subject to initial conditionsand boundary conditionswhere the symbol denotes the vector cross product in . The positive constant represents the damping parameter. The effective magnetic field depends on and is given byThe first term on the right-hand side of (4) is called the exchange magnetic field, where the positive constant is the exchange coefficient and the term parameterized by the positive constant describes nonlocal damping in ferromagnets. The -term was introduced in [2]; see also [1]. The third term is the bulk anisotropy field (which generally is taken as linear with respect to ), and is the demagnetizing field, satisfying in the stray field equationwhere is the extension by of outside the magnetic domain .

Notice that, for the sake of simplicity, the bulk uniaxial anisotropy field, generally taken as linear in , and the magnetostatic field which satisfy the stray field equation are not considered since they only induce more computations and have no mathematical influence on the results we obtain.

We define the energyThe following energy estimate holds.

Lemma 1. If is a solution of problem (1), then it satisfies, at least formally, the energy estimate

Proof. The techniques to obtain (7) are analogous to those used, for example, in [3, 4]. We rewrite LLG equation (1) in the following form:Multiplying (8) by and using the saturation constraint almost everywhere yields toIntegrating (9) on , the right-hand side of (9) becomeswhere we have used the summation convention that repeated indices are summed over and denotes the derivative with respect to space variable . Converting the two last terms of the right-hand side of (10) allows us to get (7).
Before dealing with problem (1), let us first review some previous results. We limit ourselves to mentioning a handful of references concerning existence and we refer to the survey [5] for a more detailed bibliographical account. The general framework (although without nonlocal damping, i.e., the case ) has been established in earlier papers (see, for instance, [6, 7]) using Faedo-Galerkin/Penalization (FGP) method. This method gives an approximated sequence of solutions converging to a global solution of the problem. Next results concern systems with further dissipation terms. For example, the modification considered in [8] consists in adding to the standard dissipation term in the LLG equation another higher-order term of the type . The FGP method is also used to solve the problem. In [9], a model with dry-friction dissipation, which is accounted by adding a dry-friction-like term to the standard Gilbert damping, is studied. Let us mention that a model of ferromagnetic material with hysteresis effects is studied in [10]. In this model, the magnetic moment behaviour is described by the nonlinear Landau-Lifshitz equation with an additional term modeling the hysteresis. This term takes the form of a maximal monotone operator acting on the time derivative of the magnetic moment. For this relaxed model, local existence of regular solutions is proved. Note that in the framework of current-induced magnetization dynamics the work [11] addresses global existence of weak solutions to a LLG equation where a transport term is added to effective field taking into account the effect of the injected current. A model of magnetization switching with inertial effects modeled by means of a second-order time derivative term in the effective field is considered in [12] where existence of weak solutions and their long time behaviour were established. All these proofs are based on some penalization and using various kind of regularizations.

To state the existence result, we start with the definition of weak solutions to problem (1).

Definition 2. Given that such that almost everywhere in we call a weak solution to LLG equation (1) if(1)for all , , , and satisfies the saturation constraint for almost everywhere in ;(2) in the trace sense;(3)for all , there holds(4)for all , we havewhere is given by (6).

We have the following global existence result.

Theorem 3. Let be such that almost everywhere. Then there exists a global weak solution of problem (1) in the sense of Definition 2.

Proof. The proof follows a standard scheme (Faedo-Galerkin method with a penalization of the saturation constraint) as usual in general Landau-Lifshitz-Gilbert equation; see [6, 7, 13].

The rest of the paper is divided as follows. In the next section we investigate the long time behaviour of the solutions. Section 3 discusses the sensitivity of the obtained solutions with respect to nonlocal damping factor . More precisely we characterize the limit problem for both high and small . We conclude the paper in Section 4 by giving some comments.

#### 2. The Limit as Goes to

We investigate the long time behaviour of the solutions. More precisely, we study the -limit set of the trajectories and characterize the -limit points as solutions of a suitable stationary problem. We proceed as in Carbou-Fabrie [13].

Let be a weak solution of (1). We call -limit set of the trajectory the following set:

Consider a weak solution of (1). From the energy estimate (12), the -limit set is nonempty. We denote by a point of this set. There exists a nondecreasing sequence such that and in weakly. Since is a smooth bounded domain, then tends to in strongly for ,and extracting a subsequence, we assume that tends to almost everywhere, so that the saturation constraint is satisfied almost everywhere. In addition, we remark that, for all , almost everywhere, so that . By interpolation inequalities in spaces, we obtain that, for all , tends to in strongly.

For and we define for large enoughWe have the following convergence result.

Lemma 4. The sequence satisfies the following convergences:

Proof. Following [14], we have the estimateAs lies in , one gets Since tends to in strongly, tend to in strongly. Moreover, we have obviously seen that is bounded in so there exists a subsequence still denoted by such that tends to in weakly and in strongly almost everywhere in . This ends the proof of the lemma.

Now, we consider a function such that In the weak formulation (11) we take as test function , where is a function lying in and fixed. We obtain after the change of variables We take the limit of the previous equation when tends to . In order to pass to the limit, we bound each term of the above formulation. For example, for the last term we have Since belongs to , this last term tends to zero as goes to . In the same way we pass to the limit in the other terms to obtain which impliesfor all .

We proved the following.

Theorem 5. If is a weak solution of (1), then each point in is a weak solution of the steady state systemwhich should be understood in the weak sense (21).

Remark 6. The solutions of (22) are known as harmonic maps to the unit sphere, which turn up in several equations in physics, such as Ginzburg-Landau equation, and were extensively investigated mathematically in the past decades due to their importance both in mathematics and in many applied fields; see, for example, [15, 16]. For the links between harmonic maps and the solutions of the LLG equation we refer to the book [17, Chapter 4]. Some results on the dynamic version of (22) can be found in [18].

#### 3. Limiting Process for High and Small Nonlocal Damping

We first analyze the limiting process when goes to . We denote by a global weak solution of LLG equation (1) associated with and the initial data and satisfying the energy inequality (7). We have the estimates.

Lemma 7. There exists independent of such that the sequence satisfies the estimates

We also have the following.

Lemma 8. The sequence is compact in .

Lemmas 7 and 8 imply the following convergence results.

Lemma 9. There exists a subsequence still denoted by such thatMoreover, satisfies the saturation condition .

The above convergences allow one to conclude.

Theorem 10. Let be the limit of a subsequence of as goes to . Then, the domain is uniformly magnetized.

Remark 11. The result of Theorem 10 is of interest. In fact, the added term in the effective field can act as a control for magnetization switching.

Our aim now is to pass to the limit as goes to zero in (1). We denote by a global weak solution of LLG equation (1) associated with and the initial data and satisfying the energy inequality (7). The bound on in (see (7)) allows one to get the following.

Theorem 12. Let be the limit of a subsequence of . Then, satisfiesMoreover, satisfies the saturation condition almost everywhere.

Remark 13. From (25), we can say that, for lower values of , the magnetic nonlocal damping vanishes and the classical switching continues for ever. In other words, the magnetization precesses several times around the effective field direction before it reaches equilibrium.

#### 4. Concluding Remarks

In this paper, we have considered nonlocal damping in magnetization dynamics. The model consists of a generalized LLG equation that contains a term characterizing nonlocal damping expressed in terms of in the effective field. The long time behaviour of the solutions is characterized and the sensitivity of the model to nonlocal damping parameter is discussed. The results obtained can be applied without difficulty to the case of effective field with anisotropy and demagnetizing fields. Note that the model considered in this paper neglects the additional damping due to transversal spin currents whose form is . In this case the LLG equation takes the following form (see Nembach et al. [1]):Interestingly, the additional damping cannot be written in terms of the free energy, and therefore it cannot be derived from the functional derivative of the free energy with respect to the local magnetization. This new term can significantly change the domain-wall structure in ferromagnetic materials. It would be interesting to consider this problem from both the theoretical and the numerical points of view. In particular, global existence of weak solutions will require more detailed studies.

We finally mention that an important progress was done to design schemes constructing the weak solutions to the general LLG equation. Several schemes were proposed and their convergence to weak solutions was proved. A significant step forward in the convergence theory of numerical schemes has been done recently; see [1921]. This will be helpful to give a strategy for efficient computer implementation which may reflect the true nature of the augmentation of the LLG model considered in this paper.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was supported by the PHC Volubilis program MA/14/301 “Elaboration et analyse de modèles asymptotiques en micro-magnétisme, magnéto-élasticité et électro-élasticité” with joint financial support from the French Ministry of Foreign Affairs and the Moroccan Ministry of Higher Education and Scientific Research.

#### References

1. H. T. Nembach, J. M. Shaw, C. T. Boone, and T. J. Silva, “Mode- and size-dependent Landau-Lifshitz damping in magnetic nanostructures: evidence for nonlocal damping,” Physical Review Letters, vol. 110, no. 11, Article ID 117201, 2013.
2. Y. Tserkovnyak, E. M. Hankiewicz, and G. Vignale, “Transverse spin diffusion in ferromagnets,” Physical Review B: Condensed Matter and Materials Physics, vol. 79, no. 9, Article ID 094415, 2009.
3. H. Ammari, L. Halpern, and K. Hamdache, “Asymptotic behaviour of thin ferromagnetic films,” Asymptotic Analysis, vol. 24, no. 3-4, pp. 277–294, 2000.
4. K. Hamdache and M. Tilioua, “On the zero thickness limit of thin ferromagnetic films with surface anisotropy energy,” Mathematical Models & Methods in Applied Sciences, vol. 11, no. 8, pp. 1469–1490, 2001.
5. M. Kružík and A. Prohl, “Recent developments in the modeling, analysis, and numerics of ferromagnetism,” The SIAM Review, vol. 48, no. 3, pp. 439–483, 2006.
6. A. Visintin, “On the Landau-Lifshitz equation for ferromagnetism,” Japan Journal of Applied Mathematics, vol. 2, no. 1, pp. 69–84, 1985.
7. F. Alouges and A. Soyeur, “On global weak solutions for Landau-Lifshitz equations: existence and nonuniqueness,” Nonlinear Analysis. Theory, Methods & Applications, vol. 18, no. 11, pp. 1071–1084, 1992.
8. P. Podio-Guidugli and V. Valente, “Existence of global-in-time weak solutions to a modified Gilbert equation,” Nonlinear Analysis: Theory, Methods & Applications, vol. 47, no. 1, pp. 147–158, 2001.
9. T. Roubiček, G. Tomassetti, and C. Zanini, “The Gilbert equation with dry-friction-type damping,” Journal of Mathematical Analysis and Applications, vol. 355, no. 2, pp. 453–468, 2009.
10. G. Carbou, M. A. Efendiev, and P. Fabrie, “Relaxed model for the hysteresis in micromagnetism,” Proceedings of the Royal Society of Edinburgh: Section A Mathematics, vol. 139, no. 4, pp. 759–773, 2009.
11. M. Tilioua, “Current-induced magnetization dynamics. Global existence of weak solutions,” Journal of Mathematical Analysis and Applications, vol. 373, no. 2, pp. 635–642, 2011.
12. M. Hadda and M. Tilioua, “On magnetization dynamics with inertial effects,” Journal of Engineering Mathematics, vol. 88, pp. 197–206, 2014.
13. G. Carbou and P. Fabrie, “Time average in micromagnetism,” Journal of Differential Equations, vol. 147, no. 2, pp. 383–409, 1998.
14. M. Langlais and D. Phillips, “Stabilization of solutions of nonlinear and degenerate evolution equations,” Nonlinear Analysis. Theory, Methods & Applications, vol. 9, no. 4, pp. 321–333, 1985.
15. L. C. Evans, “Partial regularity for stationary harmonic maps into spheres,” Archive for Rational Mechanics and Analysis, vol. 116, no. 2, pp. 101–113, 1991.
16. F. Lin and C. Wang, The Analysis of Harmonic Maps and Their Heat Flows, World Scientific Publishing, Hackensack, NJ, USA, 2008.
17. B. Guo and S. Ding, Landau-Lifshitz Equations, vol. 1 of Frontiers of Research with the Chinese Academy of Sciences, World Scientific, Hackensack, NJ, USA, 2008.
18. P.-L. Sulem, C. Sulem, and C. Bardos, “On the continuous limit for a system of classical spins,” Communications in Mathematical Physics, vol. 107, no. 3, pp. 431–454, 1986.
19. F. Alouges and P. Jaisson, “Convergence of a finite element discretization for the Landau-LIFshitz equations in micromagnetism,” Mathematical Models & Methods in Applied Sciences, vol. 16, no. 2, pp. 299–316, 2006.
20. S. Bartels and A. Prohl, “Convergence of an implicit finite element method for the Landau-LIFshitz-Gilbert equation,” SIAM Journal on Numerical Analysis, vol. 44, no. 4, pp. 1405–1419, 2006.
21. S. Bartels, J. Ko, and A. Prohl, “Numerical analysis of an explicit approximation scheme for the Landau-LIFshitz-Gilbert equation,” Mathematics of Computation, vol. 77, no. 262, pp. 773–788, 2008.