Advances in Materials Science and Engineering

Volume 2012 (2012), Article ID 520967, 5 pages

doi:10.1155/2012/520967

Research Article

Integral Solution of the Interface Profile of Grain Boundary Grooving by Surface Diffusion

Caifang Wang,¹ Hui Xing,² Fengping Yao,³ Changle Chen,² Jianyuan Wang,² and Fei Fan²

¹Department of Mathematics, Shanghai Maritime University, Shanghai 201306, China

²Shaanxi Key Laboratory for Condensed Matter Structure and Properties, Northwestern Polytechnical University, Xi'an 710072, China

³Department of Mathematics, Shanghai University, Shanghai 200444, China

Received 23 June 2012; Revised 14 October 2012; Accepted 15 October 2012

Academic Editor: Pavel Lejcek

Copyright © 2012 Caifang Wang et al. 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.

In polycrystalline thin films, formation of grain boundary grooving is investigated widely due to its effects on their properties. Infinite series solution of the interface of grain-boundary grooving by surface diffusion has been obtained. In this paper, with Fourier transform, a novel integral solution of the profile of grain boundary grooving is derived. Numerical examples are taken to test the performance of the integral solution. According to the numerical examples, we find that the novel solution is in accordance with the infinite series solution.

1. Introduction

In polycrystalline thin films, grain boundary grooving through the thickness of the film is a common failure mode that strongly affects their properties. Formation and development of grain boundary grooving have received a great deal of attention over decades [1–6]. The grooving forms and develops at the point of intersection when the grain boundary ends at a free boundary in order to reduce the total free energy. In absence of the external potential field, the grooving profile migrates via surface diffusion. Figure 1 shows the AFM image of grain boundary grooving of thin film. Mass transport by surface diffusion is driven by the curvature Laplacian. For convex surface, mass flows from high curvature parts to low parts, while in reverse for concave surfaces. The theory of the capillary-driven shape evolution was first provided by Mullins [1], and this kind of problem hence is so-called Mullins problem. In his study, he considered that the grooving is symmetric between the two grains in the small slope approximation. Assume that the surface energy is isotropic and the surface flux can be expressed as

where is the surface curvature gradient, is the coefficient of surface diffusion, is the arc length of the interface profile, is the isotropic surface energy, is the atomic volume, is the number of mobile atoms in the surface layer per unit, and is with its usual meaning. Therefore the surface normal speed can be written as

Defining as the interface profile and using , we could get

where . In small slope approximation, it is reduced to

Figure 1: AFM image of grain boundary grooving of BiFeO₃/LaNiO₃ thin film.

Wong et al. [7] have shown that the film profiles at different times are found to be self-similar. The wedge contact angle will change at initial time from wedge angle to the equilibrium . Therefore, the partial differential equation can be converted into an ordinary differential equation before solving it by the integral method. In this paper, we mainly seek the integral solution of (4) with proper initial and boundary conditions. The rest of the paper is organized as follows. In Section 2, we derive the main integral solution of the interface profile with Fourier transform. Then we change the integral solution into normalized solution with some variable substitution. In Section 3, we take numerical experiments to test the novel integral solution. In Section 4, we give the conclusion of this paper.

2. Integral Solution

In the following, we consider the symmetric case of the profile . Since the symmetry is about axis, we only consider the profile for . As in [1], we assume that the initial surface is flat. Also, we assume that at the groove root, the flux is zero and there is a fixed dihedral angle. Hence, the groove profile satisfies the following initial boundary problem:

where is the slope related to the dihedral angle enclosed by the surface.

2.1. Main Solution

According to the theories of partial differential equations and Fourier transform, we can get the integral solution of the initial boundary problem (5) as follows:

In the following, we use three steps to derive the integral solution (6).

Let . Then we only need to seek the solution of the following problems:

Step 1. We first consider the solution of

By taking Fourier transform with respect to , we can obtain the solution of (11)

Since a Schwartz function's Fourier transform is still a Schwartz function, we deduce that

is a Schwarz function for the variable , which makes the double integral sense in (12) when is a polynomial function. In fact, we can refer to the paper [8] for further discussions of elementary theory of (11).

Step 2. It is easy to verify that

satisfies (7)–(10). In fact, it is trivial to prove that satisfies (7), (9), and (10). Moreover, we can verify that also satisfies (8) by choosing

Step 3. Finally, we can obtain the solution of (5)

2.2. Normalized Solution

Let , and . With the variable substitution, the function can also be expressed as

Let

It is the Fourier transform of a Schwartz function . According to the properties of Fourier transform, is still Schwartz function and is integrable. With simple calculation, we arrive at

Finally, we get the normalized solution of the initial boundary problem (5)

3. Numerical Test

In this section, we take numerical experiments to show the performance of our integral solution of the interface profile of grain boundary grooving. All the numerical tests are implemented on a Dell desktop 4311s (Intel Core2 Duo CPU E7400 2.8 GHz with 2.00 GB RAM). The software Matlab is used to calculate all the numerical solutions.

3.1. Numerical Calculation of

We first calculate Schwartz function . To facilitate the calculation, we expand the function in (18) in a power series and integrate term by term. Then we get

where is the gamma function. Since is a series form, we use finite terms

to approximate in the numerical tests. The truncated remainder for is

