PARADISCYL:Cylinder-Surface cross slip

Surface cross slip

ill Ryu and Wei Cai

This tutorial describes how to implement cross slip mechanism from the cylinder surface. The theoretical background is published in Computing Image Stress in an Elastic Cylinder(Proceedings of the National Academy of Sciences, 105, 14304 (2008)http://micro.stanford.edu/~caiwei/papers/Weinberger08PNAS-bccpillar.pdf (PDF)])

How to select slip plane of the surface segments

To implement surface cross slip in cylinder code, we change the slip plane of the surface nodes with respect to the magnitude of the force on the surface nodes.To do that,dislocation character of the surface nodes should be screw-like.

If there is a screw dislocation in the BCC crystal in left figure, we can see that there are three possible slip planes(red, blue green planes in the figure). From the MD calculation, we knew the image stress generate the force to move dislocation in the direction along which the length of dislocation shortens.Therefore, image force points downward for the front node, while it points upward for the back node(Figure 1). However, P-K force points in same direction both front node and back node. Taking the summation of these forces into account, the slip plane is selected as the one on which the summation of forces has maximum value.

(a)	(b)
Fig.1 (a) When the P-K force is dominant(blue plane is taken as a slip plane for the front node, while red plane is selected for the back node)(b) When the image force is dominant(green plane is taken as a slip plane) In this mechanism, straight edge dislocation has two different slip planes(See Figure 1-(a)). Algorithm 1. Searching for the surface node(node in Figure 2.(a)) 2. Find the neighbor node(nbr1 in Figure 2.(a)) 3. Check if the character of surface segment is similar to screw $1.0-\mathbf {b} \cdot \mathbf {\xi } <=\epsilon$ , where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon} is a tolerance. 4. Given burgers vector, there are three possible slip planes. For example, if Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{b} = a[111]} , then possible slip planes are Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{n}_1 = [1\bar{1}1],\mathbf{n}_2 = [0\bar{1}1],\mathbf{n}_3 = [\bar{1}01]} (See figure 2(b)) 5. Compute projected forces on each plane. Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{F}_{i} = \mathbf{F}_{total}-(\mathbf{F}_{total} \cdot \mathbf{n}_i)\mathbf{n}_i} 6. Based on the magnitudes of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{F}_{i}} , choose slip plane of the surface dislocation segment. (a) (b) Fig.3 (a)Schematic (b) Slip system Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{b} = a[111],\mathbf{n}_1 = [1\bar{1}1],\mathbf{n}_2 = [0\bar{1}1],\mathbf{n}_3 = [\bar{1}01]} .