Void-induced cross slip of screw dislocations in fcc copper

Pinning interaction between a screw dislocation and a void in fcc copper is investigated by means of molecular dynamics simulation. A screw dislocation bows out to undergo depinning on the original glide plane at low temperatures, where the behavior of the depinning stress is consistent with that obtained by a continuum model. If the temperature is higher than 300 K, the motion of a screw dislocation is no longer restricted to a single glide plane due to cross slip on the void surface. Several depinning mechanisms that involve multiple glide planes are found. In particular, a depinning mechanism that produces an intrinsic prismatic loop is found. We show that these complex depinning mechanisms significantly increase the depinning stress

Effect of nickel on point defects diffusion in Fe – Ni alloys

Iron-Nickel alloys are perspective alloys as nuclear energy structural materials because of their good radiation damage tolerance and mechanical properties. Understanding of experimentally observed features such as the effect of Ni content to radiation defects evolution is essential for developing predictive models of radiation. Recently an atomic-scale modelling study has revealed one particular mechanism of Ni effect related to the reduced mobility of clusters of interstitial atoms in Fe-Ni alloys. In this paper we present results of the microsecond-scale molecular dynamics study of point defects, i.e. vacancies and self-interstitial atoms, diffusion in Fe-Ni alloys. It is found that the addition of Ni atoms affects diffusion processes: diffusion of vacancies is enhanced in the presence of Ni, whereas diffusion of interstitials is reduced and these effects increase at high Ni concentration and low temperature. The role of Ni solutes in radiation damage evolution in Fe-Ni ferritic alloys is discussedPeer ReviewedPostprint (author's final draft

Point defect dynamics in bcc metals

We present an analysis of the time evolution of self-interstitial atom and vacancy (point defect) populations in pure bcc metals under constant irradiation flux conditions. Mean-field rate equations are developed in parallel to a kinetic Monte Carlo (kMC) model. When only considering the elementary processes of defect production, defect migration, recombination and absorption at sinks, the kMC model and rate equations are shown to be equivalent and the time evolution of the point defect populations is analyzed using simple scaling arguments. We show that the typically large mismatch of the rates of interstitial and vacancy migration in bcc metals can lead to a vacancy population that grows as the square root of time. The vacancy cluster size distribution under both irreversible and reversible attachment can be described by a simple exponential function. We also consider the effect of highly mobile interstitial clusters and apply the model with parameters appropriate for vanadium and $\alpha-$iron.Comment: to appear in Phys. Rev.

Solute effects on edge dislocation pinning in complex alpha-Fe alloys

Reactor pressure vessel steels are well-known to harden and embrittle under neutron irradiation, mainly because of the formation of obstacles to the motion of dislocations, in particular, precipitates and clusters composed of Cu, Ni, Mn, Si and P. In this paper, we employ two complementary atomistic modelling techniques to study the heterogeneous precipitation and segregation of these elements and their effects on the edge dislocations in BCC iron. We use a special and highly computationally efficient Monte Carlo algorithm in a constrained semi-grand canonical ensemble to compute the equilibrium configurations for solute clusters around the dislocation core. Next, we use standard molecular dynamics to predict and analyze the effect of this segregation on the dislocation mobility. Consistently with expectations our results confirm that the required stress for dislocation unpinning from the precipitates formed on top of it is quite large. The identification of the precipitate resistance allows a quantitative treatment of atomistic results, enabling scale transition towards larger scale simulations, such as dislocation dynamics or phase field.Fil: Pascuet, Maria Ines Magdalena. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Comision Nacional de Energia Atomica. Centro Atomico Constituyentes. Departamento de Materiales; ArgentinaFil: Martínez, E.. Los Alamos National High Magnetic Field Laboratory; Estados UnidosFil: Monnet, G.. EDF–R&D; FranciaFil: Malerba, L.. SCK•CEN. Structural Materials Expert Group. Nuclear Materials Institute; Bélgic

Absorbing Boundary Conditions for Molecular Dynamics and Multiscale Modeling

We present an application of differential equation based local absorbing boundary conditions to molecular dynamics. The absorbing boundary conditions result in the absorbtion of the majority of waves incident perpendicular to the bounding surface. We demonstrate that boundary conditions developed for the wave equation can be applied to molecular dynamics. Comparisons with damping material boundary conditions are discussed. The concept is extended to the formulation of an atomistic-continuum multiscale scheme with handshaking between the regions based on absorbing boundary conditions. The multiscale model is effective in minimizing spurious reflections at the interface

The MOLDY short-range molecular dynamics package

We describe a parallelised version of the MOLDY molecular dynamics program. This Fortran code is aimed at systems which may be described by short-range potentials and specifically those which may be addressed with the embedded atom method. This includes a wide range of transition metals and alloys. MOLDY provides a range of options in terms of the molecular dynamics ensemble used and the boundary conditions which may be applied. A number of standard potentials are provided, and the modular structure of the code allows new potentials to be added easily. The code is parallelised using OpenMP and can therefore be run on shared memory systems, including modern multicore processors. Particular attention is paid to the updates required in the main force loop, where synchronisation is often required in OpenMP implementations of molecular dynamics. We examine the performance of the parallel code in detail and give some examples of applications to realistic problems, including the dynamic compression of copper and carbon migration in an iron-carbon alloy

The influence of transition metal solutes on dislocation core structure and values of Peierls stress and barrier in tungsten

Several transition metals were examined to evaluate their potential for improving the ductility of tungsten. The dislocation core structure and Peierls stress and barrier of $1/2$ screw dislocations in binary tungsten-transition metal alloys (W$_{1-x}$TM$_{x}$) were investigated using first principles electronic structure calculations. The periodic quadrupole approach was applied to model the structure of $1/2$ dislocation. Alloying with transition metals was modeled using the virtual crystal approximation and the applicability of this approach was assessed by calculating the equilibrium lattice parameter and elastic constants of the tungsten alloys. Reasonable agreement was obtained with experimental data and with results obtained from the conventional supercell approach. Increasing the concentration of a transition metal from the VIIIA group, i.e. the elements in columns headed by Fe, Co and Ni, leads to reduction of the $C^\prime$ elastic constant and increase of elastic anisotropy A=$C_{44}/C^\prime$. Alloying W with a group VIIIA transition metal changes the structure of the dislocation core from symmetric to asymmetric, similar to results obtained for W$_{1-x}$Re$_{x}$ alloys in the earlier work of Romaner {\it et al} (Phys. Rev. Lett. 104, 195503 (2010))\comments{\cite{WRECORE}}. In addition to a change in the core symmetry, the values of the Peierls stress and barrier are reduced. The latter effect could lead to increased ductility in a tungsten-based alloy\comments{\cite{WRECORE}}. Our results demonstrate that alloying with any of the transition metals from the VIIIA group should have similar effect as alloying with Re.Comment: 12 pages, 8 figures, 3 table

Direct Observation of Martensitic Phase-Transformation Dynamics in Iron by 4D Single-Pulse Electron Microscopy

The in situ martensitic phase transformation of iron, a complex solid-state transition involving collective atomic displacement and interface movement, is studied in real time by means of four-dimensional (4D) electron microscopy. The iron nanofilm specimen is heated at a maximum rate of ∼10^(11) K/s by a single heating pulse, and the evolution of the phase transformation from body-centered cubic to face-centered cubic crystal structure is followed by means of single-pulse, selected-area diffraction and real-space imaging. Two distinct components are revealed in the evolution of the crystal structure. The first, on the nanosecond time scale, is a direct martensitic transformation, which proceeds in regions heated into the temperature range of stability of the fcc phase, 1185−1667 K. The second, on the microsecond time scale, represents an indirect process for the hottest central zone of laser heating, where the temperature is initially above 1667 K and cooling is the rate-determining step. The mechanism of the direct transformation involves two steps, that of (barrier-crossing) nucleation on the reported nanosecond time scale, followed by a rapid grain growth typically in ∼100 ps for 10 nm crystallites

Effect of collision cascades on dislocations in tungsten: A molecular dynamics study

Tungsten (W) is the prime candidate material for the divertor and other plasma-facing components in DEMO. The point defects (i.e. vacancies and self-interstitials) produced in collision cascades caused by incident neutrons aggregate into dislocation loops (and voids), which strongly affect the mechanical properties. The point defects also interact with existing microstructural features, and understanding these processes is crucial for modelling the long term microstructural evolution of the material under fusion conditions. In this work, we performed molecular dynamics simulations of cascades interacting with initially straight edge dislocation dipoles. It was found that the residual vacancy number usually exceeds the residual interstitial number for cascades interacting with vacancy type dipoles, but for interstitial type dipoles these are close. We observed that a cascade near a dislocation promotes climb, i.e. it facilitates the movement of point defects along the climb direction. We also observed that the dislocations move easily along the glide direction, and that kinks are formed near the centre of the cascade, which then facilitate the movement of the dipoles. Some dipoles are sheared off by the cascade, and this is dependent on PKA energy, position, direction, and the width of dipole