Skip to content
Snippets Groups Projects

Ssr mathjax

Merged Ssr mathjax
ssr_mathjax into master
Merged Anton Akhmerov requested to merge ssr_mathjax into master
Compare
and
3 files
+ 96
− 108
Compare changes
  • Side-by-side
  • Inline
Files
3
docs/10_xray_solutions.md
+ 77
− 96
@@ -6,102 +6,91 @@ from math import pi
# Solutions for lecture 10 exercises
## Warm-up exercises
1.
??? hint "Small hint"
You can make use of the [scalar triple product](https://en.wikipedia.org/wiki/Triple_product#Scalar_triple_product).
2.
If $\mathbf{k}-\mathbf{k'}\neq \mathbf{G}$, then the argument of the exponent has a phase factor dependent on the real-space lattice points.
1. Hint: you can make use of the [scalar triple product](https://en.wikipedia.org/wiki/Triple_product#Scalar_triple_product).
2. If $\mathbf{k}-\mathbf{k'}\neq \mathbf{G}$, then the argument of the exponent has a phase factor dependent on the real-space lattice points.
Because we sum over each of these lattice points, each argument has a different phase.
Summing over all these phases results in an average amplitude of 0, resulting in no intensity peaks.
3.
No, there is a single atom, and thus only one term in the structure factor.
3. No, there is a single atom, and thus only one term in the structure factor.
This results in only a single exponent being present in the structure factor, which is always nonzero.
4.
No, an increase of the unit cell size cannot create new diffraction peaks (see lecture).
4. No, an increase of the unit cell size cannot create new diffraction peaks (see lecture).
## Exercise 1: Equivalence of direct and reciprocal lattice
1.
$$
V^*=\left|\mathbf{b}_{1} \cdot\left(\mathbf{b}_{2} \times \mathbf{b}_{3}\right)\right| = \frac{2\pi}{V}\left| (\mathbf{a}_{2} \times \mathbf{a}_{3}) \cdot\left(\mathbf{b}_{2} \times \mathbf{b}_{3}\right)\right| = \frac{(2\pi)^3}{V}
$$
1. Solution
In the second equality, we used the reciprocal lattice vector definition (see notes).
In the third equality, we used the identity:
$$
V^*=\left|\mathbf{b}_{1} \cdot\left(\mathbf{b}_{2} \times \mathbf{b}_{3}\right)\right| = \frac{2\pi}{V}\left| (\mathbf{a}_{2} \times \mathbf{a}_{3}) \cdot\left(\mathbf{b}_{2} \times \mathbf{b}_{3}\right)\right| = \frac{(2\pi)^3}{V}
$$
$$
(\mathbf{a} \times \mathbf{b}) \cdot(\mathbf{c} \times \mathbf{d})=(\mathbf{a} \cdot \mathbf{c})(\mathbf{b} \cdot \mathbf{d})-(\mathbf{a} \cdot \mathbf{d})(\mathbf{b} \cdot \mathbf{c})
$$
In the second equality, we used the reciprocal lattice vector definition (see notes).
In the third equality, we used the identity:
2.
$$
(\mathbf{a} \times \mathbf{b}) \cdot(\mathbf{c} \times \mathbf{d})=(\mathbf{a} \cdot \mathbf{c})(\mathbf{b} \cdot \mathbf{d})-(\mathbf{a} \cdot \mathbf{d})(\mathbf{b} \cdot \mathbf{c})
$$
Because the relation between direct and reciprocal lattice is symmetric, so are the expressions for the direct lattice vectors through the reciprocal ones:
2. Because the relation between direct and reciprocal lattice is symmetric, so are the expressions for the direct lattice vectors through the reciprocal ones:
$$
\mathbf{a}_{i} \epsilon_{ijk} = \frac{2\pi}{V^*} (\mathbf{b}_{j} \times \mathbf{b}_{k})
$$
$$
\mathbf{a}_{i} \epsilon_{ijk} = \frac{2\pi}{V^*} (\mathbf{b}_{j} \times \mathbf{b}_{k})
$$
where $\epsilon_{ijk}$ is the [Levi-Civita tensor](https://en.wikipedia.org/wiki/Levi-Civita_symbol#Three_dimensions)
where $\epsilon_{ijk}$ is the [Levi-Civita tensor](https://en.wikipedia.org/wiki/Levi-Civita_symbol#Three_dimensions)
3.
One set of the BCC primitive lattice vectors is given by:
$$
\mathbf{a_1} = \frac{a}{2} \left(-\hat{\mathbf{x}}+\hat{\mathbf{y}}+\hat{\mathbf{z}} \right) \\
\mathbf{a_2} = \frac{a}{2} \left(\hat{\mathbf{x}}-\hat{\mathbf{y}}+\hat{\mathbf{z}} \right) \\
\mathbf{a_3} = \frac{a}{2} \left(\hat{\mathbf{x}}+\hat{\mathbf{y}}-\hat{\mathbf{z}} \right).
$$
3. One set of the BCC primitive lattice vectors is given by:
From this, we find the following set of reciprocal lattice vectrs:
$$
\mathbf{a_1} = \frac{a}{2} \left(-\hat{\mathbf{x}}+\hat{\mathbf{y}}+\hat{\mathbf{z}} \right) \\
\mathbf{a_2} = \frac{a}{2} \left(\hat{\mathbf{x}}-\hat{\mathbf{y}}+\hat{\mathbf{z}} \right) \\
\mathbf{a_3} = \frac{a}{2} \left(\hat{\mathbf{x}}+\hat{\mathbf{y}}-\hat{\mathbf{z}} \right).
$$
$$
\mathbf{b_1} = \frac{2 \pi}{a} \left(\hat{\mathbf{y}}+\hat{\mathbf{z}} \right) \\
\mathbf{b_2} = \frac{2 \pi}{a} \left(\hat{\mathbf{x}}+\hat{\mathbf{z}} \right) \\
\mathbf{b_3} = \frac{2 \pi}{a} \left(\hat{\mathbf{x}}+\hat{\mathbf{y}} \right),
$$
From this, we find the following set of reciprocal lattice vectors:
which is forms a reciprocal FCC lattice.
The opposite relation follows directly from our previous result.
$$
\mathbf{b_1} = \frac{2 \pi}{a} \left(\hat{\mathbf{y}}+\hat{\mathbf{z}} \right) \\
\mathbf{b_2} = \frac{2 \pi}{a} \left(\hat{\mathbf{x}}+\hat{\mathbf{z}} \right) \\
\mathbf{b_3} = \frac{2 \pi}{a} \left(\hat{\mathbf{x}}+\hat{\mathbf{y}} \right),
$$
4.
Because the 1st Brillouin Zone is the Wigner-Seitz cell of the reciprocal lattice, we need to construct the Wigner-Seitz cell of the FCC lattice.
For visualization, it is convenient to look at [FCC lattice](https://solidstate.quantumtinkerer.tudelft.nl/9_crystal_structure/#face-centered-cubic-lattice) introduced in the previous lecture and count the neirest neighbours of each lattice point.
We see that each lattice point contains 12 neirest neighbours and thus the Wigner-Seitz cell contains 12 sides!
which forms a reciprocal FCC lattice.
The opposite relation follows directly from our previous result.
## Exercise 2: Miller planes and reciprocal lattice vectors
4. Because the 1st Brillouin Zone is the Wigner-Seitz cell of the reciprocal lattice, we need to construct the Wigner-Seitz cell of the FCC lattice.
For visualization, it is convenient to look at [FCC lattice](https://solidstate.quantumtinkerer.tudelft.nl/9_crystal_structure/#face-centered-cubic-lattice) introduced in the previous lecture and count the neirest neighbours of each lattice point.
We see that each lattice point contains 12 neirest neighbours and thus the Wigner-Seitz cell contains 12 sides!
1.
??? hint "First small hint"
## Exercise 2: Miller planes and reciprocal lattice vectors
The $(hkl)$ plane intersects lattice at position vectors of $\frac{\mathbf{a_1}}{h}, \frac{\mathbf{a_2}}{k}, \frac{\mathbf{a_3}}{l}$.
Can you define a general vector inside the $(hkl)$ plane?
1. Hints
??? hint "Second small hint"
Whats the best vector operation to show orthogonality between two vectors?
??? hint "First small hint"
The $(hkl)$ plane intersects lattice at position vectors of $\frac{\mathbf{a_1}}{h}, \frac{\mathbf{a_2}}{k}, \frac{\mathbf{a_3}}{l}$.
Can you define a general vector inside the $(hkl)$ plane?
2.
One can compute the normal to the plane by using result from Subquestion 1:
??? hint "Second small hint"
Whats the best vector operation to show orthogonality between two vectors?
2. One can compute the normal to the plane by using result from Subquestion 1:
$\hat{\mathbf{n}} = \frac{\mathbf{G}}{|G|}$
$\hat{\mathbf{n}} = \frac{\mathbf{G}}{|G|}$
Let us consider a very simple case in which we have the miller planes $(h00)$.
For lattice planes, there is always a plane intersecting the zero lattice point (0,0,0).
As such, the distance from this plane to the closest next one is given by:
Let us consider a very simple case in which we have the miller planes $(h00)$.
For lattice planes, there is always a plane intersecting the zero lattice point (0,0,0).
As such, the distance from this plane to the closest next one is given by:
$ d = \hat{\mathbf{n}} \cdot \frac{\mathbf{a_1}}{h} = \frac{2 \pi}{|G|} $
$ d = \hat{\mathbf{n}} \cdot \frac{\mathbf{a_1}}{h} = \frac{2 \pi}{|G|} $
3.
Since $\rho=d / V$, we must maximize $d$.
To do that, we minimize must $|G|$.
Therefore the smallest possible reciprocal lattice vectors are the (100) family of planes (in terms of FCC primitive lattice vectors).
3. Since $\rho=d / V$, we must maximize $d$.
To do that, we minimize must $|G|$.
Therefore the smallest possible reciprocal lattice vectors are the (100) family of planes (in terms of FCC primitive lattice vectors).
## Exercise 3: X-ray scattering in 2D
1.
1. See figure
```python
def reciprocal_lattice(N = 7, lim = 40):
y = np.repeat(np.linspace(-18.4*(N//2),18.4*(N//2),N),N)
@@ -119,11 +108,11 @@ def reciprocal_lattice(N = 7, lim = 40):
reciprocal_lattice()
```
2.
Since we have elastic scattering, we obtain
2. Since we have elastic scattering, we obtain
$|\mathbf{k}| = |\mathbf{k}'| = \frac{2 \pi}{\lambda} = 37.9 nm^{-1}$
3. See figure
3.
```python
reciprocal_lattice()
# G vector
@@ -139,37 +128,29 @@ plt.annotate('$\mathbf{k\'}$',(15,30),fontsize=14, ha='center',color='k');
## Exercise 4: Structure factors
1.
$S(\mathbf{G}) = \sum_j f_j e^{i \mathbf{G} \cdot \mathbf{r_j}} = f(1 + e^{i \pi (h+k+l)})$
2.
Solving for $h$, $k$, and $l$ results in
$
S(\mathbf{G}) = \begin{cases}
2f, \: \text{if $h+k+l$ is even}\\
0, \: \text{if $h+k+l$ is odd}.
\end{cases}
$
1. $S(\mathbf{G}) = \sum_j f_j e^{i \mathbf{G} \cdot \mathbf{r_j}} = f(1 + e^{i \pi (h+k+l)})$
2. Solving for $h$, $k$, and $l$ results in
Thus if $h+k+l$ is odd, diffraction peaks dissapear
$
S(\mathbf{G}) = \begin{cases}
2f, \: \text{if $h+k+l$ is even}\\
0, \: \text{if $h+k+l$ is odd}.
\end{cases}
$
3.
Thus if $h+k+l$ is odd, diffraction peaks dissapear
Let $f_1 \neq f_2$, then
3. Let $f_1 \neq f_2$, then
$
S(\mathbf{G}) = \begin{cases}
f_1 + f_2, \text{if $h+k+l$ is even}\\
f_1 - f_2, \text{if $h+k+l$ is odd}
\end{cases}
$
$
S(\mathbf{G}) = \begin{cases}
f_1 + f_2, \text{if $h+k+l$ is even}\\
f_1 - f_2, \text{if $h+k+l$ is odd}
\end{cases}
$
4.
Due to bcc systematic absences, the peaks from lowest to largest angle are:
4. Due to bcc systematic absences, the peaks from lowest to largest angle are:
$(110),(200),(211), (220), (310)$
5.
$a = 2.9100 \unicode{xC5}$
5. $a = 2.9100$Å