Grain boundary diffusion in a PeierlsNabarro potential
Abstract
We investigate the diffusion of a grain boundary in a crystalline material. We consider in particular the case of a regularly spaced lowangle grain boundary schematized as an array of dislocations that interact with each other through longrange stress fields and with the crystalline PeierlsNabarro potential. The methodology employed to analyze the dynamics of the center of mass of the grain boundary and its spatiotemporal fluctuations is based on overdamped Langevin equations. The generality and the efficiency of this technique is proved by the agreement with molecular dynamics simulations.
1 Introduction
Understanding interface kinetics in materials is an important theoretical and practical problem, since this process influences the microstructure, such as the grain size, the texture, and the interface type. From the theoretical point of view, the study of processes that involve surface and interface properties gained significant interest in nonequilibrium statistical mechanics [1, 2]. In particular, dislocations [3, 4] and grain boundaries [5, 6] provide a concrete example of driven elastic manifolds in random media [7]. Other example of this general problem are domain walls in ferromagnets [8, 9], flux line in type II superconductors [10, 11], contact lines [12, 13] and crack fronts [14, 15]. From the point of view of applications, understanding grain boundary kinetics has a great importance for polycrystalline materials, since the resulting grain microstructure determines material properties such as strength, hardness, resistance to corrosion, conductivity etc. [16]. Hence the ambitious goal of these studies is to be able to control the microstructural properties of polycrystals.
Several approaches have been employed in the literature to study grain boundary kinetics. Ref. [17] employs molecular dynamics (MD) simulations with appropriate interatomic interactions to study the diffusion of grain boundaries at the atomic scale [17]. The method allows to quantify the mobility of grain boundaries and to compare the results with experiments [17]. While MD simulations provide a very accurate description of the dynamics, the method suffers from numerical limitations and it is difficult to reach the asymptotic regime. An alternative method is provided by the Langevin approach in which the grain boundary is assumed to evolve stochastically in an external potential [18]. The dynamics of the underlying crystalline medium enters in the problem only through the noise term (due to lattice vibrations) and the periodic potential (PeierlsNabarro). Hence, the equations of motion of the atoms or molecules are not directly relevant. Indeed there is experimental evidence in supporting of separation of time scales in plastic flow [19] and it is thus possible to integrate out the fast degrees of freedom (atomic vibrations) and consider only the slow ones (dislocations position).
Here we study the evolution of a grain boundary (GB) in a crystalline material by the Langevin approach. The GB is treated as an array of interacting dislocations performing a thermally activated motion in a periodic (PeierlsNabarro) potential. Similar models have been employed in the past to study the conductivity of superionic conductors [20, 21], the relaxational dynamics of rotators [22] and Josephson tunneling junctions [18]. Notice that the crucial role played by longrange stresses is often disregarded in analyzing GB deformation. On the other hand, it has been shown in Ref. [5] that a surface tension approximation for the GB stiffness is inappropriate and one has to consider explicitly nonlocal interactions. The present model incorporates this effect in the equations of motion.
We simulate the set of Langevin equations numerically to describe the GB kinetics and its fluctuations. The results are in good agreement with MD simulations [17] and allow to clarify the origin of the short time deviations from the diffusive behavior observed in Ref. [17]. In addition, a linearized version of the model can be treated analytically and the asymptotic results are found in good agreement with the simulations. The manuscript is organized as follows: in Sec. II we introduce the model, which is first studied in the flat GB limit in Sec. III. Sec. IV presents numerical simulations of the full flexible GB problems and Sec. V discusses the continuum theory. Sec. VI is devoted to conclusions.
2 The model
To study the GB dynamics we consider a phenomenological mesoscopic
approach. We consider in particular the case of a regularly spaced lowangle
grain boundary schematized as an array of straight dislocations that interact with each other
through longrange stress fields and with the crystalline PeierlsNabarro (PN) potential.
The GB is composed by dislocations where configurations are repeated ad infinitum because of periodic boundary conditions along the direction.
Each dislocation has Burger vector of modulus parallel to the axis and
the distance between two adjacent dislocations along the direction is
fixed to be .
Each straight dislocation interacts with the lattice and with others
dislocations through longrange stress fields.
The effect of the lattice over each nth dislocation can be decomposed as the
sum of three contributions:

, the PN force where is the area of the GB, is the shear modulus and the interatomic distance;

, the average effect of the lattice fluctuations where is the viscosity coefficient;

