Advances in Mathematical Physics
Volume 2018 (2018), Article ID 5370349, 15 pages
https://doi.org/10.1155/2018/5370349

Research Article

Decay Rate and Energy Gap for the Singularity Solution of the Inhomogeneous Landau-Lifshitz Equation on

Penghong Zhong

Department of Mathematics, Guangdong University of Education, Guangzhou 510640, China

Correspondence should be addressed to Penghong Zhong; zhonghong737@163.com

Received 8 March 2018; Revised 1 July 2018; Accepted 9 July 2018; Published 1 August 2018

Academic Editor: Andrei D. Mironov

Copyright © 2018 Penghong Zhong. 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.

The singularity solution for the inhomogeneous Landau-Lifshitz (ILL) equation without damping term in -dimensional space was investigated. The implicit singularity solution was obtained for the case where the target space is on . This solution can be classified into four types that cover the global and local solutions. An estimation of the energy density of one of these types indicates its exact decay rate, which allows a global solution with finite initial energy under . Analysis of the four aperiodic solutions indicates that energy gaps that are first contributions to the literature of ILL will occur for particular coefficient settings, and these are shown graphically.

1. Introduction

The nonlinear ferromagnetic chain model (FCM) has attracted the attention of physicists and mathematicians. For physicists, FCM represents the possibility of describing magnetic field density evolution with various external fields, while mathematicians study the various solutions or excitations in geometrically motive nonlinear models. Quantitative models depicting FCM were proposed by Landau and Lifshitz in 1935 [1], and the Landau-Lifshitz equation (LLE) was proposed to the dynamics of the nonequilibrium magnetism system:

where is the variational sum of magnetic anisotropy energy and magnetic field energy, denotes magnetization intensity, is the cross product of Euclidean 3-space , and and are constants. If and , then LLE is the harmonic map heat flow (HMF). Otherwise, when , the Gilbert damping vanishes, and when and , the LLE degenerates into the Schrödinger map heat flow (SMF), which is an important equation of differential geometry. SMF can be regarded as a nonlinear Schrödinger equation that contains a derivative term. Although the existence, uniqueness, and the blow-up problem of some nonlinear Schrödinger equations [2–4] are clear, the theorem of SMF becomes more complicated and some further work still needs to be done. Similarly, comparing HMF and some general harmonic system (or even biharmonic equation) [5–8], the mapping system is more complicated than the nonmapping system due to the curvature flow of the Riemannian manifolds.

For clarity, we set and as the complex structure and metric on the Kähler manifold, respectively. Different manifolds and their covariant derivatives are and , respectively, and we define SMF as

Thus, (2) defines a mapping , which can degenerate into other subcases; for example, when and , (2) is the isotropic LLE without Gilbert damping, and if the manifold and complex structure are and (), respectively, then (2) is the hyperbolic isotropic LLE without Gilbert damping.

If and , then (2) can transform into a Schrödinger-type equation using the Hashimoto transform. Under dimensional cylindrical coordinates, let , and the curvature and torsion are, respectively,

Then, employing the Hashimoto transform,

and (2) transforms to

If , (5) is the standard cubic Schrödinger equation for which the Lax pair is obtained, and explicit solutions can be constructed using the Bäcklund transformation. Hence, (5) can be regarded as a simple case of (2). Similar to the inhomogeneous Schrödinger equation, (2) (or (1)) can be extended to the inhomogeneous expression. If we add the inhomogeneous term to (2), SME can be generalized into the inhomogeneous SME (ILL) proposed by Balakrishnan [9]:

where the inhomogeneous term, , is a scalar function. Under radially symmetrical coordinates, (6) can be expressed as

Using the Hashimoto transform (4), (7) transforms into a nonlinear inhomogeneous Schrödinger equation:

Daniel et al. [10] identified the integrable model (8) by analyzing the singularity structure of its solutions and also discussed Lax pairs, Bäcklund transformation, and soliton-like solutions of (7) (or (8)). They showed that inhomogeneous terms play an important role, and Painlevé analysis indicated that (7) (or (8)) was integrable in arbitrary dimensions if and only if .

