In recent years, several magnetic Mott insulators with strong spin-orbit coupling were suggested to be proximate to the Kitaev quantum spin liquid, whose one of the most exciting features is the fractionalization of spin excitations into itinerant Majorana fermions and static Z2 fluxes. Unfortunately, the ground states of these systems cannot be easily captured by experiment, remaining featureless to conventional local probes. Here, we propose to study the signatures of fractionalized excitations by exploiting their coupling to the lattice vibrations, dubbed magnetoelastic coupling, which arises from the fact that the interaction between spins depends on the relative distance between them. We argue that the magnetoelastic coupling can lead to the distinct modification of the phonon dynamics, which can be observed by measuring renormalized phonon spectrum, the sound attenuation, and the phonon Hall viscosity. This makes the phonon dynamics a promising tool for the characterization and identification of quantum spin liquid phases. In this work, we focus on the magnetoelastic effects in the three-dimensional Kitaev model realized on the hyperhoneycomb lattice. The hyperhoneycomb Kitaev spin liquid is particularly interesting since the strong Kitaev interaction was observed in the Kitaev magnet β-Li2IrO3, for which the spin-orbit entangled Jeff = 1/2 moments of iridium ions form precisely the hyperhoneycomb lattice.
Quantum dynamics of strongly correlated systems is a challenging problem. Although the low energy fractional excitations of one-dimensional integrable models are often well understood, exploring quantum dynamics in these systems remains challenging in the gapless regime, especially at intermediate and high energies. Based on the algebraic Bethe ansatz formalism, we study spin dynamics in a representative one-dimensional strongly correlated model, i.e., the antiferromagnetic spin-
21 XXZ chain with the Ising anisotropy, via the form-factor formulas. Various excitations at different energy scales are identified crucial to the dynamic spin structure factors under the guidance of sum rules. At small magnetic polarizations, gapless excitations dominate the low energy spin dynamics arising from the magnetic-field-induced incommensurability. In contrast, spin dynamics at intermediate and high energies is characterized by the two- and three-string states, which are multiparticle excitations based on the commensurate N\'eel ordered background. Our work is helpful for experimental studies on spin dynamics in both condensed matter and cold atom systems beyond the low energy effective Luttinger liquid theory. Based on an intuitive physical picture, we speculate that the dynamic feature at high energies due to the multiparticle antibound state excitations can be generalized to nonintegrable spin systems.
The spin fluctuations parallel to the external magnetic field in the ground state of the one-dimensional (1D)
s=21 Heisenberg antiferromagnet are dominated by a two-parameter set of collective excitations. In a cyclic chain of N sites and magnetization
0<Mz<N/2, the ground state, which contains
2Mz spinons, is reconfigured as the physical vacuum for a different species of quasiparticles, identifiable in the framework of the coordinate Bethe ansatz by characteristic configurations of Bethe quantum numbers. The dynamically dominant excitations are found to be scattering states of two such quasiparticles. For $\stackrel{\ensuremath{\rightarrow}}{N}\ensuremath{\infty},$ these collective excitations form a continuum in $(q,\ensuremath{\omega})$ space with an incommensurate soft mode. Their matrix elements in the dynamic spin structure factor ${S}_{\mathrm{zz}}(q,\ensuremath{\omega})$ are calculated directly from the Bethe wave functions for finite N. The resulting line- shape predictions for $\stackrel{\ensuremath{\rightarrow}}{N}\ensuremath{\infty}$ complement the exact results previously derived via algebraic analysis for the exact two-spinon part of ${S}_{\mathrm{zz}}(q,\ensuremath{\omega})$ in the zero-field limit. They are relevant for neutron-scattering experiments on quasi-1D antiferromagnetic compounds in a strong magnetic field.
Based on the Bethe ansatz approach and inelastic neutron scattering experiments, we reveal the evolution of confinement of many-body Bethe strings in ordered regions of the quasi-one-dimensional antiferromagnet
YbAlO3. In the antiferromagnetic phase, the spin dynamics is dominated by confined length-1 Bethe strings, whose dominancy in the high-energy branch of the excitation spectrum yields to confined length-2 Bethe strings when the material is tuned to the spin-density-wave phase. In the thermal-induced disordered region, the confinement effect disappears, and the system restores the conventional quantum integrable physics of the one-dimensional Heisenberg model. Our results establish a unified picture based on a Bethe string for the spin dynamics in different magnetic phases of
YbAlO3, and thus provide profound insight into many-body quantum magnetism.
We construct a new mean-field theory for a quantum (spin-1/2) Heisenberg antiferromagnet in one (1D) and two (2D) dimensions using a Hartree-Fock decoupling of the four-point correlation functions. We show that the solution to the self-consistency equations based on two-point correlation functions does not produce any unphysical finite-temperature phase transition, in accord with the Mermin-Wagner theorem, unlike the common approach based on the mean-field equation for the order parameter. The next-neighbor spin-spin correlation functions, calculated within this approach, reproduce closely the strong renormalization by quantum fluctuations obtained via a Bethe ansatz in 1D and a small renormalization of the classical antiferromagnetic state in 2D. The heat capacity approximates with reasonable accuracy the full Bethe ansatz result at all temperatures in 1D. In 2D, we obtain a reduction of the peak height in the heat capacity at a finite temperature that is accessible by high-order
1/T expansions.
Spin-dynamical calculations on one-dimensional systems have relied heavily on classical ($s=\ensuremath{\infty}$) theories, despite abundant evidence that quantum effects can be extremely important at low temperatures. We present a new approach to the spin dynamics of the one-dimensional isotropic
s=21 Heisenberg antiferromagnetic (HB AF) which does not involve the many-body techniques usually employed. It is based on analytic Bethe ansatz calculations of excitation energies and densities of states combined with finite-chain calculations of matrix elements. An important feature of our method is the use of rigorous selection rules and the introduction of new selection rules, which are valid for macroscopic systems in a magnetic field. We show that in zero field the dynamical two-spin correlation function ${S}_{\ensuremath{\mu}\ensuremath{\mu}}(q,\ensuremath{\omega})$ at
T=0 is governed by a two-parameter continuum of spin-wave-type excitations. In nonzero field, the longitudinal component ${S}_{\mathrm{zz}}(q,\ensuremath{\omega})$ and the transverse components ${S}_{\mathrm{xx}}(q,\ensuremath{\omega})\ensuremath{\equiv}{S}_{\mathrm{yy}}(q,\ensuremath{\omega})$ behave quite differently because they are dominated by different continua of excitations. The former is characterized by a lowest excitation branch with a zero-frequency mode moving from the zone boundary ($q=\ensuremath{\pi}$) towards the zone center (
q=0) as the field increases, whereas the latter is characterized by a lowest branch with a zero frequency mode moving from $q=0 \mathrm{to} \ensuremath{\pi}$ with increasing field. The first part of our work features an approximate analytic expression for ${S}_{\ensuremath{\mu}\ensuremath{\mu}}(q,\ensuremath{\omega})$ at zero temperature and in zero field. Although our expression is not rigorous, exact sum rules are violated only by a small amount, and good agreement exists with the few known exact results. Our studies are extended to nonzero temperatures by placing major reliance on exact finite-chain calculations. Our work was stimulated by recent neutron scattering experiments and is oriented towards experimental comparisons. Our result for the
s=21 integrated intensity is in much better agreement with neutron scattering data on Cu
Cl2\ifmmode\cdot\else\textperiodcentered\fi{}2N(
C5D5) (CPC) than the corresponding semiclassical result. Moreover, the spectral-weight distribution in ${S}_{\ensuremath{\mu}\ensuremath{\mu}}(q,\ensuremath{\omega})$ shows increasing asymmetry as $q\ensuremath{\rightarrow}\ensuremath{\pi}$, a quantum effect, again in agreement with more recent neutron scattering data. The second part of our work deals with the effects of an applied magnetic field. We extend the analytic work of Ishimura and Shiba to obtain expressions for the energies and densities of states of the various excitation continua. It is shown that these continua are expected to give rise to multiple structures in the scattering intensity. Our results appear to be in quantitative agreement with preliminary results of a neutron study in CPC in a field of 70 kOe, revealing anomalous scattering intensity peaks. Our results repeatedly demonstrate the inadequacy of classical spin-wave theory for this problem. They call for additional experimental studies on quasi-one-dimensional antiferromagnets to examine other unusual and interesting phenomena predicted by our approach.