Periodic image effects in dislocation modelling

The use of periodic boundary conditions for modelling crystal dislocations is predicated on one's ability to handle the inevitable image effects. This communication deals with an often overlooked mathematical subtlety involved in dealing with the periodic dislocation arrays, that is conditional convergence of the lattice sums of image fields. By analysing the origin of conditional convergence and the numerical artefacts associated with it, we establish a mathematically consistent and numerically efficient procedure for regularization of the lattice sums and the corresponding image fields. The regularized solutions are free from the artefacts caused by conditional convergence and regain periodicity and translational invariance of the periodic supercells. Unlike the other existing methods, our approach is applicable to general anisotropic elasticity and arbitrary dislocation arrangements. The capabilities of this general methodology are demonstrated by application to a variety of situations encountered in atomistic and continuum modelling of crystal dislocations. The applications include introduction of dislocations in the periodic supercell for subsequent atomistic simulations, atomistic calculations of the core energies and the Peierls stress and continuum dislocation dynamics simulations in three dimensions.     

Dislocation modeling in bcc lithium: A comparison between continuum and atomistic predictions in the modified embedded atoms method

In this study, the modified embedded-atom method (MEAM) was applied to compare the predictions of dislocation core properties obtained by molecular statics with the continuum predictions obtained in the framework of the simplified 1D-Peierls–Nabarro model. To this end, a set of four fictive Li potentials in the MEAM framework was proposed with the condition that all four potentials reproduce the same elastic constants, the same transition energies between bcc and fcc crystal structures, and between bcc and hcp crystal structures, while the unstable stacking fault energy on the plane {110} in the direction <111> was varied around the value predicted by first-principles. Within these potentials, direct atomistic calculations were performed to evaluate dislocation core properties such as dislocation half width and Peierls stress and the results were compared with continuum predictions. We found that the trends predicted by the Peierls–Nabarro model, i.e. (i) a decrease of the dislocation half width with increasing unstable stacking fault energy, and (ii) an increase of the Peierls stress with increasing the magnitude of the unstable stacking fault energy, were recovered using atomic calculations in the MEAM framework. Moreover, the magnitude of the dislocation half width and the Peierls stress calculated in the MEAM framework are in good agreement with the Peierls–Nabarro predictions when the dislocation half width is determined using a generic strategy. Specifically, the dislocation half width is defined as the distance for which the disregistery is included between b/4 and 3b/4. It was, therefore, demonstrated herein that the set of fictive potentials could be parameterized in the MEAM framework to validate or to disprove the continuum theory using atomistic methods.     

Free energy change of a dislocation due to a Cottrell atmosphere

The free energy reduction of a dislocation due to a Cottrell atmosphere of solutes is computed using a continuum model. We show that the free energy change is composed of near-core and far-field components. The far-field component can be computed analytically using the linearized theory of solid solutions. Near the core the linearized theory is inaccurate, and the near-core component must be computed numerically. The influence of interactions between solutes in neighbouring lattice sites is also examined using the continuum model. We show that this model is able to reproduce atomistic calculations of the nickel–hydrogen system, predicting hydride formation on dislocations. The formation of these hydrides leads to dramatic reductions in the free energy. Finally, the influence of the free energy change on a dislocation’s line tension is examined by computing the equilibrium shape of a dislocation shear loop and the activation stress for a Frank–Read source using discrete dislocation dynamics.     

A generalized Peierls–Nabarro model for kinked dislocations

We present a generalized Peierls–Nabarro model for curved dislocations incorporating directly the Peierls energies for both straight dislocations and dislocation kinks. In our model, the anisotropic elastic energy is calculated efficiently using the discrete Fourier transform on the discrete lattice sites of the slip plane, and the discreteness in both the elastic energy and the misfit energy is included. We have used our model to calculate the kink migration and nucleation energies of the 30° dislocations in silicon. The results agree well with those obtained using atomistic potentials and first principles calculations, and the experimental results.     

