We consider a quantum many-body system on a lattice which exhibits a spontaneous symmetry breaking in its infinite volume ground states, but in which the corresponding order operator does not commute with the Hamiltonian. Typical examples are the Heisenberg antiferromagnet with a Néel order, and the Hubbard model with a (superconducting) off-diagonal long range order. In the corresponding finite system, the symmetry breaking is usually “obscured ” by “quantum fluctuation ” and one gets a symmetric ground state with a long range order. In such a situation, Horsch and von der Linden proved that the finite system has a low-lying eigenstate whose excitation energy is not more than of order N −1, where N denotes the number of the sites in the lattice. Here we study the situation where the broken symmetry is a continuous one. For a particular set of states (which are orthogonal to the ground state and with each other), we prove bounds for their energy expectation values. The bounds establish that there exist ever increasing numbers of low-lying eigenstates whose excitation energies are bounded by a constant times N −1. A crucial feature of the particular low-lying
In frustrated magnetic systems with a subextensive number of classical ground states, quantum zero-point fluctuations can select a unique long-range ordered state, a celebrated phenomenon referred to as order by quantum disorder (ObQD). For frustrated spin-1/2 models, unbiased numerical methods able to expose ObQD are necessary. We show that ObQD can be identified from exact diagonalization calculations through an analysis akin to the Anderson tower of states associated with spontaneous symmetry breaking. By defining an effective quantum rotor model, we describe the competition between ObQD-induced localization of the rotor and its tunneling between symmetry-related ground states, identifying the crossover length scale from the finite-size regime where the rotor is delocalized to the infinite system-size limit where it becomes localized. This rotor model relates the characteristic splittings in the exact diagonalization energy spectrum to the ObQD selection energy scale, providing an estimate that can be compared to spin wave calculations. We demonstrate the general applicability of this approach in one-, two-, and three-dimensional frustrated spin models that exhibit ObQD.
In magnetically ordered systems, the breaking of SU(2) symmetry in the thermodynamic limit is associated with the appearance of a special type of low-lying excitations in finite-size energy spectra, the so-called tower of states (TOS). In the present work, we numerically demonstrate that there is a correspondence between the SU(2) tower of states and the lower part of the ground-state entanglement spectrum (ES). Using state-of-the-art density matrix renormalization group (DMRG) calculations, we examine the ES of the 2D antiferromagnetic
J1-
J2 Heisenberg model on both the triangular and kagome lattice. At large ferromagnetic
J2, the model exhibits a magnetically ordered ground state. Correspondingly, its ES contains a family of low-lying levels that are reminiscent of the energy tower of states. Their behavior (level counting, finite-size scaling in the thermodynamic limit) sharply reflects TOS features, and is characterized in terms of an effective entanglement Hamiltonian that we provide. At large system sizes, TOS levels are divided from the rest by an entanglement gap. Our analysis suggests that (TOS) entanglement spectroscopy provides an alternative tool for detecting and characterizing SU(2)-broken phases using DMRG.
The ground-state parameters of the two-dimensional
S=1/2 antiferromagnetic Heisenberg model are calculated using the stochastic series expansion quantum Monte Carlo method for $L\ifmmode\times\else\texttimes\fi{}L$ lattices with
L up to
16. The finite-size results for the energy
E, the sublattice magnetization
M, the long-wavelength susceptibility ${\ensuremath{\chi}}_{\ensuremath{\perp}}(q=2\ensuremath{\pi}/L)$, and the spin stiffness ${\ensuremath{\rho}}_{s},$ are extrapolated to the thermodynamic limit using fits to polynomials in
1/L, constrained by scaling forms previously obtained from renormalization-group calculations for the nonlinear $\ensuremath{\sigma}$ model and chiral perturbation theory. The results are fully consistent with the predicted leading finite-size corrections, and are of sufficient accuracy for extracting also subleading terms. The subleading energy correction $(\ensuremath{\sim}{1/L}^{4})$ agrees with chiral perturbation theory to within a statistical error of a few percent, thus providing numerical confirmation of the finite-size scaling forms to this order. The extrapolated ground- state energy per spin is $E=\ensuremath{-}0.669437(5)$. The result from previous Green's function Monte Carlo (GFMC) calculations is slightly higher than this value, most likely due to a small systematic error originating from ``population control'' bias in GFMC. The other extrapolated parameters are
M=0.3070(3), ${\ensuremath{\rho}}_{s}=0.175(2)$, ${\ensuremath{\chi}}_{\ensuremath{\perp}}=0.0625(9)$, and the spin-wave velocity
c=1.673(7). The statistical errors are comparable with those of previous estimates obtained by fitting loop algorithm quantum Monte Carlo data to finite-temperature scaling forms. Both
M and ${\ensuremath{\rho}}_{s}$ obtained from the finite-
T data are, however, a few error bars higher than the present estimates. It is argued that the
T=0 extrapolations performed here are less sensitive to effects of neglected higher-order corrections, and therefore should be more reliable.