1 Introduction

In order to describe the continuum evolution of the film-free surface, the authors in [5] proposed the following classical surface diffusion equation:

where and T are the normal surface velocity, the surface diffusivity, the number of atoms per unit area on the surface, the atomic volume, the molar volume of lattice cites in the film, the universal gas constant and the absolute temperature, respectively.. Moreover, ∆_S is the surface Laplace operator, ν represents the regularization coefficient that measures the energy of edges and corners, C_αβ means the surface curvature tensor, and µ_w being an exponentially decaying function of u that has a singularity at u → 0 (see [5]).

Particularly, in the small-slop approximation, the cases of a crystal, which has cubic symmetry and high-symmetry orientations, arise an evolution equation in the following form for the film thickness:

where are three smooth functions, respectively. From the physical consideration, Eq. (1) is supplemented with the following boundary value conditions:

and initial condition

We remark that Eq. (1) is related to the formation of quantum dots in epitaxially grown thin solid films. This formation has been attracting attention as a very promising area of nanotechnology that can lead to a new generation of electronic devices. According to the mechanism for the formation, the substrate induces the film growth in a certain crystallographic orientation. In the absence of wetting interactions with the substrate, due to a large surface-energy anisotropy, this orientation would be thermodynamically forbidden, and the surface would undergo a long-wave faceting (spinodal decomposition) instability. In [5], the authors show that wetting interactions between the film and the substrate can suppress this instability and qualitatively change its spectrum, leading to the damping of long-wave perturbations and the selection of the preferred wavelength at the instability threshold. This creates a possibility for the formation of stable regular arrays of quantum dots even in the absence of epitaxial stresses.

In [15], on the basis of the Schauder-type estimates and the techniques in Campanato spaces, the author assumed that the smooth functions satisfy

and studied the existence of classical global solutions.

The dynamic properties of diffusion systems such as the global asymptotical behaviors of solutions and global attractors are important for the study of diffusion models, which ensure the stability of diffusion phenomena and provide the mathematical foundation for the study of diffusion dynamics. For the higher-order diffusion equation, there are many classical results related to its global attractors, see, for example, Dlotko [3], Li and Zhong [7], Schimperna [10], Alouini [1], Cheskidov and Dai [2], Huang, Yang, Lu and Wang [6], Duan and Xu [4] and so on. To the best of our knowledge, the existence of the global attractor for problem (1)–(3) has not been addressed yet, which is the main goal of this article.

In this paper, by using the regularity estimates for the semigroups, iteration technique and the classical existence theorem of global attractors we consider the global attractor of solutions for the initial boundary value problem for Eq. (1). The main results are Theorem 1 in the next section and Theorem 2 in Section 3.

This article is organized as follows. In Section 2, we give some preparations for our consideration and prove the existence of global attractors for problem (1)–(3) in the Sobolev space H³(Ω). In Section 3, we prove the existence of global attractors for problem (1)–(3) in the Sobolev space H^k(Ω) with any k ≥ 0

2 Global attractor in H³

In order to consider the global attractors of problem (1)–(3), we select a closed metric space and prove the existence of global attractors of problem (1)–(3) in this space, where.

For convenience, we give the following lemma on global existence and uniqueness of solution to problem (1)–(3).

Lemma 1. Assume that and the functions satisfy (4). Then (1)–(3) admits a unique global solution, which satisfies

The proof of existence and regularity of solutions is based on the Galerkin method and a priori estimates in the following. Thanks to the above existence lemma, we know that there exists a continuous operator semigroup.

Then, by the classical existence theorem of global attractors (see [14]), we give the following theorem on the existence of the global attractor of problem (1)–(3) in H³(Ω).

Theorem 1. Assume and the three smooth functions satisfy (4). Then the solution u of problem (1)–(3) possesses a global attractor in the space H³(Ω), which attracts all bounded set in the space H³(Ω).

Similar to [14], we assume that the semigroup is generated by the solutions of Eq. (1) with initial conditions . Then we give two lemmas to prove Theorem 1.