Probing the structure and energetics of dislocation cores in SiGe alloys through Monte Carlo simulations

We present a methodology for the investigation of dislocation energetics in segregated alloys based on Monte Carlo simulations which equilibrate the topology and composition of the dislocation core and its surroundings. An environment-dependent partitioning of the system total energy into atomic contributions allows us to link the atomistic picture to continuum elasticity theory. The method is applied to extract core energies and radii of 60 degrees glide dislocations in segregated SiGe alloys which are inaccessible by other methods.     

Can vacancies lubricate dislocation motion in aluminum?

The interaction of vacancy with dislocations in Al is studied using the semidiscrete variational Peierls-Nabarro model with ab initio determined gamma surface. For the first time, we confirm theoretically the so-called vacancy lubrication effect on dislocation motion in Al, a discovery that can settle a long-standing controversy in dislocation theory for fcc metals. We provide insights into the lubrication effect by exploring the connection between dislocation mobility and its core width. We predict an increased dislocation splitting in the presence of vacancy. We find that on average there is a weak repulsion between vacancies and dislocations which is independent of dislocation character.     

二相介质中一种位错构型的弹性场及所受的力

本文建立了二相介质中与相界斜交位错的一种理想构型。根据位错各向异性弹性理论的普遍方法,以及关于界面的Gebbia-Eshelby定理,计算得到此位错在两相中的弹性场,以及位错线段所受“像力”的数值解。本文所提出的模型及算法具有一般适用性。结果可用于进一步计算此位错与其它缺陷的相互作用。对于考虑这种斜交位错在无扩散相变中的作用,本工作也是一个有意义的起点。 关键词：     

Atomistic Modeling of Point Defects and Diffusion in Copper Grain Boundaries

The atomic structure of several symmetrical tilt grain boundaries (GBs) in Cu and their interaction with vacancies and interstitials as well as self-diffusion are studied by molecular statics, molecular dynamics, kinetic Monte Carlo (KMC), and other atomistic simulation methods. Point defect formation energy in the GBs is on average lower than in the lattice but variations from site to site within the GB core are very significant. The formation energies of vacancies and interstitials are close to one another, which makes the defects equally important for GB diffusion. Vacancies show interesting effects such as delocalization and instability at certain GB sites. They move in GBs by simple vacancy-atom exchanges or by long jumps involving several atoms. Interstitial atoms can occupy relatively open positions between atoms, form split dumbbell configurations, or form highly delocalized displacement zones. They diffuse by direct jumps or by the indirect mechanism involving a collective displacement of several atoms. Diffusion coefficients in the GBs have been calculated by KMC simulations using defect jump rates determined within the transition state theory. GB diffusion can be dominated by vacancies or interstitials, depending on the GB structure. The diffusion anisotropy also depends on the GB structure, with diffusion along the tilt axis being either faster or slower than diffusion normal to the tilt axis. In agreement with Borisov's correlation, the activation energy of GB diffusion tends to decrease with the GB energy.     

Plastic deformation under high-rate loading: The multiscale approach

A two-level approach has been proposed for describing the plastic deformation under high-rate loading of metals. The characteristics of the motion of dislocations under shear stresses have been investigated at the atomistic level by using the molecular dynamics simulation. The macroscopic motion of a material has been described at the continuum level with the use of the model of continuum mechanics with dislocations, which uses information obtained at the atomistic level on the dislocation dynamics. The proposed approach has been used to study the evolution of the dislocation subsystem under shock-wave loading of an aluminum target. The behavior of the dynamic yield stress with an increase in the temperature has been analyzed. The results of the calculations are in good agreement with experimental data.     

Elastic energy of a straight dislocation and contribution from core tractions

