1. Sep 7, 2011

mehdi86

Hello,I am doing a research related to Ising Model
in m research, evolution of a multi-qubit Ising system with the initial and final
Hamiltonian is given by:
(1) Hi=(-1/2)$\sum\sigma^{(i)}_{x}$
(2) Hf=(-1/2)$\sumhi\sigma^{(i)}_{z}$+(1/2)$\sumJ_{ij}\sigma^{(i)}_{z}\sigma^{(j)}_{z}$
(3) Hs(t)=(1-s(t))Hi+s(t)Hf
and in Hf must i>j in the second summation, and hi and Jij from{-1,0,1}. please help me how can I diagonalize the Hamiltonian (3) with (1) and (2) with random instances for 20 qubits by randomly choosing hi and Jij. mostly I have problem with multiplying two sigmas in Hf.
Thanks a lot.

2. Sep 7, 2011

DrDu

Maybe you should first repair the display of your formulas so that they can be read by everybody in the forum.
I don't think it is a good idea to diagonalize the hamiltonian (which is a matrix of size 2^20 times 2^20 or 10^6 times 10^6) if you are only interested in the time evolution of selected states.
Your initial hamiltonian is diagonal (or practically so, if you interchange the roles of sigma_x and sigma_z) and it should be easy to follow the time evolution of a single state.

3. Sep 10, 2011

mehdi86

excuse me,our Hamiltonian is:
Hi=(-1/2)$\sum\sigma^{(i)}_{x}$
Hf=(-1/2)$\sum\sigma^{(i)}_{z}$hi
+(1/2)$\sum\sigma^{(i)}_{z}*sigma^{(j)}_{z}$Jij
Moreover if we want diagonalize hs(t) for two first states(i.e. 0 &1 states), what shall we do?
thank you very much.

Last edited: Sep 10, 2011
4. Sep 10, 2011

DrDu

Finding the lowest states for fixed t is quite easy: It is easy to see that applying $\exp(-\tau H)$ on any starting function which is not orthogonal to the ground state. will converge to the eigenstate corresponding to the lowest eigenvalue in the limit tau to infinity. For small tau, you can write
$\exp(-\tau H)\approx (1-\tau H)$. So you just have to apply the last expression with a small value of tau (smaller than the norm of H) repeatedly to some starting function. Eventually you have to rescale the function every now an then to avoid it to become to small. Once you have found the ground state, you can start again from a state orthogonal to it and repeat the procedure to find the first excited state. If you want to find more states, there exist more elaborate methods. Consider a book on numerical mathematics, e.g. "Numerical recipes" is a good starting point.