As all know, the bifurcation structures of the general integrable systems and autonomous differential systems [11–14] are clear due to the properties of Hamiltonian systems [15–18]. However, the structure for the unintegrable system is not so clear. If takes other forms, what would be the solution of (7)? The regularity of the unintegrable case has been somewhat clarified [19, 20] in regard of the singularity behavior of the equation. Two different finite time blowup solutions of (7) were constructed in [20], one being an explicit form and the other being an implicit solution. Careful estimation of the two spatial dimension implicit solutions indicates that the energy density decay rate is . However, if the spatial dimension is any integer , does a global smooth (or blowup) solution of (7) exist? Does (7) also present similar regularity aspects under specific smooth initial data? The regularity of dimensional ILL remains open.

Many physical phenomena develop singular behavior, for example, boundary layers or blowup solutions. Liu [21] analyzed the concentration set of the stationary weak solutions to LLE for the ferromagnetic spin chain. Based on the Ginzburg-Landau approximation, Wang [22] established the existence of a global weak solution for LLE, , with respect to smooth initial boundary data, which is smooth from a closed set with locally finite dimensional parabolic Hausdorff measure. In 2008, Huh [23] constructed infinite energy explicit blowup solutions for the modified LLE. Ding [24] constructed an infinite energy blowup solution for LLE on a hyperbolic target.

The finite energy blowup solution was firstly studied by perturbation methods [25, 26], and it was found that the collapse of a symmetric case with large initial energy formed a singularity where the blowup rate could be estimated by scale invariance. To obtain the required resolution for evolving similar solutions, adaptive mesh refinement (AMR) is required [26–28], which dynamically and locally adds numerical resolution. Although Van Den Berg et al. [26] present a formal analysis and AMR simulation for LLE with the Gilbert term, rigorous proof of blowup for 2-dimensional (2D) LLE remains open. However, Ding et al. [29] proved that or dimensional LLE with the Gilbert term will lead to a finite time blowup under specific initial boundary conditions. Some regularity and blowup results for LLE were derived considering the Maxwell field [30, 31].

Although the blowup problem is clear for specific settings (initial state or initial boundary conditions), the blowup problem for ILL is unclear. The one-dimensional ILL will not form a singularity as the inhomogeneous term is periodic [32]. However, if the inhomogeneous term is in some other specific format [33], blowup occurs for the inhomogeneous HMF.

We studied the blowup and energy gap for dimensional aperiodic ILL on target and investigated what happens as tends to infinity when the initial data is smooth and sufficiently large, in particular whether the solution develops distinct behaviors (finite or global time singularity) under the large data. Global smooth (or blowup) theory for ILL was not established, but we discuss some special solutions that form a singularity in finite time and classify these solutions and analyses based on their energy density.

In Section 2, we obtain a blowup solution for the ILL and derive the exact decay rate of the energy density for that solution. In Section 3, we obtain another local blowup solution that contains an energy gap under the inhomogeneous term and classify this solution into four types. In Section 4, we obtain the decay rate and prove the decreasing (or increasing) property of the local blowup solution for one of these types under some specific coefficient settings. Finally, in Section 5, we summarize the paper and present our conclusions.

2. Blowup Solution on a Sphere

2.1. Blowup Solution Derivation

Selecting the appropriate solution to construct the blowup solution is difficult. Group invariance (under or ) can be applied to search for exact solutions of LLE [34–36], but ILL will not admit any group invariance. Following [19, 20, 37], we use cylindrical coordinates,

where , , and are functions to be determined.

For simplicity, consider the case where does not contain and denote , , and as , , and , respectively. Then, from (9), (7) will be greatly simplified and maybe is transformed into

and

which are nonlinear partial differential equations. To solve these equations, let us assume that is a variable separation function:

where and are dependent on and , respectively; and is a constant. Since does not contain , and are denoted as and , respectively. Substituting (12) into (10) and (11),

and

where .

Solving (13)–(15),

and

where and are constants, is a function of , and satisfies

where is a constant.

The derivation of (20) is

