Liuyun Dao,Hui Shao,Anders W. Sandvik
Abstract:We study the excitation spectrum of the $S=1/2$ Heisenberg model on the randomly diluted square lattice by analytic continuation of QMC data. At dilution fractions $p=1/16$ and $p=1/8$, the dynamic structure factor $S({\bf q},\omega)$ exhibits a damped magnon peak with anomalous dispersion near ${\bf q}=(0,0)$ and $(\pi,\pi)$, a non-dispersive low-energy localization peak, and a second dispersive peak between these two features. A magnon with anomalous dispersion, close to our result, was predicted in spin wave and $T$-matrix theory [A. Chernyshev et al., PRB {\bf 65}, 104407 (2002)], above the localization energy. However, no intermediate dispersive mode was predicted. Analyzing spectral functions in real space for individual vacancy realizations by energy tomography, we find that these excitations are concentrated on a small subset of the spins adjacent to vacancies. We argue that the low-energy excitations are those of a sparse random network of effective moments at a fraction of the vacancies. There is a shift in magnon spectral weight distribution, from the spins away from vacancies at high energy to those adjacent to vacancies at lower energy. We also analyze the Anderson quantum rotor excitation at $\omega \propto N^{-1}$ (with $N=L^2$ the system size), which in the clean system is visible in $S({\bf q},\omega)$ only at ${\bf q}=(\pi,\pi)$ but spreads through the Brillouin zone when $p>0$. Weight close to ${\bf q}=(0,0)$ and $(\pi,\pi)$ is explained by local sublattice imbalance within a dimer-monomer model but there is also structure arising from correlated singlet fluctuations, which we demonstrate by enhancing said fluctuations with four-spin couplings. All spectral features found here should be observable by elastic neutron scattering experiments on layered quantum antiferromagnets doped with nonmagnetic impurities.
What problem does this paper attempt to address?
This paper attempts to address the problem of how the excitation spectrum of a 2D S=1/2 Heisenberg antiferromagnet changes when the lattice is randomly diluted (i.e., some lattice sites are replaced by non-magnetic impurities). Specifically, the authors study the dynamic structure factor \( S(q, \omega) \) by analyzing the stochastic analytic continuation of imaginary-time correlation functions obtained through quantum Monte Carlo simulations in the low-temperature limit. The study mainly focuses on relatively low dilution ratios (\( p = 1/16 \) and \( p = 1/8 \)), exploring the following aspects:
1. **Characteristics of the Excitation Spectrum**:
- Strongly damped magnon peaks with anomalous dispersion, especially near momenta \( q = (0,0) \) and \( (\pi, \pi) \).
- Low-energy non-dispersive peaks from localized excitations.
- A second dispersive peak within the continuum connecting the above two features.
2. **Comparison of Theoretical Predictions and Experimental Results**:
- Previous studies (such as spin-wave theory and T-matrix theory) predicted strongly damped spin waves in the low-energy region, which significantly reduce in energy near \( q = (0,0) \) and \( (\pi, \pi) \), and the quasiparticle description fails below a certain critical energy \( \omega_{\text{loc}} \).
- However, previous theories did not predict the intermediate dispersive mode, which is clearly observed in the numerical calculations of this paper.
3. **Nature of Localized Excitations**:
- Energy tomography analysis reveals that these excitations are mainly concentrated on spins adjacent to vacancies, and the number of spins involved in the excitations decreases as the energy lowers.
- It is suggested that low-energy excitations are composed of effective moments in a sparse random network.
4. **Anderson Quantum Rotor Excitations**:
- In a clean system, such excitations are only visible at \( q = (\pi, \pi) \), but in a diluted system, due to the breaking of translational symmetry, they extend throughout the Brillouin zone.
5. **Complex Structure in Momentum Space**:
- Besides the weight distribution near \( q = (0,0) \) and \( (\pi, \pi) \), a complex structure caused by correlated singlet fluctuations is also found.
- The impact of these fluctuations is further verified by enhancing the four-spin coupling term.
6. **Possibility of Experimental Observation**:
- All \( \omega > 0 \) excitations can be observed through inelastic neutron scattering experiments in layered quantum antiferromagnets doped with non-magnetic impurities.
- Particularly in copper oxides, substituting copper with zinc or magnesium can serve as candidate materials for experimental studies.
Overall, this paper reveals the complex excitation spectrum characteristics in diluted 2D Heisenberg antiferromagnets through detailed numerical calculations and theoretical analysis, providing important theoretical support for further experimental research.