"Numerics for an XXZ spin chain out of equilibrium - revisiting Antal's quench"

We investigate the dynamics of a critical XXZ spin-1/2 chain prepared in an inhomogeneous initial state with different magnetizations on the left and right halves. We simulate the real-time evolution using the time-evolving block decimation (TEBD) method. We follow the front propagation by measuring the magnetization and entanglement entropy profiles, and we focus on the situation where the initial state is not fully polarized. At long times, as in the free fermion case (T. Antal et al. 1999), a large central region develops where correlations become time-independent and translation invariant. The shape and speed of the fronts is studied numerically and we evaluate the stationary current as a function of initial magnetic field and as a function of the anisotropy \Delta. We compare the results with the conductance of a Tomonaga-Luttinger liquid, and with the exact free-fermion solution at \Delta=0. We also investigate the two-point correlations in the stationary region and find a good agreement with the "twisted" form obtained by J. Lancaster and A. Mitra (2010) using bosonization. Some deviations are nevertheless observed for strong currents.