Lemma 2. There exists a bounded set whose size depends only on Ω such that for all , there exists

Proof. It suffices to prove that there is a positive constant C such that for large t, there holds. Now we begin to prove the lemma.

Let

Integrating by parts and by (2) we have

Hence, F (t) ≤ F (0), then

Then by Poincaré’s inequality and (2) we derive that

Combining (5) and (6) together gives

Notice that, then there exists such that

It then follows from (6) and (7) that

Note that . By Poincaré’s inequality and (8) we deduce that for,

Hence, by (7), (8) and (9) we get sup and . Multiplying (1) with D6u and integrating it over Ω, we obtain

By Nirenberg’s inequality we have

Then

Summing up, we have

Assuming that ε′ is small enough and it satisfies 1 − 8ε′ > 0, we have

Applying Gronwall’s inequality, we have

Therefore, for initial data in any bounded set , there is a uniform time t₂(B) depending on B such that for ,

Sobolev embedding theorem gives.Adding (7), (8), (9) and (10), we have , then the lemma is proved.

From the above lemma we know that {S(t)}_t≥0 has a bounded absorbing set in H³(Ω). In what follows, we prove the precompactness of the orbit in H³(Ω).

Lemma 3. For any initial data u0 in any bounded set, there is a T (B) > 0 such that

Proof. The uniform boundedness of H³(Ω) norm of u(t) has been achieved in Lemma 2. In what follows, we give the estimate on H⁴-norm.

Differentiating (1) with respect to x, multiplying the resultant by D⁷u and integrating on Ω, using the boundary conditions, we have

By Nirenberg’s inequality we have

Hence, we have

Summing up, we derive that

Assuming that ε′ is small enough and it satisfies 1 − 6ε′ > 0, we have

By (10) we have

Integrating (14) between t and t + 1 and by (11), we have

By Poincaré’s inequality we have

Therefore, by (13), (15) and the uniform Gronwall inequality we have

The lemma is proved.

Now we give the proof of Theorem 1.

Proof of Theorem 1. Combining Lemma 2 with Lemma 3, by [14] we have completed the proof of Theorem 1.

3 Global attractor in Hk

In order to consider the global attractors for Eq. (1) in the Hk space, we introduce the definition as follows:

In this article, we let

be a nonlinear function, and we assume that the linear operator L = D⁶ : H₁ → H in (16)–(18) is a sectorial operator, which generates an analytic semigroup e^tL, and L induces the fractional-power operators and fractional-order spaces as follows:

where Hα = D(L^α) is the domain of L^α. By the semigroup theory of linear operators, H_β ⊂ H_α is a compact inclusion for any β > α.

Now we give the main theorem of this article, which provides the existence of global attractors of Eq. (1) in any kth space H^k.

Theorem 2. Assume and the smooth functions satisfy (4), then the solution u of problem (1)–(3) possesses a global attractor in the space H^k, which attracts all the bounded set of H^k in the H^k-norm, where k ∈ [0, ∞).

On the basis of Ma and Wang [9], it is well known that the solution of problem (1)–(3) can be expressed as

where L = D⁶ and

Then (20) means

In order to prove Theorem 2, we first prove the following lemma.

Lemma 4. For any bounded set U ∈ Hα, there exists a constant C > 0 such that

Proof. For α = 1/2, this follows from Lemma 2, i.e., for any bounded set U ⊂ H_1/2, there exists a constant C, C > 0, such that

Then we only need to prove (21) for any α ≥ 1/2.

Step 1. We prove that for any bounded set U ⊂ Hα (1/2 ≤ α < 2/3), there exists a constant C > 0 such that

We claim that g : H1/2 → H is bounded, by Sobolev embedding theorem we have

Then we obtain

which means that g : H_1/2 → H is bounded. By (19), (20) and (23) we find that

where β = 1/3 + α (0 < β < 1). Hence, (22) is valid.

Step 2. We prove that for any bounded set U ⊂ H_α (2/3 ≤ α < 5/6), there exists a constant C > 0 such that

