281 research outputs found
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 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 screw dislocations in binary
tungsten-transition metal alloys (WTM) were investigated using
first principles electronic structure calculations. The periodic quadrupole
approach was applied to model the structure of 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 elastic constant and
increase of elastic anisotropy A=. Alloying W with a group
VIIIA transition metal changes the structure of the dislocation core from
symmetric to asymmetric, similar to results obtained for WRe
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
- …