If and are nonzero constants, regardless of the value of , the left sections of (20) and (21) are continuous functions, and (21) is nonzero. Therefore, according to the implicit function theorem, there exists a function that satisfies (20), and we have the following conclusion.

Theorem 1. Assuming that takes the form of (18) and satisfies (20), there is a solution of (7):

where , , , and are arbitrary constants; and constant .

On , from (9),

and of (23) is

which indicates that (22) is a solution of (7) which blows up at .

The energy density of (22) will also blow up, as can be seen from the energy density:

2.2. Decay Rate of the Blowup Solution

From the solutions and the energy density of the upper section, it is difficult to see the decay behavior of energy density. The decay characteristics of the solution determine whether solution energy can be finite value across the whole space. Therefore, we investigate the decay rate of the solution using the implicit solution from Section 2.1.

If we take the derivative of (20) with respect to ,

and, from (17), (24), and (25),

From (26), , so if and ,

where, in this case, is a strictly increasing function in the interval .

Substituting   () into the left side of (20),

If and , then, from (20),

and if ,

Equations (30) and (31) mean that determines the limit to be a positive or a negative value of . Hence,

Combining and (27),

From (33), as , the decay rate of the energy density is

where is function of ; and as ,  .

Thus, we may obtain the decay rate of the case. If and , then the decay rate of the energy density can be similarly proven. Therefore,

where satisfies

and from (27), if , we may estimate the solution

where is a function of , where if and only if , then  , which has the limit

Thus, is a continuous function with the limit at of . Hence, whether the energy is a finite value depends on the convergence of at . From (20), if ,

where is constant. Thus,

and if , then

We can obtain the initial energy density from (41) by setting :

From (40), the initial energy density is a finite value and if and only if ,   will be in infinite value. Combining this and the decay behavior as , the energy of the finite spatial area can be estimated as

where is a constant and if and only if . Furthermore, the total energy of the spatial region is

From the above analysis, we may conclude that, assuming that takes the form of (18) and satisfies (20), there is the following solution of (7).

Theorem 2. If , takes the form of (18) and satisfies (20), the energy density of (22) is

as

and the energy in satisfies

where is constant; and depends on , such that as ,  .

Similar to the proof of [20], the decay rate of the norm of gradient can now be estimated. Combining (18) and (26),

and to estimate the decay rate as , we substitute   ():

Hence,

The gradient morn and energy density are related:

so that, combining Theorem 1 and , as ,

where is function of , such that if and only if ,  .

3. Multiple Branches of the Blowup Solutions

We investigate the classification of the blowup solution. Since (7) may contain many solutions different from (22), we extend (22) to a general form to find other solutions. In particular, we use (9) as an undetermined solution and assume that   (), , and are constant. We set

where

and

Therefore,

If we set

then we can obtain a solution of (56) which satisfies the initial condition :

The derivative of the left of (58) is

We set , , and . Thus, (59) is if and only if , and

Regardless of value, (59) and the left section of (58) are continuous functions, and (59) is not zero. According to the existence theorem of implicit functions, for any

there is an that satisfies (58).

From the above process, we obtain the following conclusion.

Theorem 3. Assuming that takes the form of (18) and satisfies (58), then there is a solution of (7):

where , , , , and are any constants; and is constant.

The energy density of the solution from Theorem 3 is

Although (62) is similar to (22), their respective energies ((63) and (27)) are different. In (62), satisfies (58) which contains , and this defines the range of :

If

then

For (56), if takes the form of (66), the right of (56) is , and is .

We must solve

which can be expressed as

If ,

and

Solving

we obtain

Similarly, solving

we obtain

If and , then the inequality is

and if and , then

Thus, if ,  ,  , and

then (68) holds, and since ,

If ,  ,  ,  , and

then

Similarly, if , we can solve from , and we can solve for the (or ) case. From (58) and (56), if , the change of is shown in Table 1.

Table 1: Evolution of where and .

,  
,   

If ,
If ,

If ,

If ,
If ,

If ,

Similar to Table 1, if , the change of is shown in Table 2.

Table 2: Evolution of with and .

,
,

If ,
If ,

