nature Vol 440127 April 2006|doi:10.1038/nature04658 LETTERS Dislocation multi-junctions and strain hardening Vasily V.Bulatov,Luke L.Hsiung',Meijie Tang',Athanasios Arsenlis',Maria C.Bartelt'Wei Cai2 Jeff N.Florando',Masato Hiratani',Moon Rhee',Gregg Hommes',Tim G.Piercel&Tomas Diaz de la Rubial At the microscopic scale,the strength of a crystal derives from the Owing to their fixed ends,the lines merge only partially (when motion,multiplication and interaction of distinctive line defects b2=-b the lines will partially annihilate).Bounded at its ends by called dislocations.First proposed theoretically in 1934(refs 1-3) two triple nodes,the resulting junction zips along the [111]direction to explain low magnitudes of crystal strength observed experi- because each of the two parent dislocations is allowed to move only mentally,the existence of dislocations was confirmed two decades on its glide plane(the plane containing both the Burgers vector and later45.Much of the research in dislocation physics has since the dislocation line itself)with normal vectors n=(011)and focused on dislocation interactions and their role in strain hard- n2=(101),respectively.Because most of the interacting dislocations ening,a common phenomenon in which continued deformation move on non-parallel glide planes,attractive collisions zip junctions increases a crystal's strength.The existing theory relates strain of limited length along the intersection lines of the glide planes. hardening to pair-wise dislocation reactions in which two inter- The frequency and strength of such pair-wise dislocation reactions secting dislocations form junctions that tie the dislocations tying dislocations together are believed to control the physics of together7.Here we report that interactions among three dislo- strain hardening:a common phenomenon in which continued cations result in the formation of unusual elements of dislocation deformation increases a crystal's strength. network topology,termed 'multi-junctions'.We first predict The existence and the important role of dislocations junctions in the existence of multi-junctions using dislocation dynamics and strain hardening has been confirmed by numerous theoretical?and atomistic simulations and then confirm their existence by experimental studies and,more recently,by DD simulations transmission electron microscopy experiments in single-crystal Very recently,analysing previously publishedand our own new molybdenum.In large-scale dislocation dynamics simulations, simulations,we observed the formation of complex topological multi-junctions present very strong,nearly indestructible, connections in which more than two dislocation lines merge obstacles to dislocation motion and furnish new sources for together.Curious about possible causes for such anomalous for- dislocation multiplication,thereby playing an essential role in mations,we proceeded to investigate whether attractive reactions the evolution of dislocation microstructure and strength of among three or more dislocations are possible.For b.c.c.crystals,one deforming crystals".Simulation analyses conclude that multi- such candidate reaction is: junctions are responsible for the strong orientation dependence 1/2T11+1/2[1T)+1/2[11T=1/2[111] (1) of strain hardening in body-centred cubic crystals. ba The amount of slip produced by a propagating dislocation is quantified by its Burgers vector b,which is equal to one of the Given that the elastic energy stored in a dislocation's strain field (typically smallest)repeat vectors of the crystal lattice.Exactly what is proportional to llbl,the Frank estimate of the energy reduction happens when two dislocations collide depends on the lengths of two lines,their collision geometry and applied stress.Most significantly, the collision outcomes are affected by the Burgers vectors of two colliding dislocations.Given that a dislocation's energy is pro- portional to the square of its Burgers vector,the approximate Frank energy criterion predicts that when (b+b)2<b+b or, equivalently,bb2<0,the two lines will attract and merge into a product dislocation-a'junction'-with Burgers vector bi=b+b2, thereby reducing the internal energy of the system.In particular, when b2=-b,the two dislocations can annihilate completely, leaving no product. Figure la and b shows a junction-forming collision of two dislocation lines in a dislocation dynamics (DD)simulation (see the Methods section and Supplementary Discussion 1 for more details).The initial configuration consists of two straight dislocation lines of equal length made to intersect at their midpoints,while their Figure 1 Formation(zipping)and yielding of dislocation junctions in the endpoints are rigidly fixed (Fig.1a).Expressed in the units of the DD simulations.a,Two dislocation lines are initially brought to lattice constant,the Burgers vectors of two lines are b=1/2[111] intersection at their midpoints.b,Once the interaction between two lines is and b2=1/2[111],typical of the body-centred cubic(b.c.c.)crystals. turned on,two lines zip a binary junction,J.c,A snapshot showing a binary junction unzipping under stress.d,A third line is brought to intersect the In the DD simulation,the elastic interaction between two lines causes binary junction.e,The interaction among three lines makes them zip a them to merge into a junction dislocation with Burgers vector long multi-junction,4.f,A snapshot showing a multi-junction acting as a bi=b+b2=[001](Fig.1b,see also Supplementary Video 1). Frank-Read source of dislocation multiplication. Lawrence Livermore National Laboratory,University of California,Livermore,California94550,USA.2Department of Mechanical Engineering.Stanford University,Stanford, California 94305,USA. 1174 2006 Nature Publishing Group
© 2006 Nature Publishing Group Dislocation multi-junctions and strain hardening Vasily V. Bulatov1 , Luke L. Hsiung1 , Meijie Tang1 , Athanasios Arsenlis1 , Maria C. Bartelt1 , Wei Cai1,2, Jeff N. Florando1 , Masato Hiratani1 , Moon Rhee1 , Gregg Hommes1 , Tim G. Pierce1 & Tomas Diaz de la Rubia1 At the microscopic scale, the strength of a crystal derives from the motion, multiplication and interaction of distinctive line defects called dislocations. First proposed theoretically in 1934 (refs 1–3) to explain low magnitudes of crystal strength observed experimentally, the existence of dislocations was confirmed two decades later4,5. Much of the research in dislocation physics has since focused on dislocation interactions and their role in strain hardening, a common phenomenon in which continued deformation increases a crystal’s strength. The existing theory relates strain hardening to pair-wise dislocation reactions in which two intersecting dislocations form junctions that tie the dislocations together6,7. Here we report that interactions among three dislocations result in the formation of unusual elements of dislocation network topology, termed ‘multi-junctions’. We first predict the existence of multi-junctions using dislocation dynamics and atomistic simulations and then confirm their existence by transmission electron microscopy experiments in single-crystal molybdenum. In large-scale dislocation dynamics simulations, multi-junctions present very strong, nearly indestructible, obstacles to dislocation motion and furnish new sources for dislocation multiplication, thereby playing an essential role in the evolution of dislocation microstructure and strength of deforming crystals8 . Simulation analyses conclude that multijunctions are responsible for the strong orientation dependence of strain hardening in body-centred cubic crystals. The amount of slip produced by a propagating dislocation is quantified by its Burgers vector b, which is equal to one of the (typically smallest) repeat vectors of the crystal lattice. Exactly what happens when two dislocations collide depends on the lengths of two lines, their collision geometry and applied stress. Most significantly, the collision outcomes are affected by the Burgers vectors of two colliding dislocations. Given that a dislocation’s energy is proportional to the square of its Burgers vector, the approximate Frank energy criterion9 predicts that when (b1 þ b2) 2 , b1 2 þ b2 2 or, equivalently, b1b2 , 0, the two lines will attract and merge into a product dislocation—a ‘junction’—with Burgers vector bj ¼ b1 þ b2, thereby reducing the internal energy of the system. In particular, when b2 ¼ 2b1, the two dislocations can annihilate completely, leaving no product. Figure 1a and b shows a junction-forming collision of two dislocation lines in a dislocation dynamics (DD) simulation (see the Methods section and Supplementary Discussion 1 for more details). The initial configuration consists of two straight dislocation lines of equal length made to intersect at their midpoints, while their endpoints are rigidly fixed (Fig. 1a). Expressed in the units of the lattice constant, the Burgers vectors of two lines are b1 ¼ 1/2[111] and b2 ¼ 1/2[111], typical of the body-centred cubic (b.c.c.) crystals. In the DD simulation, the elastic interaction between two lines causes them to merge into a junction dislocation with Burgers vector bj ¼ b1 þ b2 ¼ [001] (Fig. 1b, see also Supplementary Video 1). Owing to their fixed ends, the lines merge only partially (when b2 ¼ 2b1 the lines will partially annihilate). Bounded at its ends by two triple nodes, the resulting junction zips along the [111] direction because each of the two parent dislocations is allowed to move only on its glide plane (the plane containing both the Burgers vector and the dislocation line itself) with normal vectors n1 ¼ (011) and n2 ¼ (101), respectively. Because most of the interacting dislocations move on non-parallel glide planes, attractive collisions zip junctions of limited length along the intersection lines of the glide planes. The frequency and strength of such pair-wise dislocation reactions tying dislocations together are believed to control the physics of strain hardening: a common phenomenon in which continued deformation increases a crystal’s strength. The existence and the important role of dislocations junctions in strain hardening has been confirmed by numerous theoretical6,7 and experimental10,11 studies and, more recently, by DD simulations12,13. Very recently, analysing previously published14 and our own new simulations, we observed the formation of complex topological connections in which more than two dislocation lines merge together. Curious about possible causes for such anomalous formations, we proceeded to investigate whether attractive reactions among three or more dislocations are possible. For b.c.c. crystals, one such candidate reaction is: b1 1=2½111 þ b2 1=2½111 þ b3 1=2½111 ¼ b4 1=2½111 ð1Þ Given that the elastic energy stored in a dislocation’s strain field is proportional to kbk2 , the Frank estimate of the energy reduction LETTERS Figure 1 | Formation (zipping) and yielding of dislocation junctions in the DD simulations. a, Two dislocation lines are initially brought to intersection at their midpoints. b, Once the interaction between two lines is turned on, two lines zip a binary junction, J. c, A snapshot showing a binary junction unzipping under stress. d, A third line is brought to intersect the binary junction. e, The interaction among three lines makes them zip a long multi-junction, 4. f, A snapshot showing a multi-junction acting as a Frank–Read source of dislocation multiplication. 1 Lawrence Livermore National Laboratory, University of California, Livermore, California 94550, USA. 2 Department of Mechanical Engineering, Stanford University, Stanford, California 94305, USA. Vol 440|27 April 2006|doi:10.1038/nature04658 1174
NATUREVol 44027 April 2006 LETTERS in the above reaction,(Ilbl2+b2ll+llbll2)/lbill2 =3,is much junction depends on multiple factors,including dislocation line greater than in the binary reaction shown in Fig.1b,(lbll2+b2l2)/ lengths and orientations,direction of applied stress and exactly llbill2=1.5.To see whether such a reaction is indeed feasible,we now how the junction is incorporated in the dislocation network.To overlay on top of the binary junction a third dislocation(Fig.Id) assess qualitatively the relative holding power of binary and multi- with Burgers vector ba=1/2[111]and glide plane na=(110).The junctions,we performed a large series of DD simulations in which DD simulation result is shown in Fig.le,where the third dislocation both binary and multi-junctions are formed in the same geometries with Burgers vector ba has reacted with the junction dislocation and and subjected to the same straining conditions(see Supplementary transformed it into a 1/2(111 dislocation,exactly in accord with the Discussion 1 for the detailed results).Depending on the Burgers proposed reaction (1)(Supplementary Video 2).The transformed vectors and line orientations,the dislocations may mutually repel or junction segment with Burgers vector ba of 1/2-type connects attract each other.In cases when they attract,the lines either zip together three dislocations at its ends and is defined as a multi- junctions or stay crossed'.The crossed attractive configurations do junction.Remarkably,the resulting multi-junction extends well not hold dislocations together appreciably and are destroyed by beyond the original length of the binary junction,corroborating application of a small stress,whereas junction unzipping(Fig.1c) our expectation that the ternary reaction (1)can result in a much necessitates a significantly higher stress. greater energy reduction than the classical binary reaction. The superior strength of multi-junctions compared to binary The nodes at each end of the multi-junction that tie together four junctions manifests itself in several ways.First,the multi-junctions dislocation lines we term 'multi-nodes'or'4-nodes.These 4-nodes form and exhibit measurable strength over a wider range of initial are distinct from simple crosses of two dislocations in that all four line orientations than the binary junctions.Second,in the collision lines entering the node have different Burgers vectors.These 4-nodes geometries where both binary and multi-junctions zip,the latter are beautifully symmetric:all four distinct Burgers vectors of require much higher stress to unzip and release the dislocations. 1/2-type enter the 4-node exactly once.Therefore,it is the Finally,whereas the binary junctions could eventually be unzipped only possible 4-node of this kind in b.c.c.crystals.Curiously,the under all tested line and stress orientations,the multi-junctions were same symmetric 4-node can be formed through four different reac- found to be indestructible across a wide range ofline orientations and tions among three lines,for example b +b2+(-b)=(-ba).This stress states.In such cases,rather than unzip and release the lines to non-uniqueness brings about an interesting feature of the dislocation network topology in b.c.c.crystals that is not present in the conventional dislocation networks consisting solely of binary junctions.In the latter,it is possible to trace every single line with 1/2Burgers vector through each junction it enters.It is even possible to uniquely deconstruct the entire network into individual 1/2lines.However,the topology of dislocation networks containing 4-nodes is different in principle:it is now impossible to specify which of the four dislocations in a particular 4-node are the parents and which is the product,and to 'trace'a given 1/2line through the network.While it is still possible to deconstruct the network into constituent lines,there are a combinatorially large number of different ways to do this.The formation of multi-junctions results in the topological irreversibility or untraceability of dislocation networks. Although the DD simulations provide credible support for the existence of multi-junctions,it is desirable to verify this finding with a discrete atomistic model that does not rely on the continuum theory of dislocations.Figure 2a shows the result of one such simulation in which three different dislocations (left-hand side of equation (1))were introduced into a small fragment of the b.c.c. single crystal and then allowed to relax the lattice distortions produced by the inserted dislocations(see the Methods section for details).In the relaxed configuration (Fig.2a),two distinct 4-nodes are instantly recognizable,as is the junction dislocation with the 1/2[111]Burgers vector (right-hand side of equation(1)).Taken together,the DD and the atomistic simulations appear to be convincing:multi-junctions should exist. Yet,for definitive verification,we rely on transmission electron microscopy (TEM)of molybdenum single crystals.Figure 2b-d shows three different views of a single fragment of the dislocation network containing a binary junction and a 4-node.Four lines entering the 4-node are numbered from 1 to 4.A unique TEM signature of the symmetric 4-node is that,in certain contrast con- ditions,one of four dislocations must appear out of contrast while the other three lines remain visible.The appearance ofa symmetric 4-node 0.2Hm in each of the four TEM frames shown in Fig.2c and d is unmistakable. Similar 4-node dislocation configurations were also found in other Figure 2 Atomistic simulations and experiments confirm that multi- regions of the TEM foil,leading us to believe that their occurrence junctions exist in b.c.c.molybdenum.a,A multi-junction formed in an may not be rare. atomistic simulation.b,A TEM micrograph containing a symmetric 4-node As discussed above,multi-junctions appear to hold dislocations In this view all four dislocations(1-4)entering the multi-node are visible. together more tightly than binary junctions.It is difficult to quantify c,View in which dislocation 1 becomes invisible.The length of the scale bar this difference exactly because the stress required to unzip a given is 0.2 um.d,Another view in which dislocation 3 is invisible. 1175 2006 Nature Publishing Group
© 2006 Nature Publishing Group in the above reaction, (kb1k2 þ kb2k2 þ kb3k2 )/kb4k2 ¼ 3, is much greater than in the binary reaction shown in Fig. 1b, (kb1k2 þ kb2k2 )/ kbj k2 ¼ 1.5. To see whether such a reaction is indeed feasible, we now overlay on top of the binary junction a third dislocation (Fig. 1d) with Burgers vector b3 ¼ 1/2[111] and glide plane n3 ¼ (110). The DD simulation result is shown in Fig. 1e, where the third dislocation with Burgers vector b3 has reacted with the junction dislocation and transformed it into a 1/2[111] dislocation, exactly in accord with the proposed reaction (1) (Supplementary Video 2). The transformed junction segment with Burgers vector b4 of 1/2,111.-type connects together three dislocations at its ends and is defined as a multijunction. Remarkably, the resulting multi-junction extends well beyond the original length of the binary junction, corroborating our expectation that the ternary reaction (1) can result in a much greater energy reduction than the classical binary reaction. The nodes at each end of the multi-junction that tie together four dislocation lines we term ‘multi-nodes’ or ‘4-nodes’. These 4-nodes are distinct from simple crosses of two dislocations in that all four lines entering the node have different Burgers vectors. These 4-nodes are beautifully symmetric: all four distinct Burgers vectors of 1/2,111.-type enter the 4-node exactly once. Therefore, it is the only possible 4-node of this kind in b.c.c. crystals. Curiously, the same symmetric 4-node can be formed through four different reactions among three lines, for example b1 þ b2 þ (2b4) ¼ (2b3). This non-uniqueness brings about an interesting feature of the dislocation network topology in b.c.c. crystals that is not present in the conventional dislocation networks consisting solely of binary junctions. In the latter, it is possible to trace every single line with 1/2,111. Burgers vector through each ,001. junction it enters. It is even possible to uniquely deconstruct the entire network into individual 1/2,111. lines. However, the topology of dislocation networks containing 4-nodes is different in principle: it is now impossible to specify which of the four dislocations in a particular 4-node are the parents and which is the product, and to ‘trace’ a given 1/2,111. line through the network. While it is still possible to deconstruct the network into constituent lines, there are a combinatorially large number of different ways to do this. The formation of multi-junctions results in the topological irreversibility or untraceability of dislocation networks. Although the DD simulations provide credible support for the existence of multi-junctions, it is desirable to verify this finding with a discrete atomistic model that does not rely on the continuum theory of dislocations. Figure 2a shows the result of one such simulation in which three different dislocations (left-hand side of equation (1)) were introduced into a small fragment of the b.c.c. single crystal and then allowed to relax the lattice distortions produced by the inserted dislocations (see the Methods section for details). In the relaxed configuration (Fig. 2a), two distinct 4-nodes are instantly recognizable, as is the junction dislocation with the 1/2[111] Burgers vector (right-hand side of equation (1)). Taken together, the DD and the atomistic simulations appear to be convincing: multi-junctions should exist. Yet, for definitive verification, we rely on transmission electron microscopy (TEM) of molybdenum single crystals. Figure 2b–d shows three different views of a single fragment of the dislocation network containing a binary junction and a 4-node. Four lines entering the 4-node are numbered from 1 to 4. A unique TEM signature of the symmetric 4-node is that, in certain contrast conditions, one of four dislocations must appear out of contrast while the other three lines remain visible. The appearance of a symmetric 4-node in each of the four TEM frames shown in Fig. 2c and d is unmistakable. Similar 4-node dislocation configurations were also found in other regions of the TEM foil, leading us to believe that their occurrence may not be rare. As discussed above, multi-junctions appear to hold dislocations together more tightly than binary junctions. It is difficult to quantify this difference exactly because the stress required to unzip a given junction depends on multiple factors, including dislocation line lengths and orientations, direction of applied stress and exactly how the junction is incorporated in the dislocation network. To assess qualitatively the relative holding power of binary and multijunctions, we performed a large series of DD simulations in which both binary and multi-junctions are formed in the same geometries and subjected to the same straining conditions (see Supplementary Discussion 1 for the detailed results). Depending on the Burgers vectors and line orientations, the dislocations may mutually repel or attract each other. In cases when they attract, the lines either zip junctions or stay crossed12. The crossed attractive configurations do not hold dislocations together appreciably and are destroyed by application of a small stress, whereas junction unzipping (Fig. 1c) necessitates a significantly higher stress. The superior strength of multi-junctions compared to binary junctions manifests itself in several ways. First, the multi-junctions form and exhibit measurable strength over a wider range of initial line orientations than the binary junctions. Second, in the collision geometries where both binary and multi-junctions zip, the latter require much higher stress to unzip and release the dislocations. Finally, whereas the binary junctions could eventually be unzipped under all tested line and stress orientations, the multi-junctions were found to be indestructible across a wide range of line orientations and stress states. In such cases, rather than unzip and release the lines to Figure 2 | Atomistic simulations and experiments confirm that multijunctions exist in b.c.c. molybdenum. a, A multi-junction formed in an atomistic simulation. b, A TEM micrograph containing a symmetric 4-node. In this view all four dislocations (1–4) entering the multi-node are visible. c, View in which dislocation 1 becomes invisible. The length of the scale bar is 0.2 mm. d, Another view in which dislocation 3 is invisible. NATURE|Vol 440|27 April 2006 LETTERS 1175
LETTERS NATUREVol 440127 April 2006 glide under increasing stress(Supplementary Video 3),the multi- simulation of uniaxial straining along (001]direction exhibits a high junction yields by repetitively emitting concentric dislocation loops rate of strain hardening,as shown by a pronounced slope of the and returning to its zipped configuration(Fig.Ifand Supplementary stress-strain curve beyond the initial yield,whereas there is virtually Video 4),thus forming a new regenerative dislocation source of the no hardening exhibited in the [011]case (Fig.3a).Because dislo- Frank-Read type's.In contrast,under no stress state was a binary cations of at least three distinct Burgers vectors must be present to junction observed to act as a regenerative dislocation source. form multi-junctions,they are expected to form more frequently in We now attempt to determine whether multi-junctions play a role the [001]straining simulation where all four distinct Burgers vectors in strain hardening,that is,when continued straining demands of 1/2type are active and not in the [011]case where only two increasingly higher stress.The small-scale DD simulations reveal of the four are active.Observed both in simulations and in experi- the great strength of multi-junctions and their propensity to act as ment,the strong contrast in the orientation dependence of strain regenerative sources and imply that these dislocation tangles may hardening makes it tempting to assert that the difference is somehow play an important part in strain hardening,but alone are incon- related to multi-junctions.This assertion is further supported by the clusive.At the same time,transmission electron micrographs have diminishing multiplication rate observed in the 011]case compared demonstrated the existence of multi-junctions in b.c.c.molybdenum to the rapid and steady dislocation multiplication observed in the and support a connection between multi-junctions and strain hard- [001]case (Fig.3b),and explained by the propensity of multi- ening,but because they were taken several days after the straining junctions to form additional dislocation sources. experiments,are not conclusive.Large-scale DD simulations present To test further the assertion that multi-junctions strongly influ- the opportunity to observe the formation of multi-junctions during ence the orientation dependence of strain hardening,we now per- straining and to quantify their effects on strain hardening,via in situ form a computational experiment that is impossible to reproduce in computational experiments. the laboratory.Specifically,we repeat the same [001]straining The extreme computational cost of DD simulations has made it simulation but start with a modified initial configuration in which unfeasible to simulate dislocation ensembles large enough to com- the Burgers vectors b;and b are converted into b and bz,respect- pute strain hardening directly from the underlying motion and ively,such that the total density of dislocations remains initially interactions of many dislocation lines.Recently,we developed a unchanged.By 'doctoring'the initial structure in such a way,we new DD code ParaDiS(for Parallel Dislocation Simulator)specifi- preclude any possibility of multi-junction formation during the cally designed to take full advantage of massively parallel super- course of this simulation.The resulting stress-strain and stress- computers'5.ParaDis runs efficiently on two of the world's most density behaviours differ markedly from the full [001]straining powerful computers,Thunder and Blue Gene/L,both at the Univer- simulation (with all four Burgers vectors included).Instead,the sity of California Lawrence Livermore National Laboratory (http:/ behaviour is similar to that observed in the [011]straining simu- www.llnl.gov/computing/hpc/resources/OCF_resources.html).Here lation (Fig.3a and b).This observation further reinforces our ParaDiS enables direct simulations of large strains and strain assertion that the high hardening rate in straining along the high hardening in statistically representative dislocation ensembles. symmetry (for example,[001])directions is related to the formation Figure 3 shows the data obtained in a series of three DD of the multi-junctions(see Supplementary Discussion 2). simulations of single crystal molybdenum subjected to uniaxial To clarify how the multi-junctions affect strain hardening,we compression along two different directions (see also Supplementary investigate the relationship between the evolving dislocation micro- Fig.1).The simulation parameters were chosen to represent molyb- structure and the instantaneous flow stress.As shown in Fig.3c,over denum single crystals at an elevated temperature of 450 K.In the range of dislocation densities common to all three simulations, accordance with experimentss(see also Supplementary Fig.2),the the flow stress appears to be determined by the total dislocation 3×1013 200 150 2×1013 (edw) 100 101 0 0.01 0.02 0.01 0.02 Strain Strain 300 (edw) 0.3 d 200 0.2 100 0.1 1013 2×1013 3×1013 1013 2×1013 Dislocation density (m-) Dislocation density (m-2) Figure 3|The results of virtual straining experiments on b.c.c. Burgers vectors are absent.a,Flow stress as a function of strain. molybdenum.The black lines correspond to the full [001]straining b,Dislocation line density as a function of strain.c,Flow stress versus total simulation,the red lines are for the [011]straining,and the green lines are dislocation density.d,The fraction of lines involved in multi-junction for the 'doctored'[001]straining simulation in which two of the four configurations as a function of the total line density. 1176 @2006 Nature Publishing Group
© 2006 Nature Publishing Group glide under increasing stress (Supplementary Video 3), the multijunction yields by repetitively emitting concentric dislocation loops and returning to its zipped configuration (Fig. 1f and Supplementary Video 4), thus forming a new regenerative dislocation source of the Frank–Read type15. In contrast, under no stress state was a binary junction observed to act as a regenerative dislocation source. We now attempt to determine whether multi-junctions play a role in strain hardening, that is, when continued straining demands increasingly higher stress. The small-scale DD simulations reveal the great strength of multi-junctions and their propensity to act as regenerative sources and imply that these dislocation tangles may play an important part in strain hardening, but alone are inconclusive. At the same time, transmission electron micrographs have demonstrated the existence of multi-junctions in b.c.c. molybdenum and support a connection between multi-junctions and strain hardening, but because they were taken several days after the straining experiments, are not conclusive. Large-scale DD simulations present the opportunity to observe the formation of multi-junctions during straining and to quantify their effects on strain hardening, via in situ computational experiments. The extreme computational cost of DD simulations has made it unfeasible to simulate dislocation ensembles large enough to compute strain hardening directly from the underlying motion and interactions of many dislocation lines. Recently, we developed a new DD code ParaDiS (for Parallel Dislocation Simulator) specifi- cally designed to take full advantage of massively parallel supercomputers16. ParaDiS runs efficiently on two of the world’s most powerful computers, Thunder and Blue Gene/L, both at the University of California Lawrence Livermore National Laboratory (http:// www.llnl.gov/computing/hpc/resources/OCF_resources.html). Here ParaDiS enables direct simulations of large strains and strain hardening in statistically representative dislocation ensembles. Figure 3 shows the data obtained in a series of three DD simulations of single crystal molybdenum subjected to uniaxial compression along two different directions (see also Supplementary Fig. 1). The simulation parameters were chosen to represent molybdenum single crystals at an elevated temperature of 450 K. In accordance with experiments8 (see also Supplementary Fig. 2), the simulation of uniaxial straining along [001] direction exhibits a high rate of strain hardening, as shown by a pronounced slope of the stress–strain curve beyond the initial yield, whereas there is virtually no hardening exhibited in the [011] case (Fig. 3a). Because dislocations of at least three distinct Burgers vectors must be present to form multi-junctions, they are expected to form more frequently in the [001] straining simulation where all four distinct Burgers vectors of 1/2,111. type are active and not in the [011] case where only two of the four are active. Observed both in simulations and in experiment, the strong contrast in the orientation dependence of strain hardening makes it tempting to assert that the difference is somehow related to multi-junctions. This assertion is further supported by the diminishing multiplication rate observed in the [011] case compared to the rapid and steady dislocation multiplication observed in the [001] case (Fig. 3b), and explained by the propensity of multijunctions to form additional dislocation sources. To test further the assertion that multi-junctions strongly influence the orientation dependence of strain hardening, we now perform a computational experiment that is impossible to reproduce in the laboratory. Specifically, we repeat the same [001] straining simulation but start with a modified initial configuration in which the Burgers vectors b3 and b4 are converted into b1 and b2, respectively, such that the total density of dislocations remains initially unchanged. By ‘doctoring’ the initial structure in such a way, we preclude any possibility of multi-junction formation during the course of this simulation. The resulting stress–strain and stress– density behaviours differ markedly from the full [001] straining simulation (with all four Burgers vectors included). Instead, the behaviour is similar to that observed in the [011] straining simulation (Fig. 3a and b). This observation further reinforces our assertion that the high hardening rate in straining along the high symmetry (for example, [001]) directions is related to the formation of the multi-junctions (see Supplementary Discussion 2). To clarify how the multi-junctions affect strain hardening, we investigate the relationship between the evolving dislocation microstructure and the instantaneous flow stress. As shown in Fig. 3c, over the range of dislocation densities common to all three simulations, the flow stress appears to be determined by the total dislocation Figure 3 | The results of virtual straining experiments on b.c.c. molybdenum. The black lines correspond to the full [001] straining simulation, the red lines are for the [011] straining, and the green lines are for the ‘doctored’ [001] straining simulation in which two of the four Burgers vectors are absent. a, Flow stress as a function of strain. b, Dislocation line density as a function of strain. c, Flow stress versus total dislocation density. d, The fraction of lines involved in multi-junction configurations as a function of the total line density. LETTERS NATURE|Vol 440|27 April 2006 1176
NATUREVol 44027 April 2006 LETTERS Figure 4|Snapshots of dislocation network evolution obtained from a DD yield strain(0.2%),dislocations multiply and their collisions produce the simulation of [001]straining.a,Initially,dislocation motion results in first few multi-junctions (white).c,d,Continued dislocation multiplication binary collisions only so that the network remains all green (see Methods results in increasingly frequent dislocation collisions leading to strain section for a more detailed description of the colour scheme).b,Near the hardening and growth of the (white)sub-network of multi-junctions. density alone,irrespective of the orientation of the tensile axis and It remains to be seen what role multi-junctions play in the intricate the number of active Burgers vectors,suggesting that the weaker but dislocation patterns formed during straining of high-symmetry more numerous binary interactions predominantly define the overall crystals:it seems that multi-junctions can serve as strong anchors plastic strength.At the same time,the frequency of binary collisions is for dislocation tangles,braids,walls,etc.As to whether multi- determined by the dislocation density that increases much more junctions form in crystals other than b.c.c.,we predict that a variety rapidly in the high-symmetry,[001]straining due to proliferation of of strong multi-junctions should exist in the face-centred-cubic and new dislocation sources,multi-junctions. related high-symmetry crystals.Finally,theoretical analysis,DD As shown in Fig.3d,the topological composition of the dislocation simulations and TEM observations all suggest that dislocation microstructure differs significantly for different straining directions tangles even more complex than ternary junctions exist,but are at the same level of total dislocation density(see also Supplementary rare and their stability is likely to be marginal. Fig.3).Whereas in the [011]straining simulation both the total dislocation density and the fraction of lines involved in the multi- METHODS junctions saturate,the [001]straining simulation with the higher Dislocation dynamics simulations.In a DD simulation,dislocations are fraction of multi-junction configurations continues to evolve to represented by piece-wise straight segments interacting with each other accord- considerably higher total densities and higher fractions of lines ing to the equations of the continuum elasticity theory In asimulation,each involved in the multi-junctions.Taken together,these observations dislocation segment moves with velocity v=M.fproportional to the net force f indicate that the number of active Burgers vectors affects the rate of exerted on the segment by external loads and all other dislocation segments(here formation of new dislocation sources (multi-junctions),leading to M is the mobility tensor).A single simulation step includes:(1)calculating the significant differences in the accumulation of dislocation density forces acting on the segments,(2)advancing the segments to new positions according to their velocities,and(3)performing changes in the line topology and,thus,in the strain hardening rates (connectivity)when collisions or node instabilities are encountered.The force In situ visual observations reveal that,even in the [001]straining on a dislocation segment is calculated as the negative derivative of the system's simulations,the binary junctions form much more frequently than the energy with respect to the segment's position.The elastic constants and multi-junctions.In fact,for as long as dislocations move on different dislocation mobility function were chosen to capture the behaviour of b.c.c. planes they must intersect,making formation of binary junctions molybdenum above its athermal threshold (~450 K).The mobility was inde- unavoidable.In contrast,multi-junctions form infrequently,mostly pendent of the local line direction in the glide plane:that is,the plane containing by attachment of a third line to an existing binary junction.At the both line direction and Burgers vector.Mobility in the direction normal to the same time,while the binary junctions are observed to dynamically glide plane was a small fraction(10)of the glide mobility.Screw dislocations, unzip and reform elsewhere,the multi-junctions,once formed,are whose line directions are parallel to their Burgers vectors,were free to glide in any plane containing their line direction.All small-scale DD simulations were observed to endure.Zipping of new multi-junctions takes place conducted as though the configurations were in an infinite medium,whereas preferentially near the existing ones(see Supplementary Video 5), the larger-scale straining simulations were performed in a periodic cube 5 um on gradually building up a sub-network of multi-junctions,a strong and the side. mostly static backbone of the growing microstructure(Fig.4). Analysis of dislocation network topology.The entire network is comprised of 1177 2006 Nature Publishing Group
© 2006 Nature Publishing Group density alone, irrespective of the orientation of the tensile axis and the number of active Burgers vectors, suggesting that the weaker but more numerous binary interactions predominantly define the overall plastic strength. At the same time, the frequency of binary collisions is determined by the dislocation density that increases much more rapidly in the high-symmetry, [001] straining due to proliferation of new dislocation sources, multi-junctions. As shown in Fig. 3d, the topological composition of the dislocation microstructure differs significantly for different straining directions at the same level of total dislocation density (see also Supplementary Fig. 3). Whereas in the [011] straining simulation both the total dislocation density and the fraction of lines involved in the multijunctions saturate, the [001] straining simulation with the higher fraction of multi-junction configurations continues to evolve to considerably higher total densities and higher fractions of lines involved in the multi-junctions. Taken together, these observations indicate that the number of active Burgers vectors affects the rate of formation of new dislocation sources (multi-junctions), leading to significant differences in the accumulation of dislocation density and, thus, in the strain hardening rates. In situ visual observations reveal that, even in the [001] straining simulations, the binary junctions form much more frequently than the multi-junctions. In fact, for as long as dislocations move on different planes they must intersect, making formation of binary junctions unavoidable. In contrast, multi-junctions form infrequently, mostly by attachment of a third line to an existing binary junction. At the same time, while the binary junctions are observed to dynamically unzip and reform elsewhere, the multi-junctions, once formed, are observed to endure. Zipping of new multi-junctions takes place preferentially near the existing ones (see Supplementary Video 5), gradually building up a sub-network of multi-junctions, a strong and mostly static backbone of the growing microstructure (Fig. 4). It remains to be seen what role multi-junctions play in the intricate dislocation patterns formed during straining of high-symmetry crystals17: it seems that multi-junctions can serve as strong anchors for dislocation tangles, braids, walls, etc. As to whether multijunctions form in crystals other than b.c.c., we predict that a variety of strong multi-junctions should exist in the face-centred-cubic and related high-symmetry crystals. Finally, theoretical analysis, DD simulations and TEM observations all suggest that dislocation tangles even more complex than ternary junctions exist, but are rare and their stability is likely to be marginal. METHODS Dislocation dynamics simulations. In a DD simulation, dislocations are represented by piece-wise straight segments interacting with each other according to the equations of the continuum elasticity theory18,19. In a simulation, each dislocation segment moves with velocity v ¼ Mzf proportional to the net force f exerted on the segment by external loads and all other dislocation segments (here M is the mobility tensor). A single simulation step includes: (1) calculating the forces acting on the segments, (2) advancing the segments to new positions according to their velocities, and (3) performing changes in the line topology (connectivity) when collisions or node instabilities are encountered. The force on a dislocation segment is calculated as the negative derivative of the system’s energy with respect to the segment’s position. The elastic constants and dislocation mobility function were chosen to capture the behaviour of b.c.c. molybdenum above its athermal threshold (,450 K). The mobility was independent of the local line direction in the glide plane: that is, the plane containing both line direction and Burgers vector. Mobility in the direction normal to the glide plane was a small fraction (1026 ) of the glide mobility. Screw dislocations, whose line directions are parallel to their Burgers vectors, were free to glide in any plane containing their line direction. All small-scale DD simulations were conducted as though the configurations were in an infinite medium, whereas the larger-scale straining simulations were performed in a periodic cube 5 mm on the side. Analysis of dislocation network topology. The entire network is comprised of Figure 4 | Snapshots of dislocation network evolution obtained from a DD simulation of [001] straining. a, Initially, dislocation motion results in binary collisions only so that the network remains all green (see Methods section for a more detailed description of the colour scheme). b, Near the yield strain (,0.2%), dislocations multiply and their collisions produce the first few multi-junctions (white). c, d, Continued dislocation multiplication results in increasingly frequent dislocation collisions leading to strain hardening and growth of the (white) sub-network of multi-junctions. NATURE|Vol 440|27 April 2006 LETTERS 1177
LETTERS NATUREVol 440127 April 2006 nodes and links.A node is where three or more lines merge together,and a link is 5.Hirth,J.P.Lothe,J.Theory of Dislocations 2nd edn,9-13 (Wiley,New York, any line segment connecting two nodes of the network.For the present analysis, 1982). we label any 3-node that bounds a regular binary junction as a 'normal'or 6 Saada,G.Sur le durcissement do a recombinaison des dislocations.Acta N-node.Likewise,any 4-node formed by two dislocations crossing each other Metoll..8.841-852(1960). is also labelled an N-node.All other nodes are regarded as multi-nodes or 心 Franciosi,P.Berveiller,M.Zaoui,A.Latent hardening in copper and aluminium single crystals.Acto Metall.28,273-283 (1980). M-nodes,including the symmetric 4-nodes shown in Fig.le as well as 3-nodes 甲 Guiu,F.Pratt,P.L.The effect of orientation on the yielding and flow of produced by dissociation of the symmetric 4-nodes.Three types of links can now molybdenum single crystals.Phys.Status Solidi 15,539-552 (1966). be defined with respect to the types of nodes they connect:NN links,NM links 9 Hirth,J.P.Lothe,J.Theory of Dislocations 2nd edn,269-270 (Wiley,New and MM links.Here,notation NM is used for any link that connects an N-node York.1982). to an M-node.To compute the fraction of lines involved in the multi-junctions 10.Basinski,S.J.Basinski,Z.S.in Dislocations in Solids Vol.4 (ed.Nabarro, shown in Fig.3d,we summed the lengths of all MM links with half the lengths of F.R.N.)261-362,(North Holland,Amsterdam,1979). all NM links and divided this sum by the total length of all links.The colour Nabarro,F.R.N.,Basinski,Z.S.Holt,D.B.The plasticity of pure single scheme used in Fig.4 is as follows:all MM links are shown in white,whereas the crystals.Adk.Phy5.50,193-233(1964). 12.Madec.S.R..Devincre.B.Kubin.L.P.From dislocation junctions to forest colour of NM links is graded from white at the M-nodes to green at the N-nodes. hardening.Phys.Rev.Lett.89,255508 (2002). All NN links,including binary junctions,are shown in green. Atomistic simulation.The simulation volume was a small cube-shaped block of 13. Hahner,P.Zaiser,M.Dislocation dynamics and work hardening of fractal dislocation cell structures.Mater.Sci.Eng.A 272,443-454 (1999). a perfect b.c.c.single crystal,17 nm on each side.The initial geometry contained 14. Abraham,F.F.et al.Simulating materials failure by using up to one billion three dislocations with Burgers vectors 1/2[111],1/2[111],and 1/2[111]inter- atoms and the world's fastest computer:Work hardening.Proc.Natl Acad.Sci. secting at the block centre.The atom positions inside the block were then relaxed U5A99.5783-5787(2002). to mechanical equilibrium using the conjugate gradient method and an 15. Frank,F.C.Read,W.T.Multiplication processes for slow moving interatomic interaction function for molybdenum2.The atoms on the block dislocations.Phys.Rev.79,722-723 (1950). surfaces were fixed throughout the simulation.To visualize crystal defects,only 16. Bulatov,V.V.et al.Scalable line dynamics in ParaDiS.Supercomputing (http:// the atoms inside the block with energies exceeding the ideal bulk value by www.sc-conference.org/sc2004/schedule/pdfs/pap206.pdf)(2004). 0.095eV are shown. 17. Mughrabi,H.Dislocation wall and cell structures and long-range internal Experiment.The experiments involved three steps:(1)compression of a single- stresses in deformed metal crystals.Acta Metall.31,1367-1379 (1983). 18. Kubin,L P.Canova,G.R.The modelling of dislocation patterns.Scr.Metall. crystal molybdenum specimen to 1%total strain along the [001]axis,(2)cutting Mater..27,957-962(1992) and thinning the deformed specimen along the (101)plane to obtain electron 19. Ghoniem,N.M,Huang J.M.Wang,Z.Q.Affine covariant-contravariant transparent foils,and (3)TEM observations using a set of reflection vectors g vector forms for the elastic field of parametric dislocations in isotropic crystals that can reveal multi-junctions.In the view shown in Fig.2b the zone axis Phl.Mag.Let.82.55-63(2002). [101]and the diffraction vector g=[020],making all four dislocations 20.Finnis,M.W.Sinclair,J.5.A simple empirical N-body potential for transition entering the 4-node visible.The views in Fig.2c and d were obtained using metals.Phil.Mag.A 50,45-55 (1984). g=[121]and g=[121]which made lines b=1/2[111]and bs=1/2[111] invisible owing to the g-b =0 condition.To access additional diffraction vectors, Supplementary Information is linked to the online version of the paper at the specimen was tilted to a new zone axis (=[201])making it possible to www.nature.com/nature. identify the Burgers vectors of two remaining dislocations in a similar manner: Acknowledgements This work was supported by the US DOE Office of Basic b2 1/2[111]and ba 1/2[111],respectively. Energy Sciences and the NNSA ASC program.We thank E.Chandler and C.Mailhiot for their encouragement and unwavering support of the ParaDiS Received 14 November 2005;accepted 14 February 2006. code development effort,D.Lassila for providing the purified molybdenum crystals,M.LeBlanc for performing the mechanical tests,R.Cook for graphic 1.Orowan,E.Zur kristallplastizitat.Z.Phys.89,605-659 (1934). design and F.Abraham and G.Campbell for critical reading and editorial 2 Taylor,G.The mechanism of plastic deformation in crystal.Part I.Theoretical. suggestions for the manuscript.This work was performed under the auspices of Pr0c.R5oc.A145.363-404(1934). the US DOE by the Lawrence Livermore National Laboratory. 3.Polanyi,M.Uber eine Art Gutterstorung,die einen kristall plastich machen konnte.Z.Phys.89,660-664(1934). Author Information Reprints and permissions information is available at 4.Hirsch,P.B.,Homne,R.W.Whelan,M.J.Direct observations of the npg.nature.com/reprintsandpermissions.The authors declare no competing arrangement and motion of dislocations in aluminium.Phil.Mag.1,677-684 financial interests.Correspondence and requests for materials should be (1956). addressed to V.V.B.(bulatov1@llnl.gov). 1178 2006 Nature Publishing Group
© 2006 Nature Publishing Group nodes and links. A node is where three or more lines merge together, and a link is any line segment connecting two nodes of the network. For the present analysis, we label any 3-node that bounds a regular binary junction as a ‘normal’ or N-node. Likewise, any 4-node formed by two dislocations crossing each other is also labelled an N-node. All other nodes are regarded as multi-nodes or M-nodes, including the symmetric 4-nodes shown in Fig. 1e as well as 3-nodes produced by dissociation of the symmetric 4-nodes. Three types of links can now be defined with respect to the types of nodes they connect: NN links, NM links and MM links. Here, notation NM is used for any link that connects an N-node to an M-node. To compute the fraction of lines involved in the multi-junctions shown in Fig. 3d, we summed the lengths of all MM links with half the lengths of all NM links and divided this sum by the total length of all links. The colour scheme used in Fig. 4 is as follows: all MM links are shown in white, whereas the colour of NM links is graded from white at the M-nodes to green at the N-nodes. All NN links, including binary junctions, are shown in green. Atomistic simulation. The simulation volume was a small cube-shaped block of a perfect b.c.c. single crystal, 17 nm on each side. The initial geometry contained three dislocations with Burgers vectors 1/2[111], 1/2[111], and 1/2[111] intersecting at the block centre. The atom positions inside the block were then relaxed to mechanical equilibrium using the conjugate gradient method and an interatomic interaction function for molybdenum20. The atoms on the block surfaces were fixed throughout the simulation. To visualize crystal defects, only the atoms inside the block with energies exceeding the ideal bulk value by 0.095 eV are shown. Experiment. The experiments involved three steps: (1) compression of a singlecrystal molybdenum specimen to 1% total strain along the [001] axis, (2) cutting and thinning the deformed specimen along the (101) plane to obtain electron transparent foils, and (3) TEM observations using a set of reflection vectors g that can reveal multi-junctions. In the view shown in Fig. 2b the zone axis <[101] and the diffraction vector g ¼ [020], making all four dislocations entering the 4-node visible. The views in Fig. 2c and d were obtained using g ¼ [121] and g ¼ [121] which made lines b1 ¼ 1/2[111] and b3 ¼ 1/2[111] invisible owing to the gzb ¼ 0 condition. To access additional diffraction vectors, the specimen was tilted to a new zone axis (<[201]) making it possible to identify the Burgers vectors of two remaining dislocations in a similar manner: b2 ¼ 1/2[111] and b4 ¼ 1/2[111], respectively. Received 14 November 2005; accepted 14 February 2006. 1. Orowan, E. Zur kristallplastizitat. Z. Phys. 89, 605–-659 (1934). 2. Taylor, G. The mechanism of plastic deformation in crystal. Part I. Theoretical. Proc. R. Soc. A 145, 363–-404 (1934). 3. Polanyi, M. Uber eine Art Gutterstorung, die einen kristall plastich machen konnte. Z. Phys. 89, 660–-664 (1934). 4. Hirsch, P. B., Horne, R. W. & Whelan, M. J. Direct observations of the arrangement and motion of dislocations in aluminium. Phil. Mag. 1, 677–-684 (1956). 5. Hirth, J. P. & Lothe, J. Theory of Dislocations 2nd edn, 9–-13 (Wiley, New York, 1982). 6. Saada, G. Sur le durcissement duˆ a` recombinaison des dislocations. Acta Metall. 8, 841–-852 (1960). 7. Franciosi, P., Berveiller, M. & Zaoui, A. Latent hardening in copper and aluminium single crystals. Acta Metall. 28, 273–-283 (1980). 8. Guiu, F. & Pratt, P. L. The effect of orientation on the yielding and flow of molybdenum single crystals. Phys. Status Solidi 15, 539–-552 (1966). 9. Hirth, J. P. & Lothe, J. Theory of Dislocations 2nd edn, 269–-270 (Wiley, New York, 1982). 10. Basinski, S. J. & Basinski, Z. S. in Dislocations in Solids Vol. 4 (ed. Nabarro, F. R. N.) 261–-362, (North Holland, Amsterdam, 1979). 11. Nabarro, F. R. N., Basinski, Z. S. & Holt, D. B. The plasticity of pure single crystals. Adv. Phys. 50, 193–-233 (1964). 12. Madec, S. R., Devincre, B. & Kubin, L. P. From dislocation junctions to forest hardening. Phys. Rev. Lett. 89, 255508 (2002). 13. Ha¨hner, P. & Zaiser, M. Dislocation dynamics and work hardening of fractal dislocation cell structures. Mater. Sci. Eng. A 272, 443–-454 (1999). 14. Abraham, F. F. et al. Simulating materials failure by using up to one billion atoms and the world’s fastest computer: Work hardening. Proc. Natl. Acad. Sci. USA 99, 5783–-5787 (2002). 15. Frank, F. C. & Read, W. T. Multiplication processes for slow moving dislocations. Phys. Rev. 79, 722–-723 (1950). 16. Bulatov, V. V. et al. Scalable line dynamics in ParaDiS. Supercomputing khttp:// www.sc-conference.org/sc2004/schedule/pdfs/pap206.pdfl (2004). 17. Mughrabi, H. Dislocation wall and cell structures and long-range internal stresses in deformed metal crystals. Acta Metall. 31, 1367–-1379 (1983). 18. Kubin, L. P. & Canova, G. R. The modelling of dislocation patterns. Scr. Metall. Mater. 27, 957–-962 (1992). 19. Ghoniem, N. M., Huang, J. M. & Wang, Z. Q. Affine covariant-contravariant vector forms for the elastic field of parametric dislocations in isotropic crystals. Phil. Mag. Lett. 82, 55–-63 (2002). 20. Finnis, M. W. & Sinclair, J. S. A simple empirical N-body potential for transition metals. Phil. Mag. A 50, 45–-55 (1984). Supplementary Information is linked to the online version of the paper at www.nature.com/nature. Acknowledgements This work was supported by the US DOE Office of Basic Energy Sciences and the NNSA ASC program. We thank E. Chandler and C. Mailhiot for their encouragement and unwavering support of the ParaDiS code development effort, D. Lassila for providing the purified molybdenum crystals, M. LeBlanc for performing the mechanical tests, R. Cook for graphic design and F. Abraham and G. Campbell for critical reading and editorial suggestions for the manuscript. This work was performed under the auspices of the US DOE by the Lawrence Livermore National Laboratory. Author Information Reprints and permissions information is available at npg.nature.com/reprintsandpermissions. The authors declare no competing financial interests. Correspondence and requests for materials should be addressed to V.V.B. (bulatov1@llnl.gov). LETTERS NATURE|Vol 440|27 April 2006 1178