Simulations of copper single crystals subjected to rapid shear

Abstract

We report on nonequilibrium molecular dynamics simulations of single crystals of copper experiencing rapid shear strain. A model system, with periodic boundary conditions, which includes a single dislocation dipole is subjected to a total shear strain of close to 10% on time-scales ranging from the instantaneous to 50 ps. When the system is strained on a time-scale short compared with a phonon period, the initial total applied shear is purely elastic, and the eventual temperature rise in the system due to the subsequent plastic work can be determined from the initial elastic strain energy. The rate at which this plastic work occurs, and heat is generated, depends on the dislocation velocity, which itself is a function of shear stress. A determination of the stress-dependence of the dislocation velocity allows us to construct a simple analytic model for the temperature rise in the system as a function of strain rate, and this model is found to be in good agreement with the simulations. For the effective dislocation density within the simulations, 7:8 Â 10 11 cm À2 , we find that applying the total shear strain on time-scales of a few tens of picoseconds greatly reduces the final temperature. We discuss these results in the context of the growing interest in producing high pressure, solid-state matter, by quasi-isentropic (rather than shock) compression.

The mechanics and physics of high-speed dislocations: a critical review

“…As explained in the Introduction, this mobility law is one of the main quantities through which this atomisticcontinuum divide is bridged. In the absence of experimental measurements, an obvious strategy to obtain such a relation is to perform MD simulations of gliding dislocations from which a quantitative velocity-stress relation to be used in DD simulations may be obtained [158,235,345,346,351,352,[355][356][357][358][359]. Associated instabilities affecting high speed dislocations concerning dislocation interactions at high speed have also been explored [360,361].…”

Bouncing window for colliding nanoparticles: Role of dislocation generation

“…where C is a proportionality constant [30][31][32]. For a qualitative discussion, we may assume that v d is constant for the same impact velocity, and the dissipated energy is hence proportional to the area under the dislocation length curves in Fig.…”

Dislocation kinematics: a molecular dynamics study in Cu

“…Such studies have identified novel phenomena such as the existence of stable transonic dislocation motion under high stress in metals [13]. These studies vary in the choice of the material crystal lattice, temperatures, dislocation type, range of applied stress, loading technique (strain rate versus NPT ensemble) and boundary conditions (periodic boundary conditions [14][15][16][17][18] versus thin films [13,[19][20][21][22][23][24]). Despite this considerable amount of work there remain both physical and computational questions unresolved, e.g.…”