If ,

If ,
If ,

If ,

From Tables 1 and 2, satisfying (58) can be subdivided into the following categories:

(1),   and   (or ,   and ).  is a continuous curve. If , solving (58),

 and (56) can be expressed as

 If the right section of (82) is , then . Hence, as . Thus, the tangent of is perpendicular to the axis at . If , then from (56). Hence, the tangent of is parallel to the axis at . Specifically, in the fourth quadrant of the coordinate system, is monotonically increasing and . In the first quadrant, is monotonically decreasing and . For example, if ,  , and , then (58) is

 and the evolution of , determined by (83), is shown in Figure 1.(2),   and (or ,   and ). The evolution of in the first and the fourth quadrant are opposite that of case . In the fourth quadrant, is monotonically decreasing, and in the first quadrant, is monotonically increasing. The tangent of is perpendicular to the axis at . If we set

 the maximum value of in first quadrant will not exceed . Specifically, if ,  ,  , and , then (58) is

 and the evolution of from (85) is shown in Figure 2.(3),   and (or , and ).  is completely in the fourth quadrant, and is monotonically increasing. The maximum value of will not exceed , and the curve in this quadrant has lower bound . For example, when ,  ,  ,  , and , (58) is

 and Figure 3 shows the evolution of from (86) with and .(4),   and (or ,   and ).  is monotonically decreasing in the fourth quadrant, with upper bound .

Figure 1: Evolution of from (83) with and .

Figure 2: Evolution of from (85) with and .

Figure 3: Evolution of from (86) with and .

Cases and are multi-value, i.e., for each , There are two values of . We call these multi-branch functions. However, there is only a single monotonic function for cases and .

Let us select some specific parameters and investigate the characteristic solutions. If ,  ,  , and , then (58) is

and the evolution of (87) is shown in Figure 4.

Figure 4: Evolution of from (87) with and .

In (56), we may regard as a function of and and decompose (56) as

and

and cannot be expressed as a Hamiltonian system; that is, there is no function such that

Therefore, (90) are noncompatible systems. Solving the first equation of (90),

where is a function that is only dependent on ; and solving the second equation of (90),

where is a function that is only dependent on .

Since (91) and (92) are contradictory expressions, (90) is an incompatible system.

Equations (88)-(89) contain more solutions than (58). The evolution of these equations can also be classified into four types as follows:

1.Multibranch elliptic case I:  ,  , and (or ,  , and ).2.Multibranch hyperbolic case I: ,  , and (or ,  , and ).3.Multibranch hyperbolic case II: ,  , and (or ,  , and ).4.Multibranch elliptic case II: ,  , and (or ,  , and ).

To see the evolution of these four solutions directly, we set and , , and as follows:

1..2. and  .3.,  , and  .4. and  .

Figures 5–8 show the orientation fields for these four cases in the coordinate, respectively.

Figure 5: Multibranch elliptic case I: evolution of ().

Figure 6: Multibranch hyperbolic case I: distribution of ( and  ).

Figure 7: Multibranch hyperbolic case II: evolution of (,  , and  ).

Figure 8: Multibranch elliptic case II: distribution of ( and  ).

The direction of vector field in Figure 5 is opposite to that of Figure 8, and the direction of Figure 6 is opposite to that of Figure 7. Figures 5 and 8 (or Figures 6 and 7) are very similar. However, in Section 4, we show that the energy densities suggest that all of these four multibranch cases are distinct.

4. Blowup Solution and Energy Gap

The solutions presented in Sections 2 and 3 will blow up in finite time, but the energy of the solution presented in Section 2 is infinite on all spatial regions. However, can (7) (or (6)) contain some finite energy solutions on the sphere, since solution (62) presented in Section 3 contains some finite energy subcases? In this section, we investigate the blowup solution (62) and particularly the finite energy subcases.

Let us analyze the properties of the energy density (63) of solution (62). Since (63) contains more parameters than (27), the evolution of is more complex for (56) than (26), which leads to more complex evolutionary behavior of (63).

