benchmarks.wiki / Public workspace
CritPt / Challenge_6_main / For twisted bilayer MoTe_2, we can build the following simplified continuum…
Problem
Answer published by the source. Consult the official source to check your work against its answer.
problem description
# Problem setup:
For twisted bilayer MoTe
$_2$, we can build the following simplified continuum model at valley:
\begin{equation}
\mathcal{H} = \int d^2 r\ ( c^\dagger_{\boldsymbol{r},b}, c^\dagger_{\boldsymbol{r},t})\left(
\begin{array}{cc}
\frac{\hbar^2 \nabla^2}{2 m^*} +2 V \sum_{i=1}^3 \cos(\boldsymbol{g}_i\cdot \boldsymbol{r}- \,\psi) & w \sum_{i=1}^3 \,e^{-i\,\boldsymbol{q}_i\cdot \boldsymbol{r}} \\
w \sum_{i=1}^3 \,e^{i\,\boldsymbol{q}_i\cdot \boldsymbol{r}} & \frac{\hbar^2 \nabla^2}{2 m^*} + 2 V \sum_{i=1}^3 \cos(\boldsymbol{g}_i\cdot \boldsymbol{r} + \,\psi)
\end{array}
\right) \left(\begin{matrix} c_{\boldsymbol{r},b} \\ c_{\boldsymbol{r},t}\end{matrix}\right)\ ,
\end{equation}
where , with the three-fold rotation symmetry, , ,
\begin{equation}
a_M = \frac{a_0}{2 \sin\left( \frac{\theta}{2} \right)} \ ,
\end{equation}
is the twist angle, and $a_0=3.52 \text{\AA }$ is the lattice constant of the monolayer MoTe$_2$.
We choose $m^* = 0.6 m_e$ with the mass of an electron, meV, , and meV.
Consider a generic isolated set of bands with projector constructed by the periodic part of the Bloch states.
We can define the quantum metric as
\begin{equation}
g_{ij}(\boldsymbol{k}) = \frac{1}{2}\mathrm{Tr}[\partial_{k_i} P_{\boldsymbol{k}} \partial_{k_j} P_{\boldsymbol{k}}]\ .
\end{equation}
The gauge-invariant part of the Wannier spread of the isolated set of bands is proportional to
\begin{equation}
\mathop{\mathrm{Tr}}\mathcal{G} = \int d^2 k\ \mathop{\mathrm{Tr}}[g(\boldsymbol{k})]\ ,
\end{equation}
where the integration ranges over the first Brillouin zone.
# Main problem:
For , what are the Chern numbers of the top three bands of the model, respectively?
The answer can have a global sign freedom owing to the definition of Chern number.
Numerically evaluate $\mathop{\mathrm{Tr}}\mathcal{G}$ for the top electron band (to two decimal places) using the following conventions.
The set of Bloch momenta is
\begin{equation}
\{ (l_1/L-1/2) \boldsymbol{b}_1 + (l_2/L-1/2) \boldsymbol{b}_2 | l_1,l_2 = 0,1,2,...,L-1 \},
\end{equation}
with , , and .
Choose $\hbar/(2 m_e) = 7619.96423 \text{ meV} \cdot \AA^2$.
Choose the following Fourier transformation convention
\begin{equation}
c^\dagger_{\boldsymbol{r},l} = \frac{1}{\sqrt{ \mathcal{V}}} \sum_{\boldsymbol{k},\boldsymbol{Q}} e^{-\mathrm{i} (\boldsymbol{k}-\boldsymbol{Q})\cdot \boldsymbol{r} } c^\dagger_{\boldsymbol{k}-\boldsymbol{Q},l}\ ,
\end{equation}
where is the volume of the whole sample.
Choose all for to satisfy (i) is a reciprocal lattice vector and (ii) .
Choose all for to satisfy (i) is a reciprocal lattice vector and (ii) .Plain-text mathematical notation (without MathML)
# Problem setup:
For twisted bilayer MoTe$_2$, we can build the following simplified continuum model at K valley:
\begin{equation}
\mathcal{H} = \int d^2 r\ ( c^\dagger_{\boldsymbol{r},b}, c^\dagger_{\boldsymbol{r},t})\left(
\begin{array}{cc}
\frac{\hbar^2 \nabla^2}{2 m^*} +2 V \sum_{i=1}^3 \cos(\boldsymbol{g}_i\cdot \boldsymbol{r}- \,\psi) & w \sum_{i=1}^3 \,e^{-i\,\boldsymbol{q}_i\cdot \boldsymbol{r}} \\
w \sum_{i=1}^3 \,e^{i\,\boldsymbol{q}_i\cdot \boldsymbol{r}} & \frac{\hbar^2 \nabla^2}{2 m^*} + 2 V \sum_{i=1}^3 \cos(\boldsymbol{g}_i\cdot \boldsymbol{r} + \,\psi)
\end{array}
\right) \left(\begin{matrix} c_{\boldsymbol{r},b} \\ c_{\boldsymbol{r},t}\end{matrix}\right)\ ,
\end{equation}
where g₁=(4π)/(√(3)a_(M))(1,0)^(T), g_(i)=C₃^(i−1)g₁ with C₃ the three-fold rotation symmetry, q₁=|g₁|(0,1/√(3))^(T), q_(i)=C₃^(i−1)q₁,
\begin{equation}
a_M = \frac{a_0}{2 \sin\left( \frac{\theta}{2} \right)} \ ,
\end{equation}
θ is the twist angle, and $a_0=3.52 \text{\AA }$ is the lattice constant of the monolayer MoTe$_2$.
We choose $m^* = 0.6 m_e$ with m_(e) the mass of an electron, V=16.5 meV, ψ=−105.9^(∘), and w=−18.8 meV.
Consider a generic isolated set of N bands with projector P_(k) constructed by the periodic part of the Bloch states.
We can define the quantum metric as
\begin{equation}
g_{ij}(\boldsymbol{k}) = \frac{1}{2}\mathrm{Tr}[\partial_{k_i} P_{\boldsymbol{k}} \partial_{k_j} P_{\boldsymbol{k}}]\ .
\end{equation}
The gauge-invariant part of the Wannier spread of the isolated set of bands is proportional to
\begin{equation}
\mathop{\mathrm{Tr}}\mathcal{G} = \int d^2 k\ \mathop{\mathrm{Tr}}[g(\boldsymbol{k})]\ ,
\end{equation}
where the integration ranges over the first Brillouin zone.
# Main problem:
For θ=3.5^(∘), what are the Chern numbers of the top three bands of the model, respectively?
The answer can have a global sign freedom owing to the definition of Chern number.
Numerically evaluate $\mathop{\mathrm{Tr}}\mathcal{G}$ for the top electron band (to two decimal places) using the following conventions.
The set of Bloch momenta k is
\begin{equation}
\{ (l_1/L-1/2) \boldsymbol{b}_1 + (l_2/L-1/2) \boldsymbol{b}_2 | l_1,l_2 = 0,1,2,...,L-1 \},
\end{equation}
with L=60, b₁=g₁, and b₂=g₁+g₂.
Choose $\hbar/(2 m_e) = 7619.96423 \text{ meV} \cdot \AA^2$.
Choose the following Fourier transformation convention
\begin{equation}
c^\dagger_{\boldsymbol{r},l} = \frac{1}{\sqrt{ \mathcal{V}}} \sum_{\boldsymbol{k},\boldsymbol{Q}} e^{-\mathrm{i} (\boldsymbol{k}-\boldsymbol{Q})\cdot \boldsymbol{r} } c^\dagger_{\boldsymbol{k}-\boldsymbol{Q},l}\ ,
\end{equation}
where V is the volume of the whole sample.
Choose all Q for c_(k−Q,t)^(†) to satisfy (i) Q−q₁ is a reciprocal lattice vector and (ii) |Q|<4.1|b₁|.
Choose all Q for c_(k−Q,b)^(†) to satisfy (i) Q+q₁ is a reciprocal lattice vector and (ii) |Q|<4.1|b₁|.Original LaTeX notation
# Problem setup:
For twisted bilayer MoTe$_2$, we can build the following simplified continuum model at $\text{K}$ valley:
\begin{equation}
\mathcal{H} = \int d^2 r\ ( c^\dagger_{\boldsymbol{r},b}, c^\dagger_{\boldsymbol{r},t})\left(
\begin{array}{cc}
\frac{\hbar^2 \nabla^2}{2 m^*} +2 V \sum_{i=1}^3 \cos(\boldsymbol{g}_i\cdot \boldsymbol{r}- \,\psi) & w \sum_{i=1}^3 \,e^{-i\,\boldsymbol{q}_i\cdot \boldsymbol{r}} \\
w \sum_{i=1}^3 \,e^{i\,\boldsymbol{q}_i\cdot \boldsymbol{r}} & \frac{\hbar^2 \nabla^2}{2 m^*} + 2 V \sum_{i=1}^3 \cos(\boldsymbol{g}_i\cdot \boldsymbol{r} + \,\psi)
\end{array}
\right) \left(\begin{matrix} c_{\boldsymbol{r},b} \\ c_{\boldsymbol{r},t}\end{matrix}\right)\ ,
\end{equation}
where $\boldsymbol{g}_1 = \frac{4 \pi}{\sqrt{3} a_{M}} (1,0)^T$, $\boldsymbol{g}_i = C_3^{i-1} \boldsymbol{g}_1$ with $C_3$ the three-fold rotation symmetry, $\boldsymbol{q}_1 = |\boldsymbol{g}_1| (0, 1/\sqrt{3})^T $, $\boldsymbol{q}_i = C_3^{i-1} \boldsymbol{q}_1$,
\begin{equation}
a_M = \frac{a_0}{2 \sin\left( \frac{\theta}{2} \right)} \ ,
\end{equation}
$\theta$ is the twist angle, and $a_0=3.52 \text{\AA }$ is the lattice constant of the monolayer MoTe$_2$.
We choose $m^* = 0.6 m_e$ with $m_e$ the mass of an electron, $V=16.5$ meV, $\psi = -105.9^\circ$, and $w = -18.8$ meV.
Consider a generic isolated set of $N$ bands with projector $P_{\boldsymbol{k}}$ constructed by the periodic part of the Bloch states.
We can define the quantum metric as
\begin{equation}
g_{ij}(\boldsymbol{k}) = \frac{1}{2}\mathrm{Tr}[\partial_{k_i} P_{\boldsymbol{k}} \partial_{k_j} P_{\boldsymbol{k}}]\ .
\end{equation}
The gauge-invariant part of the Wannier spread of the isolated set of bands is proportional to
\begin{equation}
\mathop{\mathrm{Tr}}\mathcal{G} = \int d^2 k\ \mathop{\mathrm{Tr}}[g(\boldsymbol{k})]\ ,
\end{equation}
where the integration ranges over the first Brillouin zone.
# Main problem:
For $\theta = 3.5^\circ$, what are the Chern numbers of the top three bands of the model, respectively?
The answer can have a global sign freedom owing to the definition of Chern number.
Numerically evaluate $\mathop{\mathrm{Tr}}\mathcal{G}$ for the top electron band (to two decimal places) using the following conventions.
The set of Bloch momenta $\boldsymbol{k}$ is
\begin{equation}
\{ (l_1/L-1/2) \boldsymbol{b}_1 + (l_2/L-1/2) \boldsymbol{b}_2 | l_1,l_2 = 0,1,2,...,L-1 \},
\end{equation}
with $L=60$, $\boldsymbol{b}_1 = \boldsymbol{g}_1$, and $\boldsymbol{b}_2 = \boldsymbol{g}_1 + \boldsymbol{g}_2$.
Choose $\hbar/(2 m_e) = 7619.96423 \text{ meV} \cdot \AA^2$.
Choose the following Fourier transformation convention
\begin{equation}
c^\dagger_{\boldsymbol{r},l} = \frac{1}{\sqrt{ \mathcal{V}}} \sum_{\boldsymbol{k},\boldsymbol{Q}} e^{-\mathrm{i} (\boldsymbol{k}-\boldsymbol{Q})\cdot \boldsymbol{r} } c^\dagger_{\boldsymbol{k}-\boldsymbol{Q},l}\ ,
\end{equation}
where $\mathcal{V}$ is the volume of the whole sample.
Choose all $\boldsymbol{Q}$ for $c^\dagger_{\boldsymbol{k}-\boldsymbol{Q},t}$ to satisfy (i) $\boldsymbol{Q}-\boldsymbol{q}_1$ is a reciprocal lattice vector and (ii) $|\boldsymbol{Q}|< 4.1 |\boldsymbol{b}_1|$.
Choose all $\boldsymbol{Q}$ for $c^\dagger_{\boldsymbol{k}-\boldsymbol{Q},b}$ to satisfy (i) $\boldsymbol{Q}+\boldsymbol{q}_1$ is a reciprocal lattice vector and (ii) $|\boldsymbol{Q}|< 4.1 |\boldsymbol{b}_1|$.code template
Code
def answer():
r"""
Return the value of the Chern numbers of the top three bands and the gauge–invariant Wannier spread TrG.
Inputs
----------
None
Outputs
----------
chern_numbers: tuple[int, int, int]
(C1, C2, C3) – Chern numbers of the first, second and third top bands, respectively, up to an overall sign.
TrG: float
Gauge–invariant part of the Wannier spread for the top electron band, $\mathop{\mathrm{Tr}}\mathcal{G}$, rounded to 2 decimal places.
"""
# ------------------ FILL IN YOUR RESULTS BELOW ------------------
chern_numbers = ... # three integers
TrG = ... # float rounded to 2 decimal places
# ---------------------------------------------------------------
return chern_numbers, TrGDiscussion
No discussion posts on this page yet. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.
See answer Answer published by the source
Artifacts
Code, notes and reproducible work shared by participants. Files are served from a separate origin.
No artifacts on this page yet. Share reproducible code or notes in a contribution. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.
Source and history
initial import