The recursion method is used to determine the T=0 spin dynamic structure factor Szz(q,ω) in the Luttinger liquid state and in the phase-separated state of the one-dimensional t–J model. As the exchange coupling increases from zero, the dispersions and line shapes of the dominant spin excitations are observed to undergo a major metamorphosis between the free-fermion limit and the onset of phase separation. The familiar two-spinon spectrum of the Heisenberg antiferromagnetic chain emerges gradually in the strongly phase-separated state.
The recursion method is used to determine the T=0 spin dynamic structure factor Szz(q,ω) in the Luttinger liquid state and in the phase-separated state of the one-dimensional t–J model. As the exchange coupling increases from zero, the dispersions and line shapes of the dominant spin excitations are observed to undergo a major metamorphosis between the free-fermion limit and the onset of phase separation. The familiar two-spinon spectrum of the Heisenberg antiferromagnetic chain emerges gradually in the strongly phase-separated state.
The impact of the spin-flip terms on the (static and dynamic) charge and spin correlations in the Luttinger-liquid ground state of the one-dimensional (1D) t-J model is assessed by comparison with the same quantities in the 1D t-Jz model, where spin-flip terms are absent. We employ the recursion method combined with a weak-coupling or a strong-coupling continued-fraction analysis. At Jz/t=0+ we use the Pfaffian representation of dynamic spin correlations. The changing nature of the dynamically relevant charge and spin excitations on approach of the transition to phase separation is investigated in detail. At the transition point, the t-Jz ground state has zero (static) charge correlations and very short-ranged (static) spin correlations, whereas the t-J ground state is critical. The t-Jz charge excitations (but not the spin excitations) at the transition have a single-mode nature, whereas charge and spin excitations have a complicated structure in the t-J model. A major transformation of the t-J spin excitations takes place between two distinct regimes within the Luttinger-liquid phase, while the t-Jz spin excitations are found to change much more gradually. In the t-Jz model, phase separation is accompanied by N\'eel long-range order, caused by the condensation of electron clusters with an already existing alternating up-down spin configuration (topological long-range order). In the t-J model, by contrast, the spin-flip processes in the exchange coupling are responsible for continued strong spin fluctuations (dominated by two-spinon excitations) in the phase-separated state.
The exact two-spinon part of the dynamic spin structure factor ${S}_{\mathrm{xx}}(Q,\ensuremath{\omega})$ for the one-dimensional s=1/2,XXZ model at T=0 in the antiferromagnetically ordered phase is calculated using recent advances in the algebraic analysis based on (infinite-dimensional) quantum group symmetries of this model and the related vertex models. The two-spinon excitations form a two-parameter continuum consisting of two partly overlapping sheets in $(Q,\ensuremath{\omega})$ space. The spectral threshold has a smooth maximum at the Brillouin zone boundary $(Q=\ensuremath{\pi}/2)$ and a smooth minimum with a gap at the zone center (Q=0). The two-spinon density of states has square-root divergences at the lower and upper continuum boundaries. For the two-spinon transition rates, the two regimes $0<~Q<{Q}_{\ensuremath{\kappa}}$ (near the zone center) and ${Q}_{\ensuremath{\kappa}}<Q<~\ensuremath{\pi}/2$ (near the zone boundary) must be distinguished, where ${Q}_{\ensuremath{\kappa}}\ensuremath{\rightarrow}0$ in the Heisenberg limit and ${Q}_{\ensuremath{\kappa}}\ensuremath{\rightarrow}\ensuremath{\pi}/2$ in the Ising limit. In the regime ${Q}_{\ensuremath{\kappa}}<Q<~\ensuremath{\pi}/2,$ the two-spinon transition rates relevant for ${S}_{\mathrm{xx}}(Q,\ensuremath{\omega})$ are finite at the lower boundary of each sheet, decrease monotonically with increasing \ensuremath{\omega}, and approach zero linearly at the upper boundary. The resulting two-spinon part of ${S}_{\mathrm{xx}}(Q,\ensuremath{\omega})$ is then square-root divergent at the spectral threshold and vanishes in a square-root cusp at the upper boundary. In the regime $0<{Q}_{\ensuremath{\kappa}}<~\ensuremath{\pi}/2,$ in contrast, the two-spinon transition rates have a smooth maximum inside the continuum and vanish linearly at either boundary. In the associated two-spinon line shapes of ${S}_{\mathrm{xx}}(Q,\ensuremath{\omega}),$ the linear cusps at the continuum boundaries are replaced by square-root cusps. Existing perturbation studies have been unable to capture the physics of the regime ${Q}_{\ensuremath{\kappa}}<Q<~\ensuremath{\pi}/2.$ However, their line-shape predictions for the regime $0<~Q<{Q}_{\ensuremath{\kappa}}$ are in good agreement with the exact results if the anisotropy is very strong. For weak anisotropies, the exact line shapes are more asymmetric.
Spin-1/2 chains demonstrate some of the striking effects of interactions and quantum fluctuations in one-dimensional systems. The XXZ model has been used to study the unusual properties of anisotropic spin chains in an external magnetic field. The zero temperature phase diagram for this model exhibits a critical or quasi-long-range-ordered phase which is a realization of a Luttinger liquid. While many static properties of spin-1/2 chains have been explained by combinations of analytical techniques such as bosonization and Bethe ansatz, the standard approach fails in the calculation of some time-dependent correlation functions. I present a study of the longitudinal dynamical structure factor for the XXZ model in the critical regime. I show that an approximation for the line shape of the dynamical structure factor in the limit of small momentum transfer can be obtained by going beyond the Luttinger model and treating irrelevant operators associated with band curvature effects. This approach is able to describe the width of the on-shell peak and the high-frequency tail at finite magnetic field. Integrability is shown to affect the low-energy effective model at zero field, with consequences for the line shape. The power-law singularities at the thresholds of the particle-hole continuum are investigated using an analogy with the X-ray edge problem. Using methods of Bethe ansatz and conformal field theory, I compute the exact exponents for the edge singularities of the dynamical structure factor. The same methods are used to study the long-time asymptotic behavior of the spin self-correlation function, which is shown to be dominated by a high-energy excitation.
Magnon bound states emerging in one-dimensional (1D) spin chains still lack a rigorous understanding. In this Rapid Communication we show that the length-1 spin strings significantly dominate the critical properties of spinons, magnons, and free fermions in the 1D antiferromagnetic spin-1/2 chain. Using the Bethe ansatz solution, we analytically calculate the scaling functions of the thermal and magnetic properties of the model, providing a rigorous understanding of the quantum criticality of spinons. It turns out that the double maxima in specific heat elegantly mark two crossover temperatures fanning out from the critical point, indicating three quantum phases: the Tomonaga-Luttinger liquid (TLL), the quantum critical, and fully polarized ferromagnetic phases. For the TLL phase, the Wilson ratio RW=4Ks remains almost temperature independent, where Ks is the Luttinger parameter. Furthermore, by applying our results, we precisely determine the quantum scalings and critical exponents of all magnetic properties in the ideal 1D spin-1/2 antiferromagnet Cu(C4H4N2)(NO3)2, recently studied by Kono et al. [Phys. Rev. Lett. 114, 037202 (2015)]. We further find that the magnetization maximum used in experiments is not a good quantity to map out the finite-temperature TLL phase boundary.
The exact expression derived by Bougourzi, Couture, and Kacir for the two-spinon contribution to the dynamic spin structure factor S(q,\ensuremath{\omega}) of the one-dimensional s=1/2 Heisenberg antiferromagnet at T=0 is evaluated for direct comparison with finite-chain transition rates (N\ensuremath{\leqslant}28) and an approximate analytical result previously inferred from finite-N data, sum rules, and Bethe ansatz calculations. The two-spinon excitations account for 72.89% of the total intensity in S(q,\ensuremath{\omega}). The singularity structure of the exact result is determined analytically and its spectral-weight distribution evaluated numerically over the entire range of the two-spinon continuum. The leading singularities of the frequency-dependent spin autocorrelation function, static spin structure factor, and q dependent susceptibility are determined via sum rules. The impact of the non-two-spinon excitations on the integrated intensity, the susceptibility, the frequency moments, and the Euclidian time representation of S(q,\ensuremath{\omega}) is studied on the basis of finite-size data.
The physics of the one-dimensional t-J model may be determined at J/t\ensuremath{\rightarrow}0 (i.e., the U/t\ensuremath{\rightarrow}\ensuremath{\infty} limit of the repulsive Hubbard model) and at J/t=2 with the use of the Bethe ansatz. To get a full understanding of the charge and spin correlation functions for all values of J/t, one has to resort to numerical methods. We have used two numerical methods (the world-line quantum Monte Carlo algorithm and the quantum transfer-matrix algorithm) to get information on charge and spin static structure factors, pairing correlations, and critical exponents. The (1+1)-dimensional classical system corresponding to the one-dimensional t-J model, on which the Monte Carlo technique is used, is identified to a fifteen-vertex model. We show that the t-J model undergoes phase separation at large values of J/t. Before phase separation, both critical exponents and pairing correlations in the extended s-wave channel favor the onset of superconductivity. For low values of J/t, the model shows a U/t\ensuremath{\rightarrow}\ensuremath{\infty} Hubbard-like character.