The honeycomb antiferromagnet under a triaxial strain is studied using the quantum Monte Carlo simulation. The strain dimerizes the exchange couplings near the corners, thus destructs the antiferromagnetic order therein. The antiferromagnetic region
is continuously reduced by the strain. For the same strain strength, the exact numerical results give a much smaller antiferromagnetic region than the linear spin-wave theory. We then study the strained $XY$ antiferromagnet, where the magnon pseudo-magnetic field behaves quite differently. The $0$th Landau level appears in the middle of the spectrum, and the quantized energies above (below) it are proportional to $n^{frac{1}{3}} (n^{frac{2}{3}})$, which is in great contrast to the equally-spaced ones in the Heisenberg case. Besides, we find the antiferromagnetic order of the $XY$ model is much more robust to the dimerization than the Heisenberg one. The local susceptibility of the Heisenberg case is extracted by the numerical analytical continuation, and no sign of the pseudo-Landau levels is resolved. It is still not sure whether the result is due to the intrinsic problem of the numerical analytical continuation. Thus the existence of the magnon pseudo-Landau levels in the spin-$frac{1}{2}$ strained Heisenberg Hamiltonian remains an open question. Our results are closely related to the two-dimensional van der Waals quantum antiferromagnets and may be realized experimentally.
We study scaling behavior of the disorder parameter, defined as the expectation value of a symmetry transformation applied to a finite region, at the deconfined quantum critical point in (2+1)$d$ in the $J$-$Q_3$ model via large-scale quantum Monte C
arlo simulations. We show that the disorder parameter for U(1) spin rotation symmetry exhibits perimeter scaling with a logarithmic correction associated with sharp corners of the region, as generally expected for a conformally-invariant critical point. However, for large rotation angle the universal coefficient of the logarithmic corner correction becomes negative, which is not allowed in any unitary conformal field theory. We also extract the current central charge from the small rotation angle scaling, whose value is much smaller than that of the free theory.
Recent experiments [J. Guo et al., Phys. Rev. Lett.124,206602 (2020)] on thermodynamic properties of the frustrated layered quantum magnet SrCu$_2$(BO$_3$)$_2$ -- the Shastry-Sutherland material -- have provided strong evidence for a low-temperature
phase transition between plaquette-singlet and antiferromagnetic order as a function of pressure. Further motivated by the recently discovered unusual first-order quantum phase transition with an apparent emergent O(4) symmetry of the antiferromagnetic and plaquette-singlet order parameters in a two-dimensional checkerboard J-Q quantum spin model [B. Zhao et al., Nat. Phys. 15, 678 (2019)], we here study the same model in the presence of weak inter-layer couplings. Our focus is on the evolution of the emergent symmetry as the system crosses over from two to three dimensions and the phase transition extends from strictly zero temperature in two dimensions up to finite temperature as expected in SrCu$_2$(BO$_3$)$_2$. Using quantum Monte Carlo simulations, we map out the phase boundaries of the plaquette-singlet and antiferromagnetic phases, with particular focus on the triple point where these two order phases meet the paramagnetic phase for given strength of the inter-layer coupling. All transitions are first-order in the neighborhood of the triple points. We show that the emergent O(4) symmetry of the coexistence state breaks down clearly when the interlayer coupling becomes sufficiently large, but for a weak coupling, of the magnitude expected experimentally, the enlarged symmetry can still be observed at the triple point up to significant length scales. Thus, it is likely that the plaquette-singlet to antiferromagnetic transition in SrCu$_2$(BO$_3$)$_2$ exhibits remnants of emergent O(4) symmetry, which should be observable due to additional weakly gapped Goldstone modes.
Among the quantum many-body models that host anyon excitation and topological orders, quantum dimer models (QDM) provide a unique playground for studying the relation between single-anyon and multi-anyon continuum spectra. However, as the prototypica
l correlated system with local constraints, the generic solution of QDM at different lattice geometry and parameter regimes is still missing due to the lack of controlled methodologies. Here we obtain, via the newly developed sweeping cluster quantum Monte Carlo algorithm, the excitation spectra in different phases of the triangular lattice QDM. Our results reveal the single vison excitations inside the $Z_2$ quantum spin liquid (QSL) and its condensation towards the $sqrt{12}timessqrt{12}$ valence bond solid (VBS), and demonstrate the translational symmetry fractionalization and emergent O(4) symmetry at the QSL-VBS transition. We find the single vison excitations, whose convolution qualitatively reproduces the dimer spectra, are not free but subject to interaction effects throughout the transition. The nature of the VBS with its O(4) order parameters are unearthed in full scope. Our approach opens the avenue for generic solution of the static and dynamic properties of QDMs and has relevance towards the realization and detection of fractional excitations in programmable quantum simulators.
Noethers theorem is one of the fundamental laws of physics, relating continuous symmetries and conserved currents. Here we explore the role of Noethers theorem at the deconfined quantum critical point (DQCP), which is the quantum phase transition bey
ond the Landau-Ginzburg-Wilson paradigm. It was expected that a larger continuous symmetry could emerge at the DQCP, which, if true, should lead to emerged conserved current at low energy. By identifying the emergent current fluctuation in the spin excitation spectra, we can quantitatively study the current-current correlation in large-scale quantum Monte Carlo simulations. Our results reveal the conservation of the emergent current, as signified by the vanishing anomalous dimension of the current operator, and hence provide supporting evidence for the emergent symmetry at the DQCP. Our study demonstrates an elegant yet practical approach to detect emergent symmetry by probing the spin excitations, which could potentially guide the ongoing experimental search for DQCP in quantum magnets.
We study the Neel-paramagnetic quantum phase transition in two-dimensional dimerized $S=1/2$ Heisenberg antiferromagnets using finite-size scaling of quantum Monte Carlo data. We resolve the long standing issue of the role of cubic interactions arisi
ng in the bond-operator representation when the dimer pattern lacks a certain symmetry. We find non-monotonic (monotonic) size dependence in the staggered (columnar) dimerized model, where cubic interactions are (are not) present. We conclude that there is an irrelevant field in the staggered model that is not present in the columnar case, but, at variance with previous claims, it is not the leading irrelevant field. The new exponent is $omega_2 approx 1.25$ and the prefactor of the correction $L^{-omega_2}$ is large and comes with a different sign from that of the formally leading conventional correction with exponent $omega_1 approx 0.78$. Our study highlights the possibility of competing scaling corrections at quantum critical points.
Deconfined quantum critical points govern continuous quantum phase transitions at which fractionalized (deconfined) degrees of freedom emerge. Here we study dynamical signatures of the fractionalized excitations in a quantum magnet (the easy-plane J-
Q model) that realize a deconfined quantum critical point with emergent O(4) symmetry. By means of large-scale quantum Monte Carlo simulations and stochastic analytic continuation of imaginary-time correlation functions, we obtain the dynamic spin structure factors in the $S^{x}$ and $S^{z}$ channels. In both channels, we observe broad continua that originate from the deconfined excitations. We further identify several distinct spectral features of the deconfined quantum critical point, including the lower edge of the continuum and its form factor on moving through the Brillouin Zone. We provide field-theoretical and lattice model calculations that explain the overall shapes of the computed spectra, which highlight the importance of interactions and gauge fluctuations to explaining the spectral-weight distribution. We make further comparisons with the conventional Landau O(2) transition in a different quantum magnet, at which no signatures of fractionalization are observed. The distinctive spectral signatures of the deconfined quantum critical point suggest the feasibility of its experimental detection in neutron scattering and nuclear magnetic resonance experiments.