We exam the dependence of truncated error and number of terms for different in with decreasing truncated error. Table 1 illustrates the details and Table 2 lists the values of and when truncated error . It is shown that the number of terms used exploded for large as truncated error decreased. However, we do not have to calculate the value of for those far away from zero since is a Schwarz function.

Table 1: Number of terms to truncated error for different .

0.1
0.5
1
2
5
10
20

2
3
5
7
13
25
54

2
4
5
8
14
26
55

3
5
6
9
16
28
57

3
5
7
10
18
31
60

Table 2: Values of and ().

0.1
0.5
1
2
5
10
20

2
3
5
7
13
25
54

3.2. Integral Solution

In the following, we provide the numerical result of our integral solution. used in the integral solution is approximated by . The part of infinite integration is truncated to finite integral over the interval from zero to 20. To calculate numerical integration, trapezium integral formula is applied. Since the parameter only impacts the scale of the solution , here we select . We plot the normalized profile in Figure 2. To test the performance of our integral solution, we also compare the integral solution with Mullins' series solution [1] at . Here Mullins solution of the interface profile is displayed as

where is calculated by the recursion relation

Figure 2 displays the comparison results. Solid line is Mullins' solution and dash line is our integral solution. Monitoring the evolution of the depth and the width of the groove profile, we find that our solution agrees very well with Mullins' solution according to the figure.

Figure 2: Normalized profile shape due to volume diffusion.

Here we have not compared the groove profile at large because both solutions can easily become unstable due to the limited machine precision and the truncation. For Mullins' solution, we plot the profile between with different truncation of the infinite series in Figure 3. The truncation of the infinite series is chosen as , respectively. From Figure 3, we find that for , and 42, the curves of Mullins's solution grow to positive infinity and for , and 44, the curves grow to negative infinity. This phenomenon appears in a four-term cycle, and it is in agreement with the coefficients of the infinite series in (25) and (26). However, for our integral solution, with the properties of the Fourier transform, we have

Hence the groove profile at a fixed time when . Since is a Schwarz function, we can easily conclude that is stable with the properties above, and we do not have to examine the profile at large .

Figure 3: Normalized profile shape due to volume diffusion.

3.3. Properties of Integral Solution

At last, we conduct some properties of our integral solution (20). We compute

and find that the maximum of the surface is at the first occurrence of with

Combined with the data in Figure 2, the maximum of the surface is bout . The depth of the groove measured in direction from the maximum of the surface to the grain boundary is given by

whereas the separation of the two maxima is

From these expressions, we have the time independent ratio

The spatial oscillation of the groove profile at a fixed time is shown in Figure 4. The solution shows a spatial oscillation with a rapid exponential damping in the amplitude. Beyond this, we can also get the interval between the groove root and the beginning of the tail from the figure.

Figure 4: The spatial oscillation of the groove profile.

The curvature of the groove profile is defined as . However, in small slope approximation, we assume that . Hence, the curvature of the normalized profile here is . The result is plotted in Figure 5. All these properties are in agreement with Mullins' solution.

Figure 5: The curvature of the groove profile.

4. Conclusion

In this work, we mainly seek a novel solution of the interface profile of grain-boundary grooving by surface diffusion with small slope approximation. With Fourier transform theory and properties of partial differential equation, we derive an integral solution (6) for the Mullins problem. Taking variable substitution, we get the normalized solution (20). According to the numerical tests, we find that our solution is identical with Mullins' solution. The rates of growth of the depth and width remain proportional to . Besides, we can get the properties of the groove profile with theoretical analysis, and these properties are no longer affected by the machine precision and the truncation of infinite series.

Acknowledgments

C. Wang is supported by Science and Technology Program of Shanghai Maritime University no. 20110054. H. Xing, C. Chen and F. Fan are supported by the National Natural Science Foundation of China Grant no. 61078057 and 50702046.

References

W. W. Mullins, “Theory of thermal grooving,” Journal of Applied Physics, vol. 28, no. 3, p. 333, 1957.

W. W. Mullins, “Grain boundary grooving by volume diffusion,” Transactions of the Metallurgical Society of the AIME, vol. 218, pp. 354–361, 1960.

H. Zhang and H. Wong, “Coupled grooving and migration of inclined grain boundaries: regime I,” Acta Materialia, vol. 50, pp. 1983–1994, 2002.

M. Bouville, C. Dongzhi, and D. J. Srolovitz, “Grain-boundary grooving and agglomeration of alloy thin films with a slow-diffusing species,” Physical Review Letters, vol. 98, no. 8, Article ID 085503, 2007.

M. Bouville, “Effect of grain shape on the agglomeration of polycrystalline thin films,” Applied Physics Letters, vol. 90, no. 6, Article ID 061904, 3 pages, 2007.

L. Klinger, “Surface evolution in two-component system,” Acta Materialia, vol. 50, no. 13, pp. 3385–3395, 2002.

H. Wong, M. J. Miksis, P. W. Voorhees, and S. H. Davis, “Capillarity driven motion of solid film wedges,” Acta Materialia, vol. 45, no. 6, pp. 2477–2484, 1997.

F. Yao, “Schauder estimates for parabolic equation of bi-harmonic type,” Applied Mathematics and Mechanics, vol. 28, no. 1, pp. 1503–1516, 2007.