Lattice Green function for extended defect calculations: Computation and error estimation with long-range forces

Abstract

Computing the atomic geometry of lattice defects--point defects, dislocations, crack tips, surfaces, or boundaries--requires an accurate coupling of the local strain field to the long-range elastic field. Periodic boundary conditions used by classical potentials or density-functional theory may not accurately reproduce the correct bulk response to an isolated defect; this is especially true for dislocations. Recently, flexible boundary conditions have been developed to produce the correct long-range strain field from a defect--effectively "embedding" a defect in a finite cell with infinite bulk response, isolating it from either periodic images or free surfaces. Flexible boundary conditions require the calculation of the bulk response with the lattice Green function (LGF). While the LGF can be computed from the dynamical matrix, for supercell methods (periodic boundary conditions) it can only be calculated up to a maximum range. We illustrate how to accurately calculate the lattice Green function and estimate the error using a cutoff dynamical matrix combined with knowledge of the long-range behavior of the lattice Green function. The effective range of deviation of the lattice Green function from the long-range elastic behavior provides an important length scale in multiscale quasi-continuum and flexible boundary-condition calculations, and measures the error introduced with periodic boundary conditions.

35 pages, 6 figures

References in corpus (2)

• Ab Initio Study of Screw Dislocations in Mo and Ta: A new picture of plasticity in bcc transition metals
• Coarse-grained molecular dynamics: Nonlinear finite elements and finite temperature

Cited by in corpus (24)

• First-principles data for solid-solution strengthening of magnesium: From geometry and chemistry to properties
• Cu/Ag EAM Potential Optimized for Heteroepitaxial Diffusion from ab initio Data
• Deep Learning Enabled Strain Mapping of Single-Atom Defects in 2D Transition Metal Dichalcogenides with Sub-picometer Precision
• Prediction of thermal cross-slip stress in magnesium alloys from direct first principles data
• Automatic numerical evaluation of vacancy-mediated transport for arbitrary crystals: Onsager coefficients in the dilute limit using a Green function approach
• Universal correlation between electronic factors and solute-defect interactions in bcc refractory metals
• Electronic-structure study of an edge dislocation in Aluminum and the role of macroscopic deformations on its energetics
• Computing energy barriers for rare events from hybrid quantum/classical simulations through the virtual work principle
• Analysis of Boundary Conditions for Crystal Defect Atomistic Simulations
• Geometries of edge and mixed dislocations in bcc Fe from first principles calculations
• Quantized Dislocations
• Direct calculation of lattice Green function with arbitrary interactions for general crystals
• First-principles calculation of H vibrational excitations at a dislocation core of Pd
• Predictive Crystal Plasticity Modeling of Single Crystal Nickel Based on First-Principles Calculations
• Convergence rate for numerical computation of the lattice Green's function
• Asymptotic Expansion of the Elastic Far-Field of a Crystalline Defect
• Lattice Green's function for crystals containing a planar interface
• PEXSI-$\Sigma$: A Green's function embedding method for Kohn-Sham density functional theory
• Core Energies in Isolated Edge and Mixed Dislocations in BCC Fe from First-principles Energy Density Method
• Line defects in infinite networks of resistors
• On the Cauchy-Born Approximation at Finite Temperature
• The Green tensor of Mindlin's anisotropic first strain gradient elasticity
• Sound attenuation derived from quenched disorder in solids
• Surface lattice Green's functions for high-entropy alloys