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.
Developing analytical and numerical tools for strongly correlated systems is a central challenge for the condensed matter physics community. In the absence of exact solutions and controlled analytical approximations, numerical techniques have often contributed to our understanding of these systems. Exact Diagonalization (ED) requires the storage of at least two vectors the size of the Hilbert space under consideration (which grows exponentially with system size) which makes it affordable only for small systems. The Density Matrix Renormalization Group (DMRG) uses an intelligent Hilbert space truncation procedure to significantly reduce this cost, but in its present formulation is limited to quasi-1D systems. Quantum Monte Carlo (QMC) maps the Schrödinger equation to the diffusion equation (in imaginary time) and only samples the eigenvector over time, thereby avoiding the memory limitation. However, the stochasticity involved in the method gives rise to the the "sign problem" characteristic of fermion and frustrated spin systems. The first part of this thesis is an effort to make progress in the development of a numerical technique which overcomes the above mentioned problems. We consider novel variational wavefunctions, christened "Correlator Product States" (CPS), that have a general functional form which hopes to capture essential correlations in the ground states of spin and fermion systems in any dimension. We also consider a recent proposal to modify projector (Green's Function) Quantum Monte Carlo to ameliorate the sign problem for realistic and model Hamiltonians (such as the Hubbard model). This exploration led to our own set of improvements, primarily a semistochastic formulation of projector Quantum Monte Carlo. Despite their limitations, existing numerical techniques can yield physical insights into a wide variety of problems. The second part of this thesis considers one such numerical technique - DMRG - and adapts it to study the Heisenberg antiferromagnet on a generic tree graph. Our attention turns to a systematic numerical and semi-analytical study of the effect of local even/odd sublattice imbalance on the low energy spectrum of antiferromagnets on regular Cayley trees. Finally, motivated by previous experiments and theories of randomly diluted antiferromagnets (where an even/odd sublattice imbalance naturally occurs), we present our study of the Heisenberg antiferromagnet on the Cayley tree at the percolation threshold. Our work shows how to detect "emergent" low energy degrees of freedom and compute the effective interactions between them by using data from DMRG calculations.
Linear Heisenberg antiferromagnets (HAFs) are chains of spin-
S sites with isotropic exchange
J between neighbors. Open and periodic boundary conditions return the same ground-state energy per site in the thermodynamic limit, but not the same spin
SG when $S\ensuremath{\ge}1$. The ground state of open chains of
N spins has
SG=0 or
S, respectively, for even or odd
N. Density-matrix renormalization-group calculations with different algorithms for even and odd
N are presented up to
N=500 for the energy and spin densities $\ensuremath{\rho}(r,N)$ of edge states in HAFs with
S=1, 3/2, and 2. The edge states are boundary-induced spin density waves (BI-SDWs) with $\ensuremath{\rho}(r,N)\ensuremath{\propto}{(\ensuremath{-}1)}^{r\ensuremath{-}1}$ for
r=1,2,...,N. The SDWs are in phase when
N is odd, are out of phase when
N is even, and have finite excitation energy $\mathrm{\ensuremath{\Gamma}}(N)$ that decreases exponentially with
N for integer
S and faster than
1/N for half integer
S. The spin densities and excitation energy are quantitatively modeled for integer
S chains longer than $5\ensuremath{\xi}$ spins by two parameters, the correlation length $\ensuremath{\xi}$ and the SDW amplitude, with $\ensuremath{\xi}=6.048$ for
S=1 and 49.0 for
S=2. The BI-SDWs of
S=3/2 chains are not localized and are qualitatively different for even and odd
N. Exchange between the ends for odd
N is mediated by a delocalized effective spin in the middle that increases $|\mathrm{\ensuremath{\Gamma}}(N)|$ and weakens the size dependence. The nonlinear sigma model ($\mathrm{NL}\ensuremath{\sigma}\mathrm{M}$) has been applied to the HAFs, primarily to
S=1 with even
N, to discuss spin densities and exchange between localized states at the ends as $\mathrm{\ensuremath{\Gamma}}(N)\ensuremath{\propto}{(\ensuremath{-}1)}^{N}exp(\ensuremath{-}N/\ensuremath{\xi})$.
S=1 chains with odd
N are fully consistent with the $\mathrm{NL}\ensuremath{\sigma}\mathrm{M}$;
S=2 chains have two gaps $\mathrm{\ensuremath{\Gamma}}(N)$ with the same $\ensuremath{\xi}$ as predicted whose ratio is 3.45 rather than 3; the $\mathrm{NL}\ensuremath{\sigma}\mathrm{M}$ is more approximate for
S=3/2 chains with even
N and is modified for exchange between ends for odd
N.
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.
We study the properties of the ground state and the low-lying excited states of an alternating spin-1 and spin-
21 antiferromagnetic chain with ferromagnetic next-nearest-neighbor coupling with a variety of methods, namely, spin-wave analysis, density-matrix renormalization-group (DMRG) method, and exact-diagonalization method. The ground state of this model with
2N spins is ferrimagnetic with the total spin
SG=N∕2 for all parameter values. In the long chain limit, a gapless excitation is found with characteristics of goldstone mode of the ferromagnetic order with total spin $S={S}_{G}\ensuremath{-}1$. The lowest gapped excitation, however, exists in the total spin
S=SG+1 sector. Interestingly, we find that with the increase in ferromagnetic coupling, spin-wave gapped excitation becomes lower than the massless mode for the large wave numbers. Correspondingly, from DMRG analysis we find that the excitation in the total spin
S=SG+1 is lower than that in the total spin $S={S}_{G}\ensuremath{-}1$ when the size of the chain is short and the ferromagnetic interaction is large. Characteristics of these excitations are also reflected in the low temperature thermodynamic quantities, where the low-energy properties are strongly affected by the system size and the coupling strength.