We analyze real-time dynamics of the two-dimensional Bose-Hubbard model after a sudden quench starting from the Mott insulator by means of the two-dimensional tensor-network method. Calculated single-particle correlation functions are found to be in
good agreement with a recent experiment [Y. Takasu {it et al.}, Sci. Adv. {bf 6}, eaba9255 (2020)], which cross validates the experiment and the numerical simulation. By estimating the phase and group velocities from the single-particle and density-density correlation functions, we predict how these velocities vary in the moderate interaction region, which will be useful for future experiments.
We study a spin-1/2 XXZ model with a four-spin interaction on a two-leg ladder. By means of effective field theory and matrix product state calculations, we obtain rich ground-state phase diagrams that consist of eight distinct gapped phases. Four of
them exhibit spontaneous symmetry breaking with either a magnetic or valence-bond-solid (VBS) long-range order. The other four are featureless, i.e., the bulk ground state is unique and does not break any symmetry. The featureless phases include the rung singlet (RS) and Haldane phases as well as their variants, the RS* and Haldane* phases, in which twisted singlet pairs are formed between the two legs. We argue and demonstrate that Gaussian transitions with the central charge c=1 occur between the featureless phases and between the ordered phases while Ising transitions with c=1/2 occur between the featureless and ordered phases. The two types of transition lines cross at the SU(2)-symmetric point, where the criticality is described by the SU(2)_2 Wess-Zumino-Witten theory with c=3/2. The RS-Haldane* and RS*-Haldane transitions give examples of topological phase transitions. Interestingly, the RS* and Haldane* phases, which have highly anisotropic nature, appear even in the vicinity of the isotropic case. We demonstrate that all the four featureless phases are distinguished by topological indices in the presence of certain symmetries.
Motivated by recent experiments with Rydberg atoms in an optical tweezer array, we accurately map out the ground-state phase diagram of the antiferromagnetic Ising model on a square lattice with longitudinal and transverse magnetic fields using the q
uantum Monte Carlo method. For a small but nonzero transverse field, the transition longitudinal field is found to remain nearly constant. By scrutinizing the phase diagram, we uncover a narrow region where the system exhibits reentrant transitions between the disordered and antiferromagnetic phases with increasing transverse field. Our phase diagram provides a useful benchmark for quantum simulation of a Rydberg atom system.
We develop a numerical method based on matrix product states for simulating quantum many-body systems at finite temperatures without importance sampling and evaluate its performance in spin 1/2 systems. Our method is an extension of the random phase
product state (RPPS) approach introduced recently [T. Iitaka, arXiv:2006.14459]. We show that the original RPPS approach often gives unphysical values for thermodynamic quantities even in the Heisenberg chain. We find that by adding the operation of Trotter gates to the RPPS, the sampling efficiency of the approach significantly increases and its results are consistent with those of the purification approach. We also apply our method to a frustrated spin 1/2 system to exemplify that it can simulate a system in which the purification approach fails.
We study the magnetization process of the $S=1$ Heisenberg model on a two-leg ladder with further neighbor spin-exchange interaction. We consider the interaction that couples up to the next-nearest neighbor rungs and find an exactly solvable regime w
here the ground states become product states. The next-nearest neighbor interaction tends to stabilize magnetization plateaus at multiples of 1/6. In most of the exactly solvable regime, a single magnetization curve shows two series of plateaus with different periodicities.
We investigate a quantum phase transition between a Neel phase and a valence bond solid (VBS) phase, in each of which a different Z2 symmetry is broken, in a spin-1/2 two-leg XXZ ladder with a four-spin interaction. The model can be viewed as a one-d
imensional variant of the celebrated J-Q model on a square lattice. By means of variational uniform matrix product state calculations and an effective field theory, we determine the phase diagram of the model, and present evidences that the Neel-VBS transition is continuous and belongs to the Gaussian universality class with the central charge c=1. In particular, the critical exponents $beta, eta,$ and, $ u$ are found to satisfy the constraints expected for a Gaussian transition within numerical accuracy. These exponents do not detectably change along the phase boundary while they are in general allowed to do so for the Gaussian class.
Recent discovery of the half quantized thermal Hall conductivity in $alpha$-RuCl$_3$, a candidate material for the Kitaev spin liquid, suggests the presence of a highly entangled quantum state in external magnetic fields. This field induced phase app
ears between the low field zig-zag magnetic order and the high field polarized state. Motivated by this experiment, we study possible field induced quantum phases in theoretical models of the Kitaev magnets, using the two dimensional tensor network approach or infinite tensor product states. We find various quantum ground states in addition to the chiral Kitaev spin liquid occupying a small area in the phase diagram. They form a band of emergent quantum phases in an intermediate window of external magnetic fields, somewhat reminiscent of the experiment. We discuss the implications of these results in view of the experiment and previous theoretical studies.
We provide new insights into the Abelian and non-Abelian chiral Kitaev spin liquids on the star lattice using the recently proposed loop gas (LG) and string gas (SG) states [H.-Y. Lee, R. Kaneko, T. Okubo, N. Kawashima, Phys. Rev. Lett. 123, 087203 (
2019)]. Those are compactly represented in the language of tensor network. By optimizing only one or two variational parameters, accurate ansatze are found in the whole phase diagram of the Kitaev model on the star lattice. In particular, the variational energy of the LG state becomes exact(within machine precision) at two limits in the model, and the criticality at one of those is analytically derived from the LG feature. It reveals that the Abelian CSLs are well demonstrated by the short-ranged LG while the non-Abelian CSLs are adiabatically connected to the critical LG where the macroscopic loops appear. Furthermore, by constructing the minimally entangled states and exploiting their entanglement spectrum and entropy, we identify the nature of anyons and the chiral edge modes in the non-Abelian phase with the Ising conformal field theory.
There has been a great interest in magnetic field induced quantum spin liquids in Kitaev magnets after the discovery of neutron scattering continuum and half quantized thermal Hall conductivity in the material $alpha$-RuCl$_3$. In this work, we provi
de a semiclassical analysis of the relevant theoretical models on large system sizes, and compare the results to previous studies on quantum models with small system sizes. We find a series of competing magnetic orders with fairly large unit cells at intermediate magnetic fields, which are most likely missed by previous approaches. We show that quantum fluctuations are typically strong in these large unit cell orders, while their magnetic excitations may resemble a scattering continuum and give rise to a large thermal Hall conductivity. Our work provides an important basis for a thorough investigation of emergent spin liquids and competing phases in Kitaev magnets.
We provide a framework for understanding the gapless Kitaev spin liquid (KSL) in the language of tensor network(TN). Without introducing Majorana fermion, most of the features of the KSL including the symmetries, gauge structure, criticality and vort
ex-freeness are explained in a compact TN representation. Our construction reveals a hidden string gas structure of the KSL. With only two variational parameters to adjust, we obtain an accurate KSL ansatz with the bond dimension D = 8 in a compact form, where the energy is about 0.007% higher than the exact one. In addition, the opening of gap and non-Abelian phase driven by a magnetic field are naturally understood in our construction.
