EPSRC Reference: 
EP/R043612/1 
Title: 
Boundary Conditions for Atomistic Simulation of Material Defects 
Principal Investigator: 
Ortner, Professor C 
Other Investigators: 

Researcher CoInvestigators: 

Project Partners: 

Department: 
Mathematics 
Organisation: 
University of Warwick 
Scheme: 
Standard Research 
Starts: 
01 October 2018 
Ends: 
31 March 2022 
Value (£): 
442,958

EPSRC Research Topic Classifications: 
Continuum Mechanics 
Eng. Dynamics & Tribology 
Mathematical Analysis 
Numerical Analysis 

EPSRC Industrial Sector Classifications: 
No relevance to Underpinning Sectors 


Related Grants: 

Panel History: 

Summary on Grant Application Form 
Atomistic simulations are an indispensable tool of modern materials science, solid state physics and chemistry, as they allow scientists to study individual atoms and molecules in a way that is impossible in laboratory experiments. Understanding atomistic processes opens up avenues for the manipulation of matter at the atomic scale in order to achieve superior material properties for applications in science and engineering.
One of the most common tasks of atomistic materials modelling is to determine properties of crystalline defects, including their atomic structure, formation, activation and ionisation energies, from which electronic and atomistic mechanisms of chemical reactivity, charge mobility, etc., can be directly discovered, and mesoscopic material properties or coarsegrained models (e.g., employed in kinetic MonteCarlo, discrete dislocation dynamics, continuum fracture laws, transport simulations) can be derived.
Defects distort the surrounding host lattice, generating longranging elastic (and possibly also electrostatic) fields. Since practical schemes necessarily work in small computational domains they cannot explicitly resolve these farfields but must employ artificial boundary conditions (e.g., periodic boundary conditions) to emulate the elastic bulk. This approximation gives rise to a simulation error that must be controlled and ideally balanced against other model and/or discretisation errors. For example, for a wide class of defects encompassing all (neutral) point defects and straight dislocations it is shown by Ehrlacher, Ortner and Shapeev (2016) that the geometry error decays with a universal rate O(N^{1/2}) where N denotes the number of atoms in the computational cell. For a cubic scaling computational chemistry model, this slow rate is particularly severe. For cracks, it turns out that the standard models even yield schemes that are divergent in N.
This extremely slow rate of convergence or even divergence represents both a theoretical and computational challenge, which we propose to address in this project. Specifically, we will develop a hierarchy of highaccuracy boundary conditions for four common classes of defects: charge neutral point defects, dislocations, cracks, and charged defects. At its core, this research involves the development of a range of new analytical tools to describe elastic and polarisation fields in crystalline solids and how they are coupled to defect cores. The analytical results will feed directly back into materials simulation methodology through new algorithms and open source software. The effect of these new algorithms will be to enhance both the reliability and efficiency of atomistic simulation of materials, and enable simulation of particularly complex defect structures that have so far been inaccessible with conventional tools.

Key Findings 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk

Potential use in nonacademic contexts 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk

Impacts 
Description 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk 
Summary 

Date Materialised 


Sectors submitted by the Researcher 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk

Project URL: 

Further Information: 

Organisation Website: 
http://www.warwick.ac.uk 