We report the comparative study of the energy-level properties of the Ising, the Heisenberg, and the dipolar quantum Hamiltonians of 1/2 spins on a two-dimensional (2D) square lattice. We evaluate the number of constants of motion by the dimension of the null space of the von Neumann operator. In this way, we show that the three Hamiltonians are ordered according to: Ising---Heisenberg---dipolar, with respect to their integrability properties. Extra dynamical symmetries persist beyond the geometric symmetries in the 2D Heisenberg Hamiltonian as evidenced by a Poisson statistics for the level spacings of the Hamiltonian desymmetrized by the geometric symmetries. On the other hand, the dipolar Hamiltonian has no extra symmetries beyond the geometric ones because the desymmetrized Hamiltonian shows Wigner spacing statistics characteristic of orthogonal ensembles of random matrices.
We report the comparative study of the energy-level properties of the Ising, the Heisenberg, and the dipolar quantum Hamiltonians of 1/2 spins on a two-dimensional (2D) square lattice. We evaluate the number of constants of motion by the dimension of the null space of the von Neumann operator. In this way, we show that the three Hamiltonians are ordered according to: Ising---Heisenberg---dipolar, with respect to their integrability properties. Extra dynamical symmetries persist beyond the geometric symmetries in the 2D Heisenberg Hamiltonian as evidenced by a Poisson statistics for the level spacings of the Hamiltonian desymmetrized by the geometric symmetries. On the other hand, the dipolar Hamiltonian has no extra symmetries beyond the geometric ones because the desymmetrized Hamiltonian shows Wigner spacing statistics characteristic of orthogonal ensembles of random matrices.
Spectral properties of XXZ spin chains are studied intensively. The integrable XXZ spin chain is one of the most important quantum integrable systems, and its Hamiltonian is given by H = J L ∑ j=1 (S j S x j+1 + S y j S y j+1 +∆S z jS z j+1). (1) Here, L is the number of sites; S j = 1 2 σ j and σ α j (α = x, y, z) is the Pauli matrix on the jth site. Equation (1) corresponds to the Heisenberg spin chain at ∆ = 1, to the XY model at ∆ = 0, and to the Ising spin chain when ∆ → ∞. In this thesis, we study on XXZ spin chains with next-nearest-neighbor (NNN) coupling and/or with a magnetic field as well as Eq. (1). First, we discuss level crossings in the energy spectra of the integrable XXZ spin chain. Some of the level crossings are counterexamples of a theorem by von Neumann and Wigner: two energy levels with the same symmetry rarely cross each other in the flow of the spectrum. Multiple degeneracies appear at the particular crossing points where the symmetry of the system is enhanced in the spectral flow. The eigenvalues and the eigenvectors are obtained in an analytical form for 3-site and 4-site chains. The spectral flow of a 6-site chain is shown by numerical calculation. Not only level crossings but also “level bifurcations” are found in the flow of eigenvalues for the inhomogeneous transfer matrix of the XXZ spin chain. The “level bifurcations” are branches such as bifurcation or multiple branches in the flow of eigenvalues and not observed for the homogeneous transfer matrix. The observations are obtained by numerical calculation. We also confirm the appearance of some of the branches through the Bethe ansatz method. The behavior of some branching points are analyzed in the flow of eigenvalues for a particular inhomogeneous case. Level statistics is discussed for the XXZ spin chains with NNN coupling or with a random magnetic field. Evaluating the level-spacing distribution, the spectral rigidity and the number variance, we confirm the following conjecture: the level-spacing distribution should be described by the Poisson distribution for an integrable model and by the Wigner distribution for a nonintegrable model. The XXZ spin chains with the NNN coupling (J2) or with a random magnetic field (h) are nonintegrable, so that mostly the level statistics shows Wigner behavior. However, the level statistics even for a nonintegrable model sometimes shows Poisson-like behavior because of the finite-size effects. For the NNN coupled XXZ spin chain, the behavior of level statistics comes close to the Poissonian around J2 = 0 and/or ∆ = 1. As h/J becomes large, Poisson-like behavior appears for the level statistics of the XXZ spin chain under random magnetic field. When h/J 0, however, Poisson-like behavior appears for small L. In addition to the Poisson-like behavior, unexpected non-Wigner behavior of level statistics can appear due to mixed symmetry. When we discuss level statistics, we have to consider the whole symmetry of a system and perform desymmetrization. The NNN coupled XXZ chains are invariant under spin rotation around the z axis, translation, reflection. However, we have to consider spin reversal symmetry only when S z tot = 0. If spin reversal symmetry is not considered for S z tot = 0, the behavior of level statistics is described by almost the numerical average of Poisson and Wigner behavior. To observe the dynamical properties of the energy spectra, we investigate energy diffusion for NNN coupled XXZ spin chains in the presence of a periodically oscillating magnetic field. Diffusion coefficients (D) are found to obey the power law with respect to both the field strength (B0) and driving frequency (ω): D ∝ (B0ω) . We observe that β = 2 for “linear response regime” and β = 1 for “non-perturbative regime”. In the linear response regime, the diffusion coefficient can be calculated using the Kubo formula. On the other hand, B0ω is so large in the non-perturbative regime that the perturbation theory fails. The diffusion coefficient is smaller than that predicted by the Kubo formula and the perturbation theory fails in the non-perturbative regime. The ranges of the linear response and the non-perturbative regimes depend on J2. When J2 is small, we observe large-amplitude oscillations, which is attributed to a non-diffusive and ballistic nature of the underlying energy diffusion. In this way, we investigate the energy spectra of finite quantum XXZ spin chains in several different aspects. We hope that the results of the thesis will contribute to the development of the physics of low-dimensional systems and mesoscopic systems.
Spectra of several types of finite Heisenberg lattices have been studied using methods of statistical spectroscopy. In particular, spectral density distributions, discrete spectra generated from the spectral density distribution moments, and the distributions of spacings between the neighboring energy levels have been analyzed. The densities have been shown to be nearly Gaussian and the spectra generated from their several-moment Gram-Charlier approximations are in a good agreement with the exact ones, except for several extreme energy levels. Spacings between the neighboring levels exhibit clustering of the energy levels stronger than that given by the Poisson distribution, showing an essential influence of the point-group symmetries of the lattices.
Using an efficient one and two qubit gate simulator operating on graphical processing units, we investigate ergodic properties of a quantum Ising spin model on a two-dimensional lattice, which is periodically driven by a δ -pulsed transverse magnetic field. We consider three different dynamical properties: (i) level density, (ii) level spacing distribution of the Floquet quasienergy spectrum, and (iii) time-averaged autocorrelation function of magnetization components. Varying the parameters of the model, we found transitions between ordered (non-ergodic) and quantum chaotic (ergodic) phases, but the transitions between flat and non-flat spectral density do not correspond to transitions between ergodic and non-ergodic local observables. Even more surprisingly, we found good agreement of level spacing distribution with the Wigner surmise of random matrix theory for almost all values of parameters except where the model is essentially non-interacting, even in regions where local observables are not ergodic or where spectral density is non-flat. These findings question the versatility of the interpretation of level spacing distribution in many-body systems and stress the importance of the concept of locality.
The quantum phase transitions of dipoles confined to the vertices of two-dimensional lattices of square and triangular geometry is studied using path integral ground state quantum Monte Carlo. We analyze the phase diagram as a function of the strength of both the dipolar interaction and a transverse electric field. The study reveals the existence of a class of orientational phases of quantum dipolar rotors whose properties are determined by the ratios between the strength of the anisotropic dipole-dipole interaction, the strength of the applied transverse field, and the rotational constant. For the triangular lattice, the generic orientationally disordered phase found at zero and weak values of both dipolar interaction strength and applied field is found to show a transition to a phase characterized by net polarization in the lattice plane as the strength of the dipole-dipole interaction is increased, independent of the strength of the applied transverse field, in addition to the expected transition to a transverse polarized phase as the electric field strength increases. The square lattice is also found to exhibit a transition from a disordered phase to an ordered phase as the dipole-dipole interaction strength is increased, as well as the expected transition to a transverse polarized phase as the electric field strength increases. In contrast to the situation with a triangular lattice, on square lattices, the ordered phase at high dipole-dipole interaction strength possesses a striped ordering. The properties of these quantum dipolar rotor phases are dominated by the anisotropy of the interaction and provide useful models for developing quantum phases beyond the well-known paradigms of spin Hamiltonian models, implementing in particular a novel physical realization of a quantum rotor-like Hamiltonian that possesses an anisotropic long range interaction.
For a two-spin model which is (classically) integrable on a five-dimensional hypersurface in six-dimensional parameter space and for which level degeneracies occur exclusively (with one known exception) on four-dimensional manifolds embedded in the integrability hypersurface, we investigate the relations between symmetry, integrability, and the assignment of quantum numbers to eigenstates. We calculate quantum invariants in the form of expectation values for selected operators and monitor their dependence on the Hamiltonian parameters along loops within, without, and across the integrability hypersurface in parameter space. We find clear-cut signatures of integrability and nonintegrability in the observed traces of quantum invariants evaluated in finite-dimensional invariant Hilbert subspaces. The results support the notion that quantum integrability depends on the existence of action operators as constituent elements of the Hamiltonian.
We numerically study quantum chaos properties of long-range XXZ dipolar Hamiltonian spin systems. Two geometries are considered: (i) an open chain with 19 spins, (ii) a face-centered cubic lattice with 14 spins. Energy level-spacing distribution indicates that the three-dimensional geometry is highly chaotic, while the one-dimensional system is mildly chaotic for small chains, but has increasing chaoticity for larger chains. We also look at statistical properties of energy eigenvectors, and of one- and two-body local observables. Finally, we present some preliminary results on time-evolution, local spin dynamics and thermalization. Quantum chaos may have important implications for "scrambling" of quantum information, in both condensed matter systems and in astrophysical applications such as black holes.