If ,  , and (or ,  , and ), then is negative in , as shown in Figures 4 and 8. Substituting   () into the left of (58), we obtain a function , with the limit

where .

From (93), the positive and negative of can be determined by ; that is,

Substituting into (63), the decay rate of the energy density as is

where is a function of , such that as ,  .

Similar to the proof of Theorem 2, we have the following corollary.

Corollary 4. If satisfies (58), ,  , and (or ,  , and ), then ; the energy density of (62) is

and the energy in satisfies

where is constant; and is only dependent on , such that as ,  .

In Corollary 4, parameters , , and belong to one type from Table 1 (or Table 2). The other three parameters types were divided into two types: and .

If , then (63) may be expressed as

and the monotonicity of (63) and that of (98) are the same under . Substituting into (98),

Let us separate one equation from this, such as

and solve to obtain

Thus, if , the left section of (100) achieves its maximum or minimum value.

The derivative of the left of (99) is

If , then the exact relationship of and (or and ) is shown in Figure 9 (Figure 10).

Figure 9: Relationship of and with

Figure 10: Relationship of and with .

Figure 11 shows the evolution of energy density from (63) for multi ranch elliptic case I with and , and Figure 12 shows multibranch hyperbolic case I with ,  , and . In these two cases, and . The behavior is shown for , and the function .

Figure 11: at different times with and .

Figure 12: in different time under ,  , and .

The multibranch elliptic case I (Figure 11) is a local solution that terminates at some point of . The multibranch hyperbolic case I (Figure 12) is also a local solution. However, there are some energy gaps on . Figure 13 shows the evolution of multibranch hyperbolic case II with and is also local with some energy gaps.

Figure 13: for different times with and .

The situation for the multibranch elliptic case II is different from the other three cases. As predicted in Corollary 4, this solution diverges, as shown in Figure 14 with and . is a global solution and is close to zero as . Near the peaks, there are also some energy gaps.

Figure 14: for different times with and .

We present some connections of ILL with real physics in the ending of this section. Clearly, we return the continuous ILL model back to the discrete Heisenberg model in the real physics. We plan to use the inhomogeneous Heisenberg chain model to discuss some physically relevant potential then.

ILL model is the continuum limitation of the discrete inhomogeneous Heisenberg chain, where the corresponding Heisenberg Hamiltonian (representing nearest neighbor interaction between and ) is as follows:

where and are on-site spins and denotes the on-site ferromagnetic parameter which is the inhomogeneous physical interactions at the position .

Equation (7) is a continuum model that can predicate any value of the site spin vector in a mathematical simulation way. Equation (7) plays a similar role as a Schrödinger equation for the Heisenberg model in quantum physics. As we can see in Section 1, the continuum limit of this system is also a nonlinear inhomogeneous Schrödinger equation (8) under the radial coordinates. For the discrete inhomogeneous Heisenberg chain, or has been used in the description of the magnetic properties and spin states of conjugated hydrocarbons and a large variety of metallic and organometallic compounds. What does the blowup solution (terminated in a finite time) represent in this discrete lattice system? In the continuum model ((7) or (8)), the blowup solution means that the gradient will be an infinite value, while the time is a finite time T. Accordingly, this finite time means that there are distortion and destruction in every position of sites of a lattice. In this time , there is a phase transition in the ferromagnetic solid.

Theorem 2 indicates that the phase transition (distortion and destruction) is possible under some special . More clearly, if the exact quantitative scale of is as , this breaking down will happen in the ferromagnetic material. From a discrete point of view, represents a correspondence of by spatial position. The distribution means that the value of is the value of . Hence, we expect that the phase transition in real ferromagnetic structure () will begin in some blowup time under the setting of (as ).