We claim that g : Hα → H1/6 is bounded, by Sobolev embedding theorem we have

Then we obtain

which means that g : Hα → H1/6 is bounded. On the basis of Step 1 and (25), we deduce that

where β = 1/6 + α (0 < β < 1). Hence, (24) is valid.

Step 3. We prove that for any bounded set U ⊂ Hα (5/6 ≤ α < 1), there exists a constant C > 0 such that

We claim that g : Hα → H1/3 is bounded, by Sobolev embedding theorem we have

Then we obtain

which means that g : H_α → H_1/3 is bounded. On the basis of Step 2 and (27), we deduce that

Hence, (26) is valid.

In the same fashion as in the proof of (26), by iteration we can prove that for any bounded set U ⊂ H_α (α > 0), there exists a constant C > 0 such that

That is, for all α ≥ 0, the semigroup S(t) generated by problem (1)–(3) is uniformly compact in H_α. Lemma 4 is proved.

Lemma 5. For any α ≥ 0, problem (1)–(3) has a bounded absorbing set in H_α. That is, for any bounded set U ∈ H_α, there exists T > 0 and a constant C > 0 independent of u₀ such that

Proof. For α = 1/2, this follows from Lemma 3. Then we prove (28) for any α > 1/2. We proceed in the following steps.

Step 1. We prove that for any 1/2 ≤ α < 2/3, (1)–(3) has a bounded absorbing set in H_α.

By (19) we have

Assume that B is the bounded absorbing set of problem (1)–(3) and B satisfies B ⊂ H_1/2. In addition, we also assume the time t₀ > 0 such that

Note that is the first eigenvalue of the equation

Then for any given T > 0 and u₀ ∈ U ⊂ H_α (α ≥ 1/2), we can obtain

Adding (23) and (29) together, we have

where C > 0 is a constant independent of u₀. Then by (30) and (31) we have that (28) holds for all 1/2 ≤ α < 2/3.

Step 2. We prove that for any 2/3 ≤ α < 5/6, problem (1)–(3) has a bounded absorbing set in H_α.

Adding (25) and (29) together, we have

where C > 0 is a constant independent of u₀. Then by (30) and (32) we have that (28) holds for all 2/3 ≤ α < 5/6.

Step 3. We can use the same method as the above step to prove that for any 5/6 ≤ α < 1, problem (1)–(3) has a bounded absorbing set in Hα. By the iteration method we can obtain that (28) holds for all α ≥ 1/2.

Now we give the proof of Theorem 2.

Proof. Combining Lemmas 4 and 5, we have completed the proof of Theorem 2.

4 Conclusion

In the study of a mechanism for the formation of quantum dots on the surface of thin solid films, there arise a sixth-order parabolic equation (1) (see [5]). Mathematically, Eq. (1) is a nonlinear evolution equation. Studying the properties of solutions for Eq. (1) is so interesting and maybe useful for the study of the surface of thin solid films. Recently, Zhao [15] studied the existence of classical solutions of the initial boundary value problem of such equation. In order to study the long-time behavior of solutions, we prove the existence of global attractor of Eq. (1) with boundary and initial value conditions. The main idea comes from Temam [14] and Ma and Wang [9]. By using the properties of sectorial operator we define a semigroup related to the fractional-order Sobolev space Hk(0, 1) (k ∈ [0, ∞)). Then applying Sobolev embedding theorem and iteration technique, we obtain some useful a priori estimates and prove the existence of absorbing sets and asymptotic compactness of semigroup, obtain our main result on the existence of global attractor.

It is worth pointing out that there are also some other papers that studied the existence of global attractor for dissipative equations in fractional-order Sobolev spaces (see, e.g., semilinear parabolic equation [11], fourth-row Cahn–Hilliard equations [12,16], modified Swift–Hohenberg equations [13], sixth-row Cahn–Hilliard equations [8] and the reference cited therein). We also remark that although both [8] and this paper are focus on the sixth-order diffusion equations, the equation considered in [8] is a sixth-order convective Cahn–Hilliard equation, which is different from equation (1). Due to the different terms of both papers, the calculations are different from the paper of Liu and Liu [8].

