1.1) Openscienceprize find the groundstate of a Heisenberg model.
Why the ground state of a Heisenberg model? Well, the ground state of a system is it’s lowest energy configuration, and the ground state tends to be where exciting and unexpected quantum effects flourish. For example, when the element Mercury is cooled down near its ground state (below 4 Kelvin), it has no electrical resistance. None. Zero! Near its thermal ground state, a chunk of Mercury is a superconductor. Researchers have been trying to understand what materials exhibit superconductivity and how it’s even possible ever since it was first discovered. Kagome lattice is an arrangement of atoms in 2D pattern commonly found in minerals. Figure 1 shows a basic unit cell of the lattice made up of a hexagonal star. Because of its geometry, this lattice is heavily studied in the field of condensed matter physics. When atoms in the lattice are anti-aligned by spin, the triangular shape creates a frustration (right side of figure 1), that is, the atom on the right doesn’t know what to do as it is supposed to be anti-aligned with the one on the left and top at the same time. Figure 1: On the left: Kagome lattice unit cell consisting of 12 vertices. On the right, an example of a frustrated system with atoms anti aligned by spin. This type of frustrated system is believed to be of great importance in the study of spin liquids and superconducting materials at high temperature
1.2) Heisenberg spin-1/2 model on the kagome lattice
First, we define, $$ X_{i} \equiv \bigotimes_{k=1}^{\text{|vertices|}} (\sigma^{\mu})_{k} \text{ ; where $\sigma^{\mu}=\sigma^{x}$ for $k = i$} \\ $$
where $\sigma^{x}$, $\sigma^{y}$, and $\sigma^{z}$ are the usual $2 \times 2$ Pauli matrices, and $I$ is the $2 \times 2$ Identity matrix.
The subscript $k$ simply denotes the site/position on which the corresponding operator $\sigma^{\mu}$ acts.
The $XXX$ Heisenberg model Hamiltonian for this Prize is defined as $$ \begin{equation} \boxed{H = \sum_{\langle i j \rangle}^N X_i X_{j} + Y_i Y_{j} + Z_i Z_{j}} \text{ ; where $N =$ |edges|}\\ \text{The dimension of $H$ is $2^{\text{|vertices|}} \times 2^{\text{|vertices|}}$.} \end{equation} \\ $$ where $\langle i j \rangle$ is sum over neighboring spins and $X$, $Y$, and $Z$ are defined above for every site $i$.
The arrangement we are interested in is a kagome lattice, a unique lattice structure made up of corner-sharing triangles (as shown in the unit cell below).
$$ \begin{equation} \text{$\therefore$ |vertices| $= 18$ and |edges| is $12$ for the kagome lattice structure depicted above.} \end{equation} \\ $$
We now define the Pauli-vector: $$ \begin{equation} \overrightarrow{\sigma_{i}} \equiv (X_{i}, Y_{i}, Z_{i}) \\ \therefore \boxed{H = \sum \limits_{}^{18} \overrightarrow{\sigma_{i}}.\overrightarrow{\sigma_{j}}} \\ \text{$\therefore$ dim($H$) $= 2^{\text{12}} \times 2^{\text{12}}$} \end{equation} $$
We attempt to simulate the ground state of this arrangment using the Variational Quantum Eigensolver (VQE) The exact ground state of the Kagome arrangement can be analytically obtained by first analysing the triangular sub-structure. This triangular sub-structure/graph has 3 vertices (qubits) and 3 edges.
The Hamiltonian can therefore be written as: $$ \begin{equation} \boxed{H_{3} = \sum \limits_{}^{3} \overrightarrow{\sigma_{i}}.\overrightarrow{\sigma_{j}} = \overrightarrow{\sigma_{1}}.\overrightarrow{\sigma_{2}} + \overrightarrow{\sigma_{2}}.\overrightarrow{\sigma_{3}} + \overrightarrow{\sigma_{3}}.\overrightarrow{\sigma_{1}}} \end{equation} $$
Hence, the dimension of $H_{3} = 2^{3} \times 2^{3}$.
(Note that $\overrightarrow{\sigma_{i}}.\overrightarrow{\sigma_{j}} = \overrightarrow{\sigma_{j}}.\overrightarrow{\sigma_{i}}$)
We denote the total angular moment of this triangular spin arrangment by $\overrightarrow{\sigma}$.
$$ \begin{equation} \therefore \boxed{\overrightarrow{\sigma} \equiv \sum \limits_{i=1}^{|vertices|} \overrightarrow{\sigma_{i}} = \overrightarrow{\sigma_{1}} + \overrightarrow{\sigma_{2}} + \overrightarrow{\sigma_{3}}} \end{equation} $$
Now, $$ \begin{equation} \overrightarrow{\sigma}.\overrightarrow{\sigma} = (\overrightarrow{\sigma_{1}} + \overrightarrow{\sigma_{2}} + \overrightarrow{\sigma_{3}}).(\overrightarrow{\sigma_{1}} + \overrightarrow{\sigma_{2}} + \overrightarrow{\sigma_{3}}) \equiv (\overrightarrow{\sigma})^{2} \end{equation} $$
Similarly, we define, $$ \begin{equation} (\overrightarrow{\sigma_{i}})^{2} \equiv \overrightarrow{\sigma_{i}}.\overrightarrow{\sigma_{i}} \text{ ; where $1 \leq i \leq 3$} \\ \implies \boxed{H_{3} = \frac{1}{2}\bigg[(\overrightarrow{\sigma})^{2} – (\overrightarrow{\sigma_{1}})^{2} – (\overrightarrow{\sigma_{2}})^{2} – (\overrightarrow{\sigma_{3}})^{2}\bigg]} \end{equation} $$
Now, the Spin Angular Momentum operator (dimensionless) is given by: $$ \begin{equation} \boxed{\overrightarrow{S_{i}} = \frac{1}{2} \overrightarrow{\sigma_{i}}} \end{equation} $$
Therefore, in terms of the Spin Angular Momentum operator (dimensionless): $$ \begin{equation} \implies \boxed{H_{3} = \frac{1}{2}\bigg[(2\overrightarrow{S})^{2} – (2\overrightarrow{S_{1}})^{2} – (2\overrightarrow{S_{2}})^{2} – (2\overrightarrow{S_{3}})^{2}\bigg]} \end{equation} $$
$$ \begin{equation} \text{where $\overrightarrow{S} = \sum \limits_{i=1}^{\text{|vertices|}} \overrightarrow{S_{i}} = \overrightarrow{S_{1}} + \overrightarrow{S_{2}} + \overrightarrow{S_{3}}$ (since |vertices| $= 3$)} \\ \text{$(\overrightarrow{S})^{2} \equiv \overrightarrow{S}.\overrightarrow{S}$ and $(\overrightarrow{S_{i}})^{2} \equiv \overrightarrow{S_{i}}.\overrightarrow{S_{i}}$} \end{equation} $$
$$ \begin{equation} \boxed{\therefore H_{3} = 2\bigg[(\overrightarrow{S})^{2} – (\overrightarrow{S_{1}})^{2} – (\overrightarrow{S_{2}})^{2} – (\overrightarrow{S_{3}})^{2}\bigg]} \end{equation} $$