We have studied the low-energy excitation spectrum of a dimerized and frustrated antiferromagnetic Heisenberg chain. We use an analytic approach, based on a description of the excitations as triplets above a strong-coupling singlet ground state. The quasiparticle spectrum is calculated by treating the excitations as a dilute Bose gas with infinite on-site repulsion. Additional singlet
(S=0) and triplet
(S=1) modes are found as two-particle bound states of the elementary triplets. We have also calculated the contributions of the elementary and collective excitations into the spin structure factor. Our results are in excellent agreement with exact diagonalizations and dimer series expansions data as long as the dimerization parameter \ensuremath{\delta} is not too small $(\ensuremath{\delta}>0.1),$ i.e., while the elementary triplets can be treated as localized objects.
In this work the ordering processes in classical frustrated Ising models and the ground-state behavior of frustrated quantum spin models are investigated. In the first part the phase transition from a high-temperature paramagnetic phase to a antiferromagnetic ground state is analyzed for the square-lattice Ising model with competing antiferromagnetic interactions between nearest- and next-nearest-neighbor spins. In the case of spatially isotropic interactions the phase transition to a collinear state is of special interest. Using Monte-Carlo simulations on finite lattices as well as field-theoretical methods for the continuous model it is shown that depending on the frustration three different sectors can be distinguished that show (i) Ising-universal behavior (ii) Ashkin-Teller-like non-universal behavior (varying critical exponents) and (iii) non-critical first-order behavior. For spatially anisotropic interactions the emergence of an intermediate incommensurate phase is verified by means of Monte-Carlo simulations which has been predicted by earlier mean-field calculations. The nature of this phase is studied analyzing spin-spin correlation functions and its ordering wave-vector reveals a strong temperature and size dependence. In the second part the influence of quantum fluctuations on classical antiferromagnetic ground states are investigated in anisotropic quantum spin-1/2 Heisenberg models with competing interactions. Using quantum Monte-Carlo simulations, exact diagonalization and perturbation theory the ground-state phase diagrams are calculated for a square lattice with nearest- and next-nearest-neighbor interactions and a honeycomb lattice with up to third-nearest-neighbor interactions whereby the fluctuations are chosen to be non-frustrating. For both cases the stability regions of classical magnetic states are identified and an intermediate phase is detected that shows no sign of long-range magnetic order. Further calculations also rule out valence-bond-crystal or topological order for this state, thus a quantum disordered ground state is stabilized in these models.
We apply the projected Hartree-Fock theory (PHF) for approximating ground states of Heisenberg spin clusters. Spin-rotational, point-group, and complex-conjugation symmetry are variationally restored from a broken-symmetry mean-field reference, where the latter corresponds to a product of local spin states. A fermionic formulation of the Heisenberg model furnishes a conceptual connection to PHF applications in quantum chemistry and detailed equations for a self-consistent field optimization of the reference state are provided. Different PHF variants are benchmarked for ground-state energies and spin-pair correlation functions of antiferromagnetic spin rings and three different polyhedra, with various values of the local spin-quantum number
s. Although PHF is not suitable to study the thermodynamic limit (where it reduces to the conventional HF results), the low computational cost and the compact wave-function representation make PHF a promising complement to existing approaches for ground states of finite spin clusters, particularly for large local spin
s and a moderately large number of sites
N. The present work may also motivate future explorations of more accurate post-PHF methods for Heisenberg spin clusters.
We study the low-energy dynamics of
S=1/2 antiferromagnetic Heisenberg clusters constructed by diluting a square lattice at vacancy concentration
p at and below the percolation threshold ${p}^{\ensuremath{\ast}}\ensuremath{\approx}0.407$. The finite-size scaling behavior of the average excitation gap, $⟨\ensuremath{\Delta}⟩\ensuremath{\sim}{L}^{\ensuremath{-}z}$, where
L is the cluster length, is obtained using quantum Monte Carlo results for an upper bound ${\ensuremath{\Delta}}^{\ensuremath{\ast}}$ to $\ensuremath{\Delta}$, derived from sum rules. At the percolation threshold, we obtain a dynamic exponent $z=3.6\ifmmode\pm\else\textpm\fi{}0.1\ensuremath{\approx}2{D}_{f}$ for clusters with singlet
(S=0) ground state. Here
Df=91/48 is the fractal dimensionality of the percolating cluster. We argue that this large dynamic exponent---roughly twice that expected for quantum-rotor excitations---is a consequence of weakly interacting localized effective magnetic moments, which form due to local sublattice imbalance. This picture is supported by an extremal-value analysis of local spectral gaps, which delivers an exponent relation (between
z and two exponents characterizing the local-gap distribution) reproduced by our simulation data. However, the average $⟨{\ensuremath{\Delta}}^{\ensuremath{\ast}}⟩$ over all clusters, which have mostly ground-state spin
S>0, scales with a smaller exponent than for the
S=0 clusters alone; $z\ensuremath{\approx}1.5{D}_{f}$. Lanczos exact diagonalization for small clusters show that typically, $S\ensuremath{\rightarrow}S\ensuremath{-}1$ in the lowest-energy excitations while the dominant spectral weight originates from $S\ensuremath{\rightarrow}S+1$ excitations. Thus, the scaling of $⟨{\ensuremath{\Delta}}^{\ensuremath{\ast}}⟩$ for clusters with ground state
S>0 does not reflect the lowest-energy excitations but the higher $S\ensuremath{\rightarrow}S+1$ excitations. This result can be understood within a valence bond picture. To further explore the scenario of localized moments, we introduce a classical dimer-monomer aggregation model to study the distribution of nearest-neighbor sites forming dimers (which are the objects used in mapping to the quantum-rotor model) and unpaired spins (monomers). The monomers are localized and, thus, effective magnetic moments should form in the spin system. We also study the lowest triplet excitation of
S=0 clusters using quantum Monte Carlo calculations in the valence bond basis. The triplet is concentrated at some of the classical monomer regions, confirming the mechanism of moment formation. The number of spins (and moment regions) affected by the excitation scales as a nontrivial power of the cluster size. For a dimer-diluted bilayer Heisenberg model with weak interlayer coupling (where the system remains N\'eel ordered), there is no sublattice imbalance. In this case we find $z\ensuremath{\approx}{D}_{f}$, consistent with quantum-rotor excitations. For a single layer at $p<{p}^{\ensuremath{\ast}}$ we find $z\ensuremath{\approx}2=D$, which indicates that the weakly interacting localized moment mechanism is valid only exactly at the percolation point. There is a crossover behavior close to the percolation point.
We outline the recent results on the ground state for a class of one- and two-dimensional frustrated quantum spin models with competing ferro(F)- and antiferromagnetic (AF) interactions. Frustrated spin systems are known to have many interesting properties due to large quantum fluctuations. As a result of these fluctuations usual mean-field approach gives quite crude (if not false) description of these systems. Therefore, exactly solvable models are very instructive at investigations of such systems. The exact ground state wave function of proposed models has a structure of the valence-bond state (VBS) type. One of the 1D model describes the transition line between the F and AF phase. The exact singlet ground state on this line has a double-spiral ordering. Using different approximation methods we study the magnetization curve in the AF phase. The second considered set of the 1D and 2D models has an exact non-degenerate ground state with exponentially decaying spin correlations. We also proposed the 1D and 2D electronic models with exact ground state represented in terms of singlet bond functions which are the generalization of the RVB functions including ionic states.