Acknowledgments

The authors would like to express their deep thanks to the referee’s valuable suggestions for the revision and improvement of the manuscript.

Referencias

1. B. Alouini, Finite dimensional global attractor for a dissipative anisotropic fourth order Schrödinger equation, J. Differ. Equations, 266:6037–6067, 2019, https://doi.org/ 10.1016/j.jde.2018.10.044.

2. A. Cheskidov, M. Dai, The existence of a global attractor for the forced critical surface quasi- geostrophic equation in .., J. Math. Fluid Mech., 20:213–225, 2018, https://doi.org/ 10.1007/s00021-017-0324-7.

3. T. Dlotko, Global attractor for the Cahn–Hilliard equation in .. and .., J. Differ. Equations,113:381–393, 1994, https://doi.org/10.1006/jdeq.1994.1129.

4. N. Duan, X. Xu, Global dynamics of a fourth-order parabolic equation describing crystal surface growth, Nonlinear Anal. Model. Control, 24(2):159–175, 2019, https://doi. org/10.15388/NA.2019.2.1.

5. A.A. Golovin, M.S. Levine, T.V. Savina, S.H. Davis, Faceting instability in the presernce of wetting interactions: A mechanism for the formation of quantum dots, Phys. Rev. B, 70:235342, 2004, https://doi.org/10.1103/PhysRevB.70.235342.

6. L. Huang, X. Yang, Y. Lu, T. Wang, Global attractors for a nonlinear one-dimensional compressible viscous micropolar fluid model, Z. Angew. Math. Phys., 70:40, 2019, https: //doi.org/10.1007/s00033-019-1083-5.

7. D. Li, C. Zhong, Global attractor for the Cahn–Hilliard system with fast growing nonlinearity, J. Differ. Equations, 149(2):191–210, 1998, https://doi.org/10.1006/jdeq. 1998.3429.

8. C. Liu, A. Liu, The existence of global attractor for a sixth order parabolic equation, Sci. Bull., Ser. A, Appl. Math. Phys., Politeh. Univ. Buchar.,76(1):115–128, 2014.

9. T. Ma, S. Wang, Bifurcation Theory and Applications, World Sci. Ser. Nonlinear Sci., Ser. A, Vol. 53, World Scientific, Hackensack, NJ, 2005, https://doi.org/10.1142/5798.

10. G. Schimperna, Global attractors for Cahn–Hilliard equations with nonconstant mobility, Nonlinearity, 20(10):2365–2387, 2007, https://doi.org/10.1088/0951-7715/ 20/10/006.

11. L. Song, Y. He, Y. Zhang, The existence of global attractors for semilinear parabolic equation in H. space, Nonlinear Anal., Theory Methods Appl., 68(11):3541–3549, 2008, https://doi.org/10.1016/j.na.2007.03.045.

12. L. Song, Y. Zhang, T. Ma, Global attractor of the Cahn–Hilliard equation in H. spaces, J. Math. Anal. Appl., 355:53–62, 2009, https://doi.org/10.1016/j.jmaa.2009. 01.035.

13. L. Song, Y. Zhang, T. Ma, Global attractor of a modified Swift–Hohenberg equation in H. spaces, Nonlinear Anal., Theory Methods Appl., 72(1):183–191, 2010, https://doi. org/10.1016/j.na.2009.06.103.

14. R. Temam, Infinite Dimensional Dynamical Systems in Mechanics and Physics, Appl. Math. Sci., Vol 68, Springer, New York, 1988, https://doi.org/10.1088/0951-7715/ 18/5/013.

15. X. Zhao, A sixth-order parabolic equation describing continuum evolution of film free surface,Electron. J. Differ. Equ., 2014:223, 2014.

16. X. Zhao, B. Liu, The existence of global attractor for convective Cahn–Hilliard equation, J. Korean Math. Soc., 49(2):357–378, 2012, https://doi.org/10.4134/JKMS. 2012.49.2.357.

Notes