Multiple branches solution (see solution (62) in Section 3) allows a finite time blowup too. Accordingly, the discrete inhomogeneous Heisenberg chain of this system will generate a phase transition under some special setting of . In particular, if and (as ), the phase transition is a possible case in the multibranch elliptic case II according to Corollary 4. However, although (22) and (62) share a same breakdown in their discrete inhomogeneous Heisenberg chain system, the situation of multiple branches solution (62) is more complex than solution (22). For example, an energy gap can be found in the multibranch elliptic case II, while there is not any gap in solution (22). Why the energy gap will form in solution (62)? The most important reason is that there is a function in (58). If touches the edges of the definition domain of , jumps discontinuously from one value to another. In fact, if touches the edges of the definition domain of , will go to another track that belongs to the same multibranch. We can see that the difference between (22) and (62) is just a phase shift in the solution. From a point of the discrete inhomogeneous Heisenberg chain, the energy gap is inevitable due to the phase shift. This means that if we force the phase of the site spin vector to do a parallel movement, the jumping of the spin vector can be found in some location in any time before the blowup time. Furthermore, the scale of the jumping enlarges as the time goes on.

In the homogeneous Heisenberg chain (), all its lattice sites are equivalent. However, the self-consistent many-body ground states need not preserve all symmetries of the Hamiltonian in the inhomogeneous case. After the balance is broken, some singularity and jumping behaviors are possible according to the examples presented in the preceding sections. As we can see in these solutions, the Hamiltonian itself comprises inequivalent lattice sites that lead to the blowup and the energy gap. The inequivalent lattice sites are the case, for example, the Heisenberg chains with impurities or, more generally, the systems with different types of magnetic atoms. In the real problems, it is very common to find inhomogeneities rather than the homogeneous case. Actually, external nonuniform magnetic fields, anisotropic crystal fields, finite-size effects, boundary conditions, local impurities, lattice defects, and so forth can be induced for generalized inhomogeneities instances.

5. Conclusions

We investigated the evolution of the singular solution of dimensional ILL on . Two different ansatzes were applied in constructing the exact solutions in radially symmetric coordinates. Accordingly, two singular solutions were obtained ((22) and (62)). If one coefficient of (62) is set to , (62) degenerates into (22). All solutions show the same blowup rate . However, (22) is a global solution on the entire spatial domain, whereas (62) contains some local and global subcases.

Energy density plays an important role in the energy analysis of these two solutions. The energy density decay rate (22) scales as , which implies that the initial energy of this solution will be finite for . However, (62) is somewhat more complicated, with four subcases. If the coefficients of the solution (, , and in (62)) are suitably selected, the solution can be divided into four types: multibranch elliptic cases I and II and multibranch hyperbolic cases I and II. In contrast to (22), these subcases may lead to a new energy gap class of exact solutions.

One of these four cases is the global solution in the spatial domain. In this case, since the decay rate also scales as , the initial energy is finite. The other three subcases can be local solutions, which means that the initial energy will be finite for any spatial dimension. All these multibranch cases lead to energy gaps in the energy density. The evolution of energy density for the cases was shown in Section 3, illustrating the decay property and monotonic increasing and decreasing behaviors. These evolutions were consistent with the theorem and corollary provided in Section 3. Thus, the multibranching of for the solution plays a decisive role.

Integrability of the physical ILL has been investigated in [10]. The ILL discussed in this paper does not belong to the integrability system due to the inhomogeneous term. While many integrable spin systems are known, the integrability of IIL in the continuum limit has not been investigated. It will be of considerable interest to investigate the underlying singularity structures and energy gaps of such nonintegrable nonlinear evolution equations in detail. It would be also instructive to explore whether similar singular solutions can be obtained for LLE with Gilbert damping.

In conclusion, ILL and LLE display very interesting dynamic structures and further work is required to identify them.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author declares that he has no conflicts of interest.

Acknowledgments

This work is supported by NSFC (no. 11161057), National Science Foundation for Young Scientists of China (no. 11301014), TianYuan Special Fund of the National Natural Science Foundation of China (no. 11426068), Project for Young Creative Talents of Ordinary University of Guangdong Province (no. 2014KQNCX228), the PhD Start-up Fund of Natural Science Foundation of Guangdong Province (no. 2014A030310330), Research Award Fund for Outstanding Young Teachers in Guangdong Province (no. Yq20145084601), and Fund for Science and Technology of Guangzhou (no. 201607010352). The author is also grateful for the sponsorship from China Scholarship Council.

References

