I would like to find the exact eigenvalues of the following tight-binding Hamiltonian, written here in second quatization:

\begin{eqnarray} \hspace{-0.25in}{\mathcal{H}} &=& \mathcal{H}_0+ \frac{1}{2} \kappa a^2 \sum_{m, n} (m^2+n^2) \hat{a}^\dagger_{m,n} \hat{a}_{m,n} , \nonumber \\ \hspace{-0.25in}\mathcal{H}_0 &=& - J \sum_{m, n} \left( \hat{a}^\dagger_{m+1,n} \hat{a}_{m,n}+ e^{i \phi } \hat{a}^\dagger_{m, n+1} \hat{a}_{m,n} \right) + \mbox{h.c.}\label{eq:hof} \end{eqnarray}

How can I do so numerically? I know how to diagonalize the Hamiltonian matrix, but the question is how to get to that in the first place, i.e. how do I compute the matrix elements?

Edit: to clarify, I am only interested in finding the eigenvalues in the subspace of single-particle states.


You can try just working out the matrix elements between the different localised states, i.e. $a_1^{\dagger}\ldots a_m^{\dagger}|0\rangle$ etc. This may be rather tedious and require you to carefully keep track of phases due to the fermionic anti-symmetry. Probably a better way is to use the Jordan-Wigner transformation, which gives you an explicit matrix representation of the fermion operators. In 2D this representation of the Hamiltonian looks ugly, but it is still perfectly useful. You will be restricted to a fairly small lattice size, i.e. 10-20 sites depending on your computer and the efficiency of your code. It will be more efficient if you use the fermion number conservation to diagonalise within each number sector separately.

In the particularly simple case of a single particle on the lattice, the problem becomes rather easy, since exchange statistics are no longer relevant. In this case, you make the replacement $$a_{i,j}^{\dagger}a_{m,n} \to |i,j\rangle\langle m,n|,$$ where $$|i,j\rangle = a_{i,j}^{\dagger}|0\rangle.$$

  • $\begingroup$ Well I want to diagonalize this on a 2D lattice of at least 40x40 sites, reproducing the results of arxiv.org/abs/1403.6041v3 $\endgroup$ – Andrei Jan 14 '15 at 22:08
  • $\begingroup$ I basically just need an expression for $\mathcal{H}_{i,j}$ $\endgroup$ – Andrei Jan 14 '15 at 22:10
  • $\begingroup$ @Andrei Oh, but there they are just looking at the one-particle subspace. See my edit to the answer. $\endgroup$ – Mark Mitchison Jan 14 '15 at 22:15
  • $\begingroup$ Yes, its single-particle physics. However I still don't see how to get the matrix I need to diagonalize. $\endgroup$ – Andrei Jan 14 '15 at 22:21
  • $\begingroup$ @Andrei Well, you find a unique number $x_{ij}$ to label each site according to its coordinates $(i,j)$. Then the entry in row $x_{ij}$ and column $x_{mn}$ is the coefficient of $|i,j\rangle\langle m,n|$. $\endgroup$ – Mark Mitchison Jan 14 '15 at 22:31

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.