Data for these images was generated via atomistic computer simulation using embedded atom potentials. For this project, computations are carried out simulating a variety of dislocations in the B2 phase of NiAl although it is possible to apply these visualization techniques to the results of most any intermetallic alloy simulation. Calculations were done for both stoichiometric and non-stoichiometric simulation blocks with varying degrees of applied shear.
A stoichiometric pure edge dislocation l=[0-10], [100](001) with varying shear was simulated. The Peierls stress was found to be 2600MPa for the lower energy Ni-rich core and 3700MPa for the Al-rich core. The relaxed lattice along with the resulting Y2 invariant field of the lower energy configuration may be seen here. Ni atoms are blue while Al atoms are pink. The same dislocation may be viewed as a movie where increasing applied shear has been cast over time such that changes in the core's structure and Y2 strain field may be clearly seen before during and after the Peierls stress is reached. The shear has been applied such that glide to the left will occur. The animation shows a great deal of spreading in the strain field in the direction in which slip will occur. Note that the strain that has built up is periodically released as movement occurs and although the strain field is changing as stress is applied, the lattice itself does not appear to change unless movement is occurring.
The same pure edge dislocation was simulated in Ni-48Al by randomly substituting Al atoms with Ni atoms for different random seeds. Peierls stresses of 4000MPa and 3000MPa, respectively, were found for the Ni- and Al-rich cores. The relaxed and unstressed Ni-rich cores resulting from two different random seeds may be seen for seed1, seed2 and seed4.
The Ni-rich core was again simulated, this time with singly placed antisites and vacancies. Images of Ni antisites, Al vacancies and Ni vacancies (which are thought to be energetically unfavored) were generated: Ni antisite 1 , Al vacancy 1 and Ni vacancy 2 . Note that the Ni vacancy causes no visible change in the core structure and very little change in the strain field even though it was placed almost at the center of the core.
Multiple Al vacancies were placed at key positions with respect
to the core in the simulations resulting in the following images:
vacancy configuration 3 and
vacancy configuration 4.
Other dislocation geometries were simulated on the same slip plane having a mixed screw and edge character. Such dislocations have a Y2 as well as a Y1 component to the strain invariant. Images and Peierls stresses are:
A pure screw dislocation, l=[001], [001](010), was simulated. Although the shear was applied such that the dislocation would glide in the (010) plane, the dislocation moved along the (011) plane at an applied stress of 100MPa.
Numerous simulations were carried out on the (01-1) slip plane. The l=[011], [100](011) pure edge dislocation was simulated for various stoichiometries and the Peierls stresses were found. The Ni-48Al simulations were found to be not sensitive to placement of the antisites whereas the Ni-42Al simulations were found to be very sensitive to the given random seed:
The animation of this pure edge dislocation, indicates that the outer limbs of the invariant field slowly contract before the core moves. The onset of movement is quite sudden and occurs at an applied shear of 410MPa. Other than the contraction of the strain field, there is no warning that glide is about to occur. As the core moves, observe the deformation in both the lattice and in the strain field which is represented as color.
This dislocation has also been simulated with individual Al vacancies introduced in the core region. Images are placed here for vacancy position 1, vacancy position 3, vacancy position 5 and vacancy position 8.
The pure screw, l=[100], [100](01-1), was found to have a Peierls stress of 390MPa (shown) or 70MPa depending upon its initial configuration as discussed below. The animation of this pure screw dislocation shows the Y1 strain field undergoing a very slow contraction before sudden changes which occur at both 316MPa and 390MPa. In two steps, the core appears to reorient its longer axis along the direction of applied stress and then glide easily. Note that, except for the orientation, the core appears unchanged before and after this transformation.
In fact, the two orientations of the screw are crystallographically equivalent. We have found that this core glides easily along its longer axis at 110MPa if the shear is initially applied along that axis. Even though the two core configurations are energetically equivalent, an energy barrier must be overcome in order to change between the two. This barrier accounts for the difference in Peierls stresses.
More recently, we have been simulating B2 NiAl with titanium impurities. It is thought from just energy considerations that the Ti impurities will prefer to substitute themselves into the Al sublattice. Here are such simulations of Ni-47Al-3Ti and Ni-44Al-6Ti for the l=[0-10], [100](001) dislocation discussed above. The Ti atoms are represented in green. For comparison, simulations where the Ti atoms were substituted into the Ni sublattice were carried out. Here are images of 47Ni-50Al-3Ti and 44Ni-50Al-6Ti for the same dislocation where much higher strains have manifest.
Figure: The B2 crystal structure.
Click here to go to the 1994 Scientific Visualization and Multimedia class page.