L. D. Landau and E. M. Lifshitz, “On the theory of the dispersion of magnetic permeability in ferromagnetic bodies,” in Reproduced in Collected Papers of L. D. Landau, Z. Sowjetunion, Ed., vol. 8, pp. 101–114, Pergamon, NY, USA, 1965.

X. He, A. Qian, and W. Zou, “Existence and concentration of positive solutions for quasi-linear Schrodinger equations with critical growth,” Nonlinearity, vol. 26, no. 12, pp. 3137–3168, 2013.

Y. Sun, L. Liu, and Y. Wu, “The existence and uniqueness of positive monotone solutions for a class of nonlinear Schrödinger equations on infinite domains,” Journal of Computational and Applied Mathematics, vol. 321, pp. 478–486, 2017.

X. Zhang, L. Liu, Y. Wu, and Y. Cui, “Entire blow-up solutions for a quasilinear -Laplacian Schrödinger equation with a non-square diffusion term,” Applied Mathematics Letters, vol. 74, pp. 85–93, 2017.

X.-N. Ma, P.-H. Wang, and W. Wei, “Constant mean curvature surfaces and mean curvature flow with non-zero Neumann boundary conditions on strictly convex domains,” Journal of Functional Analysis, vol. 274, no. 1, pp. 252–277, 2018.

F. Sun, L. Liu, and Y. Wu, “Infinitely many sign-changing solutions for a class of biharmonic equation with p-Laplacian and Neumann boundary condition,” Applied Mathematics Letters, vol. 73, pp. 128–135, 2017.

P. Wang and X. Wang, “The geometric properties of harmonic function on 2-dimensional Riemannian manifolds,” Nonlinear Analysis: Theory, Methods & Applications, vol. 103, pp. 2–8, 2014.

P. Wang and L. Zhao, “Some geometrical properties of convex level sets of minimal graph on 2-dimensional Riemannian manifolds,” Nonlinear Analysis: Theory, Methods & Applications, vol. 130, pp. 1–17, 2016.

R. Balakrishnan, “On the inhomogeneous Heisenberg chain,” Journal of Physics C: Solid State Physics, vol. 15, no. 36, pp. L1305–L1308, 1982.

M. Daniel, K. Porsezian, and M. Lakshmanan, “On the integrability of the inhomogeneous spherically symmetric Heisenberg ferromagnet in arbitrary dimensions,” Journal of Mathematical Physics, vol. 35, no. 12, pp. 6498–6510, 1994.

M. Han, L. Sheng, and X. Zhang, “Bifurcation theory for finitely smooth planar autonomous differential systems,” Journal of Differential Equations, vol. 264, no. 5, pp. 3596–3618, 2018.

F. Li and Q. Gao, “Blow-up of solution for a nonlinear Petrovsky type equation with memory,” Applied Mathematics and Computation, vol. 274, pp. 383–392, 2016.

F. Meng and J. Shao, “Some new Volterra-Fredholm type dynamic integral inequalities on time scales,” Applied Mathematics and Computation, vol. 223, pp. 444–451, 2013.

H. Tian and M. Han, “Bifurcation of periodic orbits by perturbing high-dimensional piecewise smooth integrable systems,” Journal of Differential Equations, vol. 263, no. 11, pp. 7448–7474, 2017.

L. Li, F. Meng, and Z. Zheng, “Oscillation results related to integral averaging technique for linear Hamiltonian systems,” Dynamic Systems and Applications, vol. 18, no. 3-4, pp. 725–736, 2009.

L. Li, F. Meng, and Z. Zheng, “Some new oscillation results for linear Hamiltonian systems,” Applied Mathematics and Computation, vol. 208, no. 1, pp. 219–224, 2009.

J. Liu and Z. Zhao, “Variational approach to second-order damped Hamiltonian systems with impulsive effects,” Journal of Nonlinear Sciences and Applications. JNSA, vol. 9, no. 6, pp. 3459–3472, 2016.

Z. Zheng, “Invariance of deficiency indices under perturbation for discrete Hamiltonian systems,” Journal of Difference Equations and Applications, vol. 19, no. 8, pp. 1243–1250, 2013.

