We study the spin-excitation spectrum (dynamic structure factor) of the spin-1=2 square-lattice Heisenberg antiferromagnet and an extended model (the J-Q model) including four-spin interactions Q in addition to the Heisenberg exchange J. Using an improved method for stochastic analytic continuation of imaginary-time correlation functions computed with quantum Monte Carlo simulations, we can treat the sharp (δ-function) contribution to the structure factor expected from spin-wave (magnon) excitations, in addition to resolving a continuum above the magnon energy. Spectra for the Heisenberg model are in excellent agreement with recent neutron-scattering experiments on CuðDCOOÞ 2 · 4D 2 O, where a broad spectral-weight continuum at wave vector q ¼ ðπ; 0Þ was interpreted as deconfined spinons, i.e., fractional excitations carrying half of the spin of a magnon. Our results at ðπ; 0Þ show a similar reduction of the magnon weight and a large continuum, while the continuum is much smaller at q ¼ ðπ=2; π=2Þ (as also seen experimentally). We further investigate the reasons for the small magnon weight at ðπ; 0Þ and the nature of the corresponding excitation by studying the evolution of the spectral functions in the J-Q model. Upon turning on the Q interaction, we observe a rapid reduction of the magnon weight to zero, well before the system undergoes a deconfined quantum phase transition into a nonmagnetic spontaneously dimerized state. Based on these results, we reinterpret the picture of deconfined spinons at ðπ; 0Þ in the experiments as nearly deconfined spinons-a precursor to deconfined quantum criticality. To further elucidate the picture of a fragile ðπ; 0Þ-magnon pole in the Heisenberg model and its depletion in the J-Q model, we introduce an effective model of the excitations in which a magnon can split into two spinons that do not separate but fluctuate in and out of the magnon space (in analogy to the resonance between a photon and a particle-hole pair in the exciton-polariton problem). The model can reproduce the reduction of magnon weight and lowered excitation energy at ðπ; 0Þ in the Heisenberg model, as well as the energy maximum and smaller continuum at ðπ=2; π=2Þ. It can also account for the rapid loss of the ðπ; 0Þ magnon with increasing Q and the remarkable persistence of a large magnon pole at q ¼ ðπ=2; π=2Þ even at the deconfined critical point. The fragility of the magnons close to ðπ; 0Þ in the Heisenberg model suggests that various interactions that likely are important in many materials-e.g., longer-range pair exchange, ring exchange, and spin-phonon interactions-may also destroy these magnons and lead to even stronger spinon signatures than in CuðDCOOÞ 2 · 4D 2 O.

We investigate first-and second-order quantum phase transitions of the anisotropic quantum Rabi model, in which the rotating-and counter-rotating terms are allowed to have different coupling strength. The model interpolates between two known limits with distinct universal properties. Through a combination of analytic and numerical approaches we extract the phase diagram, scaling functions, and critical exponents, which allows us to establish that the universality class at finite anisotropy is the same as the isotropic limit. We also reveal other interesting features, including a superradiance-induced freezing of the effective mass and discontinuous scaling functions in the Jaynes-Cummings limit. Our findings are relevant in a variety of systems able to realize strong coupling between light and matter, such as circuit QED setups where a finite anisotropy appears quite naturally.

The ability to store information is of fundamental importance to any computer, be it classical or quantum. To identify systems for quantum memories, which rely, analogously to classical memories, on passive error protection ("self-correction"), is of greatest interest in quantum information science. While systems with topological ground states have been considered to be promising candidates, a large class of them was recently proven unstable against thermal fluctuations. Here, we propose two-dimensional (2D) spin models unaffected by this result. Specifically, we introduce repulsive long-range interactions in the toric code and establish a memory lifetime polynomially increasing with the system size. This remarkable stability is shown to originate directly from the repulsive long-range nature of the interactions. We study the time dynamics of the quantum memory in terms of diffusing anyons and support our analytical results with extensive numerical simulations. Our findings demonstrate that self-correcting quantum memories can exist in 2D at finite temperatures.

We investigate the exact solution of the honeycomb model proposed by Kitaev and derive an explicit formula for the projector onto the physical subspace. The physical states are simply characterized by the parity of the total occupation of the fermionic eigenmodes. We consider a general lattice on a torus and show that the physical fermion parity depends in a nontrivial way on the vortex configuration and the choice of boundary conditions. In the vortex-free case with a constant gauge field we are able to obtain an analytical expression of the parity. For a general configuration of the gauge field the parity can be easily evaluated numerically, which allows the exact diagonalization of large spin models. We consider physically relevant quantities, as in particular the vortex energies, and show that their true value and associated states can be substantially different from the one calculated in the unprojected space, even in the thermodynamic limit