> For the complete documentation index, see [llms.txt](https://3d-materials-lab.gitbook.io/3dmaterialslab-tutorials/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://3d-materials-lab.gitbook.io/3dmaterialslab-tutorials/dynamics-in-crystalline-solids/first-principles-calculations-of-atomic-diffusion-in-crystalline-solids.md).

# First-principles Calculations of Atomic Diffusion in Crystalline Solids

*Written by Cheng-Wei Lee* (clee2 \[at] mines \[dot] edu)

## Atomic diffusion mechanisms&#x20;

There are multiple textbook diffusion mechanisms and the vacancy-mediated diffusion is one of the most common diffusion mechanisms in ionic compounds.  The associated diffusion coefficient can be written as:

$$
D=\frac{1}{6}fZ\_{f}Z\_{m}l^{2}C\_{D}\Gamma
$$

where $$\frac{1}{6}$$ is for 3D system ( $$\frac{1}{2n}$$ and $$n=3$$) and $$f$$, $$Z\_{f}$$, and $$Z\_{m}$$ are geometry factors, which depend on the crystal structure. $$l$$ is the distance the vacancy travel for each jump, $$C\_{D}$$ is the concentration of participating defects, e.g. vacancies. Lastly, $$\Gamma$$ is the successful  jump frequency,

$$
\Gamma = \nu^{\*}e^{-\frac{\Delta E\_{m}}{k\mathrm{\_{B}}T}}
$$

where $$\nu^{\*}$$ is the attempted jump frequency, $$\Delta E\_{m}$$ is the migration barrier of a diffusing atom (see figure below)&#x20;

<figure><img src="/files/JlCUrSurP7yHukkyVmS0" alt=""><figcaption><p>A schematic for vacancy-mediated diffusion in crystalline solids. Dashed open circles indicate vacancy sites. </p></figcaption></figure>

Therefore, the diffusion coefficient for vacancy-mediated diffusion in crystalline solids can be calculated from first principles:

1. The attempt frequency can be approximated by the local dynamic matrix of the moving atom and the force due to smaller change in moving atom can be calculated using DFT.
2. The defect concentration can be calculated using the supercell approach (see defect calculations)
3. The migration barrier can be calculated using [nudged elastic band method ](broken://pages/uPpHi40ITzK4Lvpzi0fI)


---

# Agent Instructions
This documentation is published with GitBook. GitBook is the documentation platform designed so that both humans and AI agents can read, navigate, and reason over technical content effectively. Learn more at gitbook.com.

## Querying This Documentation
If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question.

Perform an HTTP GET request on the current page URL with the `ask` query parameter, and the optional `goal` query parameter:

```
GET https://3d-materials-lab.gitbook.io/3dmaterialslab-tutorials/dynamics-in-crystalline-solids/first-principles-calculations-of-atomic-diffusion-in-crystalline-solids.md?ask=<question>&goal=<endgoal>
```

`ask` is the immediate question: it should be specific, self-contained, and written in natural language.
`goal` is optional and describes the broader end goal you are ultimately trying to accomplish on behalf of the user. GitBook uses it to tailor the answer towards what is most useful for that goal.

The response will contain a direct answer to the question and relevant excerpts and sources from the documentation.

Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.