P. Zhong, S. Wang, and S. Chen, “Some periodic and blow-up solutions for Landau-LIFshitz equation,” Modern Physics Letters A, vol. 26, no. 32, pp. 2437–2452, 2011.

P. Zhong, S. Wang, and M. Zeng, “Two blowup solutions for the inhomogeneous isotropic Landau-LIFshitz equation,” Journal of Mathematical Analysis and Applications, vol. 409, no. 1, pp. 74–83, 2014.

X.-G. Liu, “Concentration sets of the Landau-LIFshitz system and quasi-mean curvature flows,” Calculus of Variations and Partial Differential Equations, vol. 27, no. 4, pp. 493–525, 2006.

C. Wang, “On Landau-LIFshitz equation in dimensions at most four,” Indiana University Mathematics Journal, vol. 55, no. 5, pp. 1615–1644, 2006.

H. Huh, “Blow-up solutions of modified Schrödinger maps,” Communications in Partial Differential Equations, vol. 33, no. 1-3, pp. 235–243, 2008.

Q. Ding, “Explicit blow-up solutions to the Schrödinger R2 to the hyperbolic 2-space H2,” Journal of Mathematical Physics, vol. 50, no. 10, 103507, 17 pages, 2009.

F. Merle, P. Raphaël, and I. Rodnianski, “Blowup dynamics for smooth data equivariant solutions to the critical Schrödinger map problem,” Inventiones Mathematicae, vol. 193, no. 2, pp. 249–365, 2013.

J. B. van den Berg and J. F. Williams, “(In-)stability of singular equivariant solutions to the Landau-LIFshitz-Gilbert equation,” European Journal of Applied Mathematics, vol. 24, no. 6, pp. 921–948, 2013.

W. Ren and X.-P. Wang, “An iterative grid redistribution method for singular problems in multiple dimensions,” Journal of Computational Physics, vol. 159, no. 2, pp. 246–273, 2000.

X.-P. Wang and Carlos J., “A Gauss-Seidel projection method for micromagnetics simulations,” Journal of Computational Physics, vol. 171, no. 1, pp. 357–372, 2001.

S. Ding and C. Wang, “Finite time singularity of the Landau-LIFshitz-Gilbert equation,” International Mathematics Research Notices, no. 4, Art. ID rnm012, 25 pages, 2007.

S. Ding and J. Lin, “Partially regular solution to Landau-LIFshitz-Maxwell equations in two dimensions,” Journal of Mathematical Analysis and Applications, vol. 351, no. 1, pp. 291–310, 2009.

J. Lin and S. Ding, “Blow-up solutions to Landau-LIFshitz-Maxwell systems,” Mathematical Methods in the Applied Sciences, vol. 34, no. 8, pp. 1006–1024, 2011.

J. Lin and S. Ding, “Smooth solution to the one-dimensional inhomogeneous non-automorphic Landau-LIFshitz equation,” Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 462, no. 2072, pp. 2397–2413, 2006.

Y. Li and Y. Wang, “Bubbling location for F-harmonic maps and inhomogeneous Landau-Lifshitz equations,” Commentarii Mathematici Helvetici, vol. 81, no. 2, pp. 433–448, 2006.

G. Yang, “The difference between Schrodinger equation derived from Schrödinger map and Landau-Lifshitz equation,” Physics Letters A, vol. 376, no. 4, pp. 231–235, 2012.

G. YANG, “Spherical cone symmetric families generated by Landau-Lifshitz equation and their evolution,” SCIENCE CHINA Mathematics, vol. 41, no. 2, p. 181, 2011.

G. Yang and B. Guo, “Some exact solutions to multidimensional Landau-LIFshitz equation with uprush external field and anisotropy field,” Nonlinear Analysis, vol. 71, no. 9, pp. 3999–4006, 2009.

B. Guo, Y. Han, and G. Yang, “Blow up problem for Landau-LIFshitz equations in two dimensions,” Communications in Nonlinear Science & Numerical Simulation, vol. 5, no. 1, pp. 43-44, 2000.