We derive an expression of the core traction contribution to the dislocation elastic energy within linear anisotropic elasticity theory using the sextic formalism. With this contribution, the elastic energy is a state variable consistent with the work of the Peach–Koehler forces. This contribution needs also to be considered when extracting from atomic simulations core energies. The core energies thus obtained are real intrinsic dislocation properties: they do not depend on the presence and position of other defects. This is illustrated by calculating core energies of edge dislocation in bcc iron, where we show that dislocations gliding in {110} planes are more stable than those gliding in {112} planes.     

A non-singular continuum theory of point defects using gradient elasticity of bi-Helmholtz type

In this paper, we develop a non-singular continuum theory of point defects based on a second strain gradient elasticity theory, the so-called gradient elasticity of bi-Helmholtz type. Such a generalised continuum theory possesses a weak nonlocal character with two internal material lengths and provides a mechanics of defects without singularities. Gradient elasticity of bi-Helmholtz type gives a natural and physical regularisation of the classical singularities of defects, based on higher order partial differential equations. Point defects embedded in an isotropic solid are considered as eigenstrain problem in gradient elasticity of bi-Helmholtz type. Singularity-free fields of point defects are presented. The displacement field as well as the first, the second and the third gradients of the displacement are derived and it is shown that the classical singularities are regularised in this framework. This model delivers non-singular expressions for the displacement field, the first displacement gradient and the second displacement gradient. Moreover, the plastic distortion (eigendistortion) and the gradient of the plastic distortion of a dilatation centre are also non-singular and are given in terms of a form factor (shape function) of a point defect. Singularity-free expressions for the interaction energy and the interaction force between two dilatation centres and for the interaction energy and the interaction force of a dilatation centre in the stress field of an edge dislocation are given. The results are applied to calculate the finite self-energy of a dilatation centre.     

Vibrational entropy of dislocations in Al

The region nearest to a lattice defect must be described by an atomistic model, while a continuum model suffices further away from the defect. We study such a separation into two regions for an edge dislocation. In particular we focus on the excess defect energy and vibrational entropy, when the dislocation core is described by a cluster of about 500–100?atoms, embedded in a large discrete and relaxed, but static, lattice. The interaction between the atoms is given by a potential of the embedded-atom model type referring to Al. The dynamic matrix of the vibrations in the cluster is fully diagonalized. The excess entropy ΔS near the core has positive and negative contributions, depending on the sign of the local strain. Typically, ΔS/k _B ≈ 2 per atomic repeat length along the dislocation core in fcc Al. In the elastic continuum region far from the dislocation core the excess entropy shows the same logarithmic divergence as the elastic energy. Although the work refers to a specific material and defect type, the results are of a generic nature.     

Molecular dynamics simulation of the interaction between a mixed dislocation and a stacking fault tetrahedron

A high number-density of nanometer-sized stacking fault tetrahedra are commonly found during irradiation of low stacking fault energy metals. The stacking fault tetrahedra act as obstacles to dislocation motion leading to increased yield strength and decreased ductility. Thus, an improved understanding of the interaction between gliding dislocations and stacking fault tetrahedra are critical to reliably predict the mechanical properties of irradiated materials. Many studies have investigated the interaction of a screw or edge dislocation with a stacking fault tetrahedron (SFT). However, atomistic studies of a mixed dislocation interaction with an SFT are not available, even though mixed dislocations are the most common. In this paper, molecular dynamics simulation results of the interaction between a mixed dislocation and an SFT in face-centered cubic copper are presented. The interaction results in shearing, partial absorption, destabilization or simple bypass of the SFT, depending on the interaction geometry. However, the SFT was not completely annihilated, absorbed or collapsed during a single interaction with a mixed dislocation. These observations, combined with simulation results of edge or screw dislocations, suggest that defect-free channel formation in irradiated copper is not likely by a single dislocation sweeping or destruction process, but rather by a complex mix of multiple shearing, partial absorption and defect cluster transportation that ultimately reduces the size of stacking fault tetrahedra within a localized region.     

