J1-J2 Model

[lynkeryn]

I need to do an exact diagonalization and get the eigenvectors J1-J2 model (square lattice) with different J2 (J1 is fixed and equal to 1). Is it possible to do this using SpinED?

[lynkeryn]

Okay, I figured out how to set the correct model so that the values of self-energy coincide with work DOI:10.1051/jp1:1996236

basis:
  number_spins: 16
  symmetries: []
hamiltonian:
  name: "j1-j2 Hamiltonian"
  terms:
    - matrix: [[0.25,  0,  0,  0],
               [0, -0.25,  0.5,  0],
               [0,  0.5, -0.25,  0],
               [0,  0,  0,  0.25]]
      sites: [[0, 1], [0, 4], [1, 2], [1, 5], [2, 3], [2, 6], [3, 0], [3, 7],
              [4, 5], [4, 8], [5, 6], [5, 9], [6, 7], [6, 10], [7, 4], [7, 11],
              [8, 9], [8, 12], [9, 10], [9, 13], [10, 11], [10, 14], [11, 8], [11, 15],
              [12, 13], [12, 0], [13, 14], [13, 1], [14, 15], [14, 2], [15, 12], [15, 3]]
    - matrix: [[0.25,  0,  0,  0],
               [0, -0.25,  0.5,  0],
               [0,  0.5, -0.25,  0],
               [0,  0,  0,  0.25]]
      sites: [[0, 15], [0, 13], [0, 7], [0, 5], [1, 12], [1, 14], [1, 4], [1, 6], [2, 13], [2, 15], [2, 5], [2, 7], [3, 14], [3, 12], [3, 6], [3, 4], [4, 11], [4, 9], [5, 8], [5, 10], [6, 9], [6, 11], [7, 10], [7, 8], [8, 15], [8, 13], [9, 12], [9, 14], [10, 13], [10, 15], [11, 14], [11, 12]]
observables: []
number_vectors: 2
output: "j1_j2_000.h5"
datatype: "float32"
max_primme_block_size: 4
max_primme_basis_size: 20