Sweep-tracing algorithm: in silico slip crystallography and tension-compression asymmetry in BCC metals

Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

Survey of Grain Boundary Energies in Tungsten and Beta-Titanium at High Temperature

Heat treatment is a necessary means to obtain desired properties for most of the materials. Thus, the grain boundary (GB) phenomena observed in experiments actually reflect the GB behaviors at relatively high temperature to some extent. In this work, 405 different GBs were systematically constructed for body-centered cubic (BCC) metals and the grain boundary energies (GBEs) of these GBs were calculated with molecular dynamics for W at 2400 K and β-Ti at 1300 K and by means of molecular statics for Mo and W at 0 K. It was found that high temperature may result in the GB complexion transitions for some GBs, such as the Σ11{332}{332} of W. Moreover, the relationships between GBEs and sin(θ) can be described by the functions of the same type for different GB sets having the same misorientation axis, where θ is the angle between the misorientation axis and the GB plane. Generally, the GBs tend to have lower GBE when sin(θ) is equal to 0. However, the GB sets with the <110> misorientation axis have the lowest GBE when sin(θ) is close to 1. Another discovery is that the local hexagonal-close packed α phase is more likely to form at the GBs with the lattice misorientations of 38.9°/<110>, 50.5°/<110>, 59.0°/<110> and 60.0°/<111> for β-Ti at 1300 K.

Revealing and Controlling the Core of Screw Dislocations in BCC Metals

Body-centred-cubic (BCC) transition metals (TMs) tend to be brittle at low temperatures, posing significant challenges in their processing and major concerns for damage tolerance in critical load-carrying applications. The brittleness is largely dictated by the screw dislocation core structure; the nature and control of which has remained a puzzle for nearly a century. Here, we introduce a universal model and a physics-based material index χ that guides the manipulation of dislocation core structure in all pure BCC metals and alloys. We show that the core structure, commonly classified as degenerate (D) or non-degenerate (ND), is governed by the energy difference between BCC and face-centred cubic (FCC) structures and χ robustly captures this key quantity. For BCC TMs alloys, the core structure transition from ND to D occurs when χ drops below a threshold, as seen in atomistic simulations based on nearly all extant interatomic potentials and density functional theory (DFT) calculations of W-Re/Ta alloys. In binary W-TMs alloys, DFT calculations show that χ is related to the valence electron concentration at low to moderate solute concentrations, and can be controlled via alloying. χ can be quantitatively and efficiently predicted via rapid, low-cost DFT calculations for any BCC metal alloys, providing a robust, easily applied tool for the design of ductile and tough BCC alloys.