, the impulsive effect of the lattice fluctuations assumed to be Gaussian for the central limit theorem and uncorrelated in space and time: , where is the diffusion coefficient [18].
The longrange stress field exercised by all the other dislocations over the nth, the PeachKoehler force , is computed considering the image dislocations method to comply with periodic boundary conditions along the direction. Making use of calculations in [23, 24] one can find the following expression
(1) 
where is the Poisson’s ratio, and . Finally the overdamped Langevin equation [18] for the GB reads
(2) 
for , or rather
(3) 
To indicate the amplitude of the and forces we introduce respectively the parameters and .
The key quantities that we consider in order to characterize the dynamics of the GB are:

the meansquare displacement of the center of mass, where ;

the meansquare width .
In the following, we first analyze the case of a flat GB for which a comparison with MD simulations approach [17] is made. Next we consider the full flexible description of the GB. Finally, we discuss a linearized version of the model that can be treated analytically.
3 Flat grain boundary
For many applications a good approximation is to consider a flat GB with a single degree of freedom, for which and for . In other words, the flat GB is described by the following equation
(4) 
where the correlation properties of the thermal fluctuations are: and . This type of equation, also known as the Kramers equation with periodic potential, has been extensively studied in the literature [18]. In particular, the meansquare displacement is known to display a combination of oscillatory and diffusive behavior [18, 25]. Different dynamical regimes are found as the potential strength or the friction varies [25]. In fact, we show next that this simple model allows to understand the shorttime deviations from diffusive behavior observed in MD simulations [17].
Integrating Eq. 4 with the initial condition by means of computer simulations (fitting and with the condition [18] where is the mobility and is the renormalized diffusion coefficient considered in [17]), we have compared the meansquare displacement to the one obtained from MD simulation in Ref. [17]. In Fig. 1 this comparison is displayed together with the meansquare displacement of the renormalized free Brownian motion (described by the equation: with and ). The agreement between the two simulations is extremely good. For higher times (), the meansquare displacement tends to the renormalized Brownian motion. Hence taking explicitly into account the sinusoidal PeierlsNabarro force in the Langevin equation allows to describe the meansquare displacement for early times of the dynamics.
One can deduce in a simplified intuitive way the temporal evolution of the meansquare displacement starting by the transition probability density for small times (small ) [18]
(5) 
Next we consider the transition probability to run from the point at time to the point at time and to run from at to at
(6) 
For a free Brownian motion () the condition implies and stochastic displacements are space independent. If we impose this condition in presence of a periodic force , we obtain
(7) 
and then
(8) 
This result implies that, with the initial condition , if the potential is convex (concave) the meansquare displacement curve is concave (convex). In the case of the PN potential, we find indeed that the meansquare displacement curve should display upper and lower deviations from the straight line, corresponding to a renormalized free Brownian motion, depending in . These deviations decrease with time so that for large times the curve should approach a straight line [18].
4 Flexible grain boundary
A more general description of the GB considers its internal deformation and the dynamics is described by Eq. 3. The dynamical behavior of the GB depends on the amplitude of the three terms in the righthand side of Eq. 3.
The parameters that characterize the behavior of the GB are , , , and . Varying the values of these parameters, in the longtime limit the GB can either exfoliate (when the noise, , is high enough with respect to and to ) or reach a stationary state (when the noise, , is small when compared to or to ). The asymptotic behavior can be read off from the width that keeps on increasing when the GB exfoliates and saturates when the GB remains stable. In Fig. 2 the comparison between these two typical situations is displayed in the case of , , , and for the case in which the GB remains stable, while for the case in which the GB exfoliates.
In what follows, we analyze the dynamical behavior of the stable GB for , , and . In Fig. 3 the average position of the GB center of mass is displayed with and without the PN force in LogLog scale for .
For long times in both cases we have a linear behavior , but for short times, in the presence of the PN force, there is a clear deviation from linearity. This result confirms the conclusion made in the previous section, that the PN force is the cause for the deviation from linearity of for short times observed in Ref. [17]. Next we characterize the morphology of the GB through the width . In Fig. 4 for is displayed with and without the PN force in LogLog scale. In the absence of the PN force (Fig. 4a) the time dependence of is qualitatively similar to the same case but with linearized PK force discussed in the next section, while in the presence of the PN force, for (Fig. 4b), exhibit a plateau for intermediate times.
5 Continuum Theory
It is possible to develop an analytic expression in the continuum limit (, and ) for short or long times for in absence of the PN force linearizing the PK force. The equation of motion for and linearized is
(9) 
To obtain the short time behavior is sufficient to rewrite Eq. 9 as a generalized OrnsteinUhlenbeck process [18]
(10) 
The general solution of Eq. 10 is
(11) 
with (where ). From the definition of , results
(12) 
Replacing the Taylor expansion of the matrix in Eq. 12 one obtains for short times () that and in the continuum limit . To obtain the long times behavior of , we rewrite Eq. 9 in Fourier space [5, 26]. Employing the decomposition , we obtain
(13) 
The first term in the righthand side of Eq. 13 can be rewritten as
(14) 
where . Using the following results
(15) 
we obtain
(16) 
so that Eq. 13 becomes
(17) 
In the long time limit (large , small ) the term can be neglected. Finally in the continuum limit we replace , by the continuum variables , and .
To reproduce the continuum limit by simulations of Eq. 9 we would need a very large GB, with . Thus a comparison between Eq. 19 and the simulations for small is possible only by introducing some effective parameters in Eq. 19. In Fig. 5 computed by the simulations with (for which the better statistic is available) is compared with the fitted theoretical prediction for short and long times. Eq. 19 also predicts that the saturated value of width () exhibits a logarithmic dependence on the GB length : . In the case in which is not linearized this result is also confirmed by numerical simulations, showing that increases logarithmically with when (see Fig. 6). In presence of a periodic potential (), however, we observe a deviation from the logarithmic growth at large . This suggests that the PeierlsNabarro potential may set a limit to the GB roughness.
6 Summary and Discussion
We have investigated the diffusion of a regularly spaced lowangle grain boundary in a crystalline material. A typical computational method to describe the dynamics of the grain boundary is to perform deterministic molecular dynamics simulations with appropriate interatomic interactions [17]. Here we have employed the overdamped Langevin approach to obtain a long time description of the dynamics, but in particular to perform a comparison with molecular dynamics simulations for a specific material [17]. The first results is the interpretation of the early times behavior of the meansquare displacement . The deviation for early times of by the case of the renormalized Brownian motion, that holds for long times, can be interpreted as the effect on the dislocations of the periodicity of the lattice giving rise to the PeierlsNabarro potential. Secondly the description of the dynamic () and the morphology () of the grain boundary by means of overdamped Langevin equations is in qualitatively good agreement with its behavior in real materials, so this approach can be considered an useful tool for these studies.
References
References
 [1] Barabàsi A L and Stanley H E 1995 Fractal Concepts in Surface Growth (Cambridge University Press)
 [2] Krug J 1997 Adv. Phys. 46 139
 [3] Zaiser M 2006 Adv. Phys. 54 185
 [4] Zapperi S and Zaiser M 2001 Mat. Sci. and Eng. A 348 309
 [5] Moretti P, Miguel MC, Zaiser M, and Zapperi S 2004 Phys. Rev. B 69 214103.
 [6] Moretti P, Miguel MC and Zapperi S 2005 Phys. Rev. B 72 014505.
 [7] Kardar M 1998 Phys. Rep. 301 85
 [8] Lemerle S, Ferré J, Chappert C, Mathet V, Giamarchi T and Le Doussal P 1998 Phys. Rev. Lett. 80 849
 [9] Zapperi S, Cizeau P, Durin G and Stanley H E 1998 Phys. Rev. B 58 6353
 [10] Bhattacharya S and Higgins M J 1993 Phys. Rev. Lett. 70 2617
 [11] Surdeanu R, Wijngarden R J, Visser E, Huijbregtse J M, Rector J H, Dam B and Griessen R 1999 Phys. Rev. Lett. 83 2054
 [12] Schäffer E and Wong P z 2000 Phys. Rev. E 61 5257
 [13] Rolley E, Guthmann C, Gombrowicz R and Repain V 1998 Phys. Rev. Lett. 80 2865
 [14] Bouchaud E 1997 J. Phys. C 9 4319
 [15] Schmittbuhl J and Måløy K 1997 Phys. Rev. Lett. 78 3888
 [16] Sutton A P and Balluffi R W 1995 Interfaces in Crystalline Materials (Monographs on the Physics and Chemistry of Materials, Clarendon Press)
 [17] Trautt Z T, Upmanyu M and Karma A 2006 Science 314 632
 [18] Risken H 1984 The FokkerPlank Equation (Springer Verlag)
 [19] Ananthakrishna G 2007 Phys. Rep. 440 113259
 [20] Fulde P, Pietronero L, Schneider W R and Strässler S 1975 Phys. Rev. Lett. 35 1776
 [21] Dieterich W, Fulde P and Peschel I 1980 Adv. Phys. 29 527
 [22] Marchesoni F and Vij J K 1985 Z. Phys. B 58 187
 [23] Hirth J P and Lothe J 1982 Theory of Dislocations (Wiley & Sons)
 [24] Friedel J 1964 Dislocations (Pergamon Press)
 [25] Ferrando R, Spadacini R, Tommei G E and Caratti G 1992 Physica A 195 506532
 [26] Chui ST 1983 Phys. Rev. B 28 178
 [27] Krug J and Meakin P 1991 Phys. Rev. Lett. 66 703