Molecular Dynamics Simulations

Generalized stacking fault energy and nature of slip

Using molecular dynamics simulations, the nature of slip in nanocrystalline metals cannot be described in terms of the absolute value of the stacking fault energy—a correct interpretation requires the generalized stacking fault energy curve, involving both stable and unstable stacking fault energies In conclusion, whether the dislocation activity in a simulation is dominated by extended partials or full dislocations travelling through the grains in the simulated time span must be understood in terms of the ratio γsf/γusf. Similarly, the observation of twin nucleation via an extended partial on an adjacent plane is dependent on the ratio γutf/γusf. To understand the slip activity observed in simulations, all aspects of the GPF curve have to be incorporated, and not just the absolute value of γsf.

Nature Materials 3 (2004) 400
Physical Review B 76 (2007) 054124


Error: (1) can't find MDsimulENstrain_rate.gif at LIVE_WWW.PEM

The effect of strain rate in molecular dynamics

To study deformation mechanisms using molecular dynamics, very high strain rates have to be applied. The effect of lowering the strain rate by two orders of magnitude on the deformation characteristics of nanocrystalline Al was investigated. For the highest strain rate, the onset of dislocation propagation is delayed, resulting in a stress overshoot. With decreasing strain rate the grain-averaged resolved shear stress reduces and cross-slip occurs more frequently. However, even at the lowest applied strain rate the grain boundary network can not accommodate all arriving dislocations, illustrating the challenge to determine the rate-limiting deformation mechanisms for experimental conditions.

Philosophical Magazine 89 (2009) 3465

Dislocation nucleation and propagation in nc environment

Detailed analysis of the nucleation and propagation of dislocations from grain boundaries in nanocrystalline Al using molecular dynamics reveals that: (1) the leading and trailing partial dislocations are nucleated at different regions within a grain boundary or even different grain boundaries; (2) propagation is hindered by a pinning–depinning mechanism influenced by the relative orientation between the Burgers vector and the ledge geometry of the grain boundary and (3) depinning is thermally activated and therefore not purely mechanically driven. Detailed analysis of the atomic processes allowed correlation of the stress signature of slip events with typical dislocation processes related to nucleation and propagation.

Acta Materialia 54 (2006) 1975
Acta Materialia 56 (2008) 4846

Dislocation cross-slip

Constant strain rate molecular dynamics simulations of nanocrystalline Al demonstrate that a significant amount of dislocations that have nucleated at the grain boundaries, exhibit cross-slip via the Fleischer mechanism as they propagate through the grain. The grain boundary structure is found to strongly influence when and where cross-slip occurs, allowing the dislocation to avoid local stress concentrations that otherwise can act as strong pinning sites for dislocation propagation.

Physical Review Letters 100 (2008) 235501

Dislocation transmission trough a general grain boundary

The atomistic simulation reveals that the regions of misfit act as strong pinning centers for the newly created dislocation. Furthermore, the simulations lead us to the assumption that the choice of the Burgers vector of the outgoing dislocation is not only influenced by the resolved shear stress of the outgoing slip system, but also by its influence on the local atomic structure of the grain boundary.

Applied Physics Letters 91 (2007) 111914

The influence of twins on the mechanical properties of nc-Al

It was shown that the influence of the presence of grown-in twin boundaries in nc Cu and Ni is different from the case of nc-Al and can be explained by the values of the extrema of the Generalized Planar Fault energy curves. For Cu and Ni, where the energy barrier for twin migration is not so different from that for nucleation of a single partial dislocation, twin boundary migration does not become dominant over the usual dislocation activity from grain boundaries. Earlier work has shown that for Al, a material where the energy barrier for twin boundary migration is close to a factor of three lower then for the emission of a partial dislocation on a non-adjacent plane, twin boundary migration will become dominant, especially when twin planes are along orientations with a high Schmid factors.

Acta Materialia 52 (2004) 2259
Advanced engineering materials 7 (2005) 16
Applied Physics Letters 85 (2004) 5863

Grain boundary migration in nanocrystalline environment

Coupled grain boundary motion was simulated in a three-dimensional nanocrystalline Al grain boundary network using molecular dynamics. It is shown that, in spite of the triple junction constraints, a symmetrical R75 tilt boundary can migrate during the microplastic regime with the same coupling factor as when simulated in a bicrystal configuration. After reaching the full plastic regime, dislocations start coming into play, changing the grain boundary structure and hindering further coupled motion.

The simulated Ni sample contained over 2 million atoms arranged in 200 grains of 5 nm average size. Grain boundary migration occurs in conjunction with dislocation emission, it is however not clear if dislocation emission is a necessary condition for grain boundary motion in this configuration.. In experimental work, it is often difficult to distinguish between purely stress driven grain boundary motion and thermal effects. The present simulations help resolve this controversy, since in our simulations the temperature is maintained at 300 K and the deformation process occurs at a very fast strain rate, with limited time for thermally activated motion. Our results indicate that stress driven grain boundary migration is one of the ways in which plasticity can be accommodated by grain boundaries in a 3D nanocrystalline GB network.

Scripta Materialia 55 (2006) 695
Scripta Materialia 65 (2011) 151

The influence of Oxygen in GB migration

One of the important differences between simulation and experiments in grain boundary (GB)-dominated metallic structures is the lack of impurities such as oxygen in computational samples. A modified variable-charge method is used to study coupled GB motion in an Al bicrystal with a [112] symmetrical tilt GB in the presence of substitutional O, and compared with the stick–slip process identified in pure configurations.  It is found that the critical shear stress for migration of the GB increases linearly with the number of O atoms. These observations are then rationalized in terms of the internal stress signature of the O atoms in the vicinity of the boundary.

Modelling and Simulation in Materials Science and Engineering 16 (2008) 025006
Acta Materialia 57 (2009) 1988