Dislocation density-based constitutive model for the mechanical behaviour of irradiated Cu

Performance degradation of structural steels in nuclear environments results from the formation of a high number density of nanometre-scale defects. The defects observed in copper-based alloys are composed of vacancy clusters in the form of stacking fault tetrahedra and/or prismatic dislocation loops that impede the motion of dislocations. The mechanical behaviour of irradiated copper alloys exhibits increased yield strength, decreased total strain to failure and decreased work hardening as compared to their unirradiated behaviour. Above certain critical defect concentrations (neutron doses), the mechanical behaviour exhibits distinct upper yield points. In this paper, we describe the formulation of an internal state variable model for the mechanical behaviour of such materials subject to these (irradiation) environments. This model has been developed within a multiscale materials-modelling framework, in which molecular dynamics simulations of dislocation–radiation defect interactions inform the final coarse-grained continuum model. The plasticity model includes mechanisms for dislocation density growth and multiplication and for irradiation defect density evolution with dislocation interaction. The general behaviour of the constitutive (homogeneous material point) model shows that as the defect density increases, the initial yield point increases and the initial strain hardening decreases. The final coarse-grained model is implemented into a finite element framework and used to simulate the behaviour of tensile specimens with varying levels of irradiation-induced material damage. The simulation results compare favourably with the experimentally observed mechanical behaviour of irradiated materials.     

Strain induced deep electronic states around threading dislocations in GaN

Combining through-focus high-resolution transmission electron microscopy and hierarchical multiscale simulations consisting of density-functional theory, analytical empirical potentials, and continuum elastic theory we demonstrate the existence of a new dislocation type in GaN. In contrast with all previously identified or suggested dislocation structures in GaN, all core atoms are fully coordinated; i.e., no broken bonds occur, implying that the dislocation should be electrically inactive. However, as we show, the giant local strain-field around the dislocation core, in combination with the small lattice constant of GaN, causes deep defect states and thus electrically active edge dislocations independent on the specific core structure.     

Computer studies of dislocation properties using two-dimensional model systems

The physical properties of dislocations are investigated by means of a two-dimensional Lennard-Jones solid. Using a self consistent linear theory for the elastic properties of two-dimensional systems, analytic expressions for the energy and entropy of a dislocation in the center of a circular finite model are derived. The energy and entropy formulas are compared with the results of molecular dynamics calculations on the model. The resulting atomic structure of the dislocation core is compared with the structure of bubble raft dislocations and with the results of earlier computer calculations.     

Atomistic sliding mechanisms of the Σ=5 symmetric tilt grain boundary in bcc iron

Atomistic computer simulations were performed to investigate the mechanisms of grain-boundary sliding in bcc Fe using molecular statics and molecular dynamics with embedded-atom method interatomic potentials. For this study we have chosen the Σ?=?5, (310)[001] symmetrical tilt boundary with tilt angle θ?=?36.9°. Sliding was determined to be governed by grain-boundary dislocation activity with Burgers vectors belonging to the displacement shift complete lattice. The sliding process was found to occur through the nucleation and glide of partial grain-boundary dislocations, with a secondary grain-boundary structure playing an important role in the sliding process. Interstitial impurities and vacancies were introduced into the grain boundary to study their role as nucleation sites for the grain-boundary dislocations. While vacancies and H interstitials act as preferred nucleation sites, C interstitials to not.     

Wave processes in elastic-viscoplastic media with dislocations

The mechanisms of plane harmonic wave propagation in homogeneous and interfaced elastic-viscoplastic media are considered using the field theory of defects with kinematic identities of a dislocation-containing elastic continuum and dynamic equations of the gauge theory of dislocations. The reflection and refraction coefficients were determined for displacement waves and defect field waves with the defect field characterized by the dislocation density tensor and flux density tensor. The dependence of the coefficients on the parameters of the interfaced media is analyzed.