## XXZChain antiferromagnet showing no entanglement?

How do I use this algorithm? What does that parameter do?
Posts: 1
Joined: 25 Nov 2019, 04:12

### XXZChain antiferromagnet showing no entanglement?

Hi, I just started using TeNPy. I was trying to simulate the XXZChain in Google Colab with parameters, Jxx = 1. , Jz = 2., hz = 0 with tebd.Engine but it keeps giving back zero in the entanglement_spectrum. Here is my code:

Code: Select all

tebd_params = {'order': 2, 'delta_tau_list': [0.1, 0.001, 1.e-5],
'max_error_E': 1.e-6,
'trunc_params': {'chi_max': 50, 'svd_min': 1.e-10}}

M = XXZChain({'L': 100, 'Jxx': 1., 'Jz': 2., 'hz': 0., 'bc_MPS': 'finite'})

product_state = ['up'] * M.lat.N_sites
psi = MPS.from_product_state(M.lat.mps_sites(), product_state, bc = M.lat.bc_MPS)
eng = tebd.Engine(psi, M, tebd_params)
eng.run_GS()

psi.entanglement_spectrum(False)

What am I doing wrong here? Thanks in advance .

Johannes
The all-up state is an exact eigenstate of the model and hence won't evolve in time, in particular if you use Sz-conservation (which you do in this case): the all-up state is the only state in the sector with $$S_{total} = L/2$$.
product_state = ['up', 'down'] * (M.lat.N_sites //2)  # assuming that M